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The electrolytic separation coefficient, a, has been determined for the potassium isotope 
separation in the electrolysis of aqueous solutions of potassium chloride at a flowing mercury 
cathode. Suspension of purified potassium chloride crystals in a liquid mixture of bromoform, 
n-pentanol, and n-hexanol served as the analytical procedure in determining the isotopic 
composition of samples. Densities of the liquid suspension mixture as a function of temperature 
were independently determined. The electrolytic separation factor for the potassium isotopes 
was found to be 1.0054+0.0005 for the temperature range of 15°-50°C, and, within the limits of 
error, was independent of temperature, of the fraction electrolyzed, of the concentration of the 
electrolytic solution, of the current density, and of the amount of back reaction at the cathode. 
An empirical relation between the electrolytic separation factors of elements thus far investi- 


gated and their atomic weights is presented. 





I. INTRODUCTION 


HE electrolytic separation of isotopes other 
than those of hydrogen has received com- 
paratively little attention. Washburn! showed 
that oxygen isotopes were noticeably separated 
by electrolysis and Johnston? determined the 
electrolytic separation factor, a, for O'* relative 
to O'§ to be 1.008+0.003 for the electrolysis of 
aqueous alkaline solutions between iron elec- 
trodes using a densimetric method of analysis. 
Eucken and Bratzler? were the first to in- 
vestigate the electrolytic separation of lithium 
isotopes and were able to set 1.07 as an upper 
limit for @ in the electrolysis of aqueous lithium 
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sulfate at a mercury cathode. Taylor and Urey‘ 
electrolyzed aqueous lithium hydroxide at a 
mercury cathode and obtained an a of 1.020 on 
the basis of mass spectrometer analysis by 
Brewer. Holleck® electrolyzed aqueous lithium 
chloride and obtained an a of 1.079 by a chemical 
atomic weight method of analysis. Taylor and 
Urey® reported an a of 1.039 for the electrolysis 
of aqueous lithium chloride. Employing a densi- 
metric method of determining relative atomic 
weights, Johnston and C. A. Hutchison’ obtained 
an a of 1.055+0.005 for the electrolysis of 
aqueous and alcoholic solutions of lithium 
chloride at a mercury cathode. They were able to 
show, within their limits of observation, that the 


( 037} I. Taylor and H. C. Urey, J. Chem. Phys. 5, 597 
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869 (1940). 


401 





402 DWIGHT A. 
a-value was independent of the cell temperature, 
the fraction electrolyzed, the amount of back 
reaction at the cathode, and within a small range, 
independent of current density. 

Johnston and the writer,’ employing a densi- 
metric method of isotopic analysis, electrolyzed 
aqueous sodium chloride at a platinum anode and 
obtained an a of 1.0060+0.0005 for the separa- 
tion of chlorine isotopes. This work also showed 
no significant dependence of the a-value on the 
cell temperature, the fraction electrolyzed, the 
amount of back reaction at the anode, and the 
current density in the range of observation. The 
a-value on Acheson graphite, 1.0068+0.0014, 
which was computed from the data of Yacouban® 
agrees with our own more precise value. The 
writer’s recent redetermination’® of the a-value 
for chlorine, with similar apparatus but by a 
different procedure, resulted in the value, 1.0063. 
The linear average of a’s obtained from all the 
chlorine electrolyses gives a= 1.0061+0.0005. 

The purpose of the present paper is to report 
the results of an investigation on the electrolytic 
separation of potassium isotopes at a mercury 
cathode and to show the experimental relation 
between the a-values for the elements thus far 
investigated and their atomic masses. 


II, EXPERIMENTAL DESCRIPTION 
A. The Electrolysis 


The apparatus used tocarry out the electrolyses 
was similar to that of Eucken and Bratzler* and 
to that of Johnston and C. A. Hutchison.’? The 
cell consisted of a four-liter round bottom Pyrex 
flask to the bottom of which was sealed an 
extension which supplied a means for a flowing 
mercury cathode. The cell contained a rubber 
stopper at its top through which passed a 
platinum lead to a platinum gauze anode, a 
thermometer for the measurement of the elec- 
trolyte temperature, and glass tubes for passage 
of gases from the cell which were formed during 
the electrolyses. These surfaces extending into 
the cell and the inner wall of the cell were coated 
with “biolloid’’ paraffin. The electrolyte em- 
ployed was an aqueous solution of potassium 
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chloride. The chlorine evolved during the elec- 
trolyses was swept from the cell by a stream of 
nitrogen and absorbed in four traps, two of which 
contained sodium hydroxide and the remainder, 
acetone. The nitrogen stream also served to give 
constant stirring of the electrolyte and yielded 
what was thought to be a uniform mixture of the 
solution throughout an electrolysis. The mercury 
for the flowing cathode was supplied from a 
reservoir and the rate of flow through the cell was 
adjusted at values varying from 2.5 to 10 liters 
per hour for the various electrolyses. The temper- 
ature of the cell was regulated by an external 
water bath. By this means the temperature of the 
electrolyte could be controlled to within +3°C. 
The amalgam from the electrolyses was drained 
from the cell into an eight-liter Erlenmeyer flask 
which was heavily coated with paraffin and 
treated as described below. 


B. Preparation of Reagents 


Potassium chloride employed in the electrolyses 
was prepared by six successive precipitations 
from aqueous solutions by addition of pure 
hydrogen chloride. The purification of po- 
tassium chloride, which was used in the prepara- 
tion of standard crystals of normal isotopic 
composition for use as comparison crystals in the 
densimetric analyses, was effected by a method 
similar to that employed by Richards and Wells” 
in their determination of the atomic weight of 
sodium and chlorine. The beginning material in 
all cases was Baker’s analyzed C.P. potassium 
chloride. Details of the purification procedure 
and evidence of its effectiveness have been 
published elsewhere.’* 

The amalgam formed during electrolysis was 
placed in eight-liter Erlenmeyer flasks which had 
been coated with paraffin. Pure hydrochloric 
acid“ was added in such quantity that its 
equivalent concentration was in slight excess of 
the concentration of the potassium ion to be 
liberated. Platinum foil was added to promote 
reaction of the amalgam. The resulting solution 
was boiled to dryness in a platinum container. 


1 See reference 13 for details of purification. 

12 T. W. Richards and R. C. Wells, “A Revision of the 
Atomic Weights of Sodium and Chlorine,” Carnegie Insti- 
tute of Washington, Publication No. 28 (1905). 

131), A. Hutchison, Phys. Rev. 66, 144 (1944). 
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The residue was given two fractional crystalliza- 
tions from water and then fused on platinum to 
eliminate traces of hydrochloric acid. The salt 
was then given eight more crystallizations from 
aqueous solutions. The material from each step of 
the crystallizations was centrifuged before the 
succeeding step. The resulting material was used 
in the preparation of crystals for suspension.'*!5 

The mercury used for the flowing cathode was 
purified by washing with dilute nitric acid and 
water, drying with filter paper, and finally by two 
successive distillations in vacuum. The purifi- 
cation was repeated before using the mercury 
again. 

Water used in all reagent preparations and 
purification procedures was obtained by double 
distillation from an aqueous alkaline permanga- 
nate solution in a fused quartz still. Various 
water samples tested had a specific conductivity 
of approximately 10-7 ohm. 

The suspension liquid used in the densimetric 
analyses was a mixture of organic liquids similar 
to that used by H. L. Johnston and C. A. 
Hutchison.? The mixture had the following com- 
position: bromoform, 40.00 ml; -hexanol, 16.94 
ml; m-pentanol, 16.00 ml. Baker’s bromoform 
(boiling point range, 150°-151°C) U.S.P. IX, was 
washed with water, dried with calcium chloride, 
and distilled under vacuum. The middle third of 
the distillate was collected in a receiver con- 
taining a small amount of n-hexanol and n- 
pentanol (25 drops of each) which acted as a 
negative catalyst preventing decomposition. The 
composition was then adjusted to that of the 
above suspension liquid by adding drop by drop 
Eastman practical grade alcohols which had been 
double distilled. The density of this suspension 
mixture remained constant within the experi- 
mental limits of error (+5X10-* g/ml) of the 
density determination during the period of 
densimetric analyses. 


C. The Isotopic Analyses and 
Fraction Electrolyzed 


The isotopic analyses were made in an appa- 
ratus similar to that used by Johnston and the 
Writer’ in their work on the electrolytic separa- 


—_— 
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tion of chlorine isotopes and that used by 
Johnston and C. A. Hutchison’ in their work on 
the separation of lithium isotopes. Details of the 
isotopic analyses have previously been outlined 
for the work on lithium.’ In order to get repro- 
ducible suspension temperatures, it was found 
necessary to anneal the crystals at a temperature 
50° below their melting point for a period of four 
hours with gradual cooling over a period of four 
hours and to store them over phosphorus 
pentoxide under vacuum. With these precautions 
suspension temperatures of individual crystals 
were reproducible to better than +0.002° (equiv- 
alent to 3.3X10-* g/ml in the density) except in 
step IV of electrolysis 12. 

The density of the suspension liquid, as a 
function of temperature, was determined by 
careful hydrostatic weighings of a 9.4-ml Pyrex 
glass bob in air, in pure water, and in the sus- 
pension medium. The necessary corrections for 
the buoyancy of the bob in air and the thermal 
expansion of Pyrex glass were applied. The ex- 
perimental data and details of the weighing 
procedure have been described elsewhere.” 

Eighteen determinations of the density of the 
suspension medium over an approximate temper- 
ature range of 25.0° to 29.3°C yielded a good 
straight line with a slope of — 1.883+0.004 x 107% 
g/ml/degree. Taking into account the thermal 
coefficient of expansion for potassium chloride,'* 
namely, a (cubical) = 113 10-*, it is found that 
one degree in the AT of suspension corresponds to 
1.659 X 10-* g/ml in the isothermal density of the 
potassium chloride. This corresponds to 0.0623 
unit in the atomic weight of potassium and to 
3.11 percent in the distribution of the potassium 
isotopes. The determination of the suspension 
temperatures to within +0.002°C, which is some 
better than the limit of reproducibility (actual 
reproducibility is referred to here, not the proba- 
ble error of suspension temperatures) thus corre- 
sponds to 6.2X10-* percent in the isotopic 
abundance of the potassium. 

The fraction of chemical equivalents of po- 
tassium ion electrolyzed was determined by 
analysis of the residual electrolyte. Samples of the 
residual electrolyte were boiled to dryness in a 


16 Klemm, Tilk, and Mullenheim, Zeits, f anorg. allgem, 
Chemie 176, 1 (1928), 
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TABLE I. Data for the separation of K isotopes in the single step electrolyses of aqueous KCL. 








Run number 1 2 3 

Average temperature (°C) 35 35 35 

Current density (amp./cm?) 0.139 0.099 0.060 

Concentration of KCI at be- 0.50 0.50 0.50 
ginning (mols/liter soln.) 

Equivalents of K in starting 1.48817 1.41512 1.38902 
material 

Equivalents of K liberated 0.10178 0.11714 0.12310 

Fraction of equivalents of K 0.071266 0.082777 0.088623 
electrolyzed 

Current efficiency» 0.78 0.75 0.82 

Suspension temperatures of 0.010 0.011 0.009 
different crystal fragments 0.010 0.010 0.009 
relative to standard crystals 0.009 0.010 0.008 
taken as zero* 0.009 0.010 0.008 

0.009 0.009 0.008 
I £,2,2 I 

Average suspension AT 0.0094 0.0100 0.0084 

Cm in the atomic weight —0.000586 —0.000623 —0.000523 
fs) 

a 1.0054 1.0049 1.0049 


4 5 6 7 8 
50 40 30 30 15 
0.099 0.099 0.099 0.099 0.099 
0.50 0.50 0.50 1.50 0.5 
1.50696 1.43252 1.42860 4.14090 1.49052 
0.11486 0.08567 0.11292 0.30930 0.09487 
0.076219 0.059803 0.079042 0.074693 0.063648 
0.77 0.55 0.76 0.69 0.64 
0.009 0.010 0.010 0.010 0.011 
0.009 0.010 0.010 0.009 0.010 
0.009 0.010 0.009 0.009 0.010 
0.009 0.009 0.009 0.009 0.010 
0.008 0.009 0.009 0.008 0.010 
I, } ae I oe 
0.0088 0.0096 0.0094 0.0090 0.0102 
—0.000548 —0.000598 —0.000586 —0.000561 —0.000635 
1.0050 1.0055 1.0055 1.0051 1.0058 








® I indicates imperfect crystals not used to obtain the average AT. 
b See p. 874 of reference 7. 


platinum container. The resulting potassium 
chloride was then fused to remove traces of excess 
chloride ion. Solutions of this potassium chloride 
were titrated with standard silver nitrate solu- 
tion. This furnished a satisfactory method of 
determining the number of residual equivalents 
of potassium ion, the whole procedure requiring 
about fifteen minutes. The number of equivalents 
electrolyzed was obtained by difference of the 
initial and residual number of equivalents of 
potassium ion. 


III. EXPERIMENTAL RESULTS 


Nineteen electrolyses were performed to de- 
termine the electrolytic separation factor, a.!7 
Experimental conditions were varied to deter- 
mine their effect on the a-value. In Table I there 
are presented the data for eight single step 
electrolyses. In the first three the current density 
was varied with the cell temperature held con- 
stant at 35°C. In the fourth through the eighth 
electrolyses, the temperature was varied from 15° 
to 50°C with the current density held constant at 
0.099 amp./cm.? The approximate concentration 


17 ~ is defined by the relationship 
din K® =ad In K* (a) 


where K** and K“ represent the numbers of atoms of the 
isotopes with mass 39 and 41, respectively. 





of the electrolyte at the beginning of each 
electrolysis is recorded for the electrolyses of both 
Tables I and II. In order to increase the isotopic 
separation and thus increase the precision of the 
determination of a, four stepwise electrolyses 
were performed. The data for these electrolyses, 
nine through twelve, are given in Table II. The 
Roman numerals indicate the number of the step 
of a given electrolysis. The potassium chloride, 
obtained on reaction of the K amalgam from a 
given step, served as the starting material for the 
succeeding step. The ninth and tenth electrolyses 
were made in two steps, the eleventh in three 
steps, and the twelfth in four steps. The a’s 
obtained here for steps other than the first are 
cumulative'® a’s. 

With each step of the electrolyses, five crystals 
of potassium chloride were prepared and _ sus- 
pended. If the AT of any of the five crystals 
deviated by a certain amount from the closely 
agreeing AT’s of remaining ones, then the non- 
agreeing AT was discarded on the basis that it 
was from an imperfect crystal. In this case new 
crystals were prepared and suspended until in all 
there were five crystals with AT’s with no large 
differences. Crystals were excluded as imperfect 

18 The cumulative a is that value taken as a constant 


through all steps of an electrolysis, needed to yield the 
experimental AT observed at the end of a given step. 








099 


i) 
49052 


09487 


.063648 


64 
O11 
010 
.010 
010 
010 
Tf FF 
0102 


000635 


0058 


each 
both 
topic 
f the 
lyses 
yses, 
The 
step 
ride, 
ym a 
r the 
lyses 
three 
» a's 
f are 


stals 
sus- 
stals 
osely 
non- 
at it 
new 
n all 
large 
rfect 


stant 


d the 





SEPARATION OF POTASSIUM 


ISOTOPES 


TABLE II. Data for the separation of K isotopes in the stepwise electrolyses of aqueous KCI. 














Run number 91 9II 101 1011 111 1111 11111 121 1211 12111 12IV 

Average temperature 35 35 35 35 35 35 35 35 35 35 35 
(°C) 

Current density 0.099 0.099 0.099 0.099 0.099 0.099 0.099 0.099 0.099 0.099 0.099 
(amp/cm?) 

Concentration of KCI at 2.50 0.50 3.00 0.50 4.00 1.00 0.33 3.94 4.14 1,82 0.61 
beginning (mols /liter 
soln.) 

Equivalents of K in 7.25062 0.42802 9.02098 0.51808 11.8393 2.55165 0.47121 17.7480 3.63492 0.88367 0.25813 
starting material 

Equivalents of K 0.48040 0.09152 0.58162 0.41412 2.4185 0.53201 0.10026 4.3925 1.10236 0.29831 0.10442 
liberated 

Fraction of equivalentsof 0.066256 0.21382 0.064474 0.20060 0.16963 0.20850 0.21277 0.24749 0.30327 0.33758 0.40452 
K electrolyzed 

Current efficiency> 0.62 0.74 0.63 0.56 0.69 0.79 0.83 0.57 0.59 0.91 0.91 

Suspension temperatures 0.010 0.018 0.009 0.019 0.010 0.018 0.027 0.009 0.019 0.025 0.035 
of different crystal frag- 0.009 0.018 0.009 0.018 0.009 0.018 0.026 0.009 0.017 0.025 0.033 
ments relative to stand- 0.009 0.018 0.009 0.018 0.009 0.018 0.026 0.009 0.017 0.024 0.032 
ard crystals taken as 0.009 0.018 0.008 0.018 0.009 0.018 0.025 0.008 0.017 0.024 0.031 
zero* 0.009 0.017 0.008 0.018 0.008 0.017 0.025 0.008 0.016 0.022 0.029 

II I I f,i,2 II I II if ee SF I II 

Average suspension AT’ 0.0092 0.0178 0.0086 0.0182 0.0090 0.0178 0.0258 0.0086 0.0172 0.0240 0.0320 

Change in the atomic —0.000573 —0.001109 —0.005358 —0.001134 —0.000561 —0.001109 —0.001607 —0.000536 —0.001072 —0.001495 —0.001 
weight of K 994 

a 1.0052 1.0053 1.0049 1.0054 1.0055 1.0055 1.0053 1.0055 1.0056 1.0053 1.0054 








* T indicates imperfect crystals not used to obtain the average AT. 
> See p. 874 of reference 7. 


if their AT differed by more than 0.006°C from 
the average AT of the closely agreeing crystals. 
Tables I and II record the AT’s of all crystals 
suspended. Linear averages of the A7’s were 
taken for the computation of a@ since there were 
no large or systematic differences for the AT’s of 
any step. 

a’s for each step were computed from the 
equation, 











1 dMxoAT 7 
log 1 F\1+ 
L pxci(A4 —Ax)) J 
a= (1) 
% €dMxcqAT 7 
log 1-7 i- 
sucktin~Ae) 








in which F is the fraction of chemical equivalents 
of potassium ion electrolyzed in a given step of an 
electrolysis; Mxc: is the molecular weight of 
normal potassium chloride; Ax is the atomic 
weight of normal potassium ; Ag and A, are the 
mass values of Ky and Ka, respectively; pxci is 
the density of normal potassium chloride at the 
suspension temperature for standard normal 
crystals; d is the isothermal density change for 
a AT of one degree ; and AT is the suspension tem- 
perature of crystals prepared from the potassium 


amalgam minus the suspension temperature of 
standard normal crystals. This equation is an 
exact solution of equation (a) of reference 17 in 
terms of the experimentally measured quantities, 
F and AT. Values used for the constants in 
Eq. (1) are: Ajg!®=38.975; A4®=40.974; 
Ax” = 39.096; Mxci?®=74.553; d=1.65910- 
g/ml-deg.; and pxcr* = 1.98651 g/ml. 


IV. DISCUSSION OF RESULTS 


A. The Value of a@ and Its Dependence on 
Experimental Conditions 


The least precisely determined quantity used 
in the computation of a is the AT of suspension. 
Therefore, it is worth while to consider the 
probable error of the AZ’s in some detail. In 
previous work” § the experimental uncertainty of 
the AT’s was given an arbitrary assignment of 
+0.002°C. Suppose, for example, in electrolysis 1 
we assign an uncertainty of +0.002°C and a 
weight of unity to each of the five AT’s obtained. 
The value of the probable error, p;, of the 
average AT obtained from internal consistency is 

19 E. Pollard, Phys. Rev. 57, 1186-7 (1940). 

20 Eleventh report of the committee on atomic weights of 


the International Union of Chemistry, J. Am. Chem. Soc. 
63, 845 (1941). 
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+0.0009°. The probable error, p., obtained from 
external consistency is +0.0002°. The ratio, 
p.-/ pi, is 0.22. However, theory predicts that this 
ratio should not differ from unity by more than 
22 percent. Thus, the arbitrary assignment of 
+0.002° to the A7’s is too large. Let us consider 
the original measurement of a suspension tem- 
perature. It has been obtained in the present and 
previous works by the relation 


Ts=(Tut+T1z)/2 (2) 


where 7's is the suspension temperature, and Ty 
and 7, are the high and low temperatures, re- 
spectively, that define a fall or rise of a crystal in 
the suspension medium.” Relation (2) is about 
the best that can be assumed, especially when the 
difference, (T7—T 1), is of the same order of 
magnitude as the errors in Ty and T,. If it is 
assumed that a temperature can be read with the 
aid of a magnifying glass from the Beckmann 
thermometer, which was employed, to within 
+0.001°, then by taking four measurements of 
both Ty and 7 ;, as was done in the present 
work, the probable error of either Ty or Ty is 
reduced to +0.0005°. From relation (3), there is 
produced a probable error of +0.00035° in Ts. 
The relation used to compute the AT’s of sus- 


pension is, 
AT=T,—Ty, (3) 


where 74 is the suspension temperature of 
crystals prepared from the amalgam and Ty is 
the suspension temperature for standard normal 
crystals. The error produced in each AT of 
Tables I and II by relation (3) is +0.0005°. The 
resulting probable error in the average AT for a 
single electrolysis is +0.00022° by internal con- 
sistency and +0.00017° by external consistency. 
The ratio, p./pi, is then 0.77 in good agreement 
with that predicted by theory when no system- 
atic errors are introduced with the temperature 
measurements and relation (2) is assumed cor- 
rect. Certainly, because of the manner of meas- 
uring AZ’s in this work, the assignment of 
+0.001° as the error in the average AT’s should 
well include all errors. 

The probable error of the a-value for each 
electrolysis was computed on the basis of the 
errors associated with each quantity in the right 


21 See reference 7 for the procedure of suspending crystals. 
p pe g cry 
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member of Eq. (1). The probable errors in these 
a’s were used to assign weights to the nineteen 
a’s obtained. The resulting average a is 1.0054. 
The ratio, p./p:, for the average a is 1.48 where 
theory predicts its deviation from unity should 
be only 10 percent. Hence, systematic errors 
were shown to be present in one or more stages 
of the work, and therefore, the above assignment 
of weights to the a’s was discarded. To obtain an 
average a, weights were arbitrarily assigned to 
each a numerically equal to the number of 
electrolytic steps employed in its determination, 
i.e., the single step a’s were given a weight of one, 
the double step a’s were given a weight of two, 
etc. The resulting average is 1.0054. The error 
limit was taken as the maximum deviation of 
individual a’s from this average which yielded 
+0.0005. 

The a-values for electrolyses 4-8 vary from 
1.0058 to 1.0050 over the temperature interval of 
15°-50°C. If a deviation in the average AT of 
+0.001° is assumed, the corresponding deviation 
in a is +0.0007 unit. It is doubtful if any 
significance can be attached to the above ap- 
parent variation of a with temperature since the 
variation is only slightly larger than that pro- 
duced by an experimental uncertainty of +0.00!° 
in AT. Likewise, no significance can be attached 
to the slight variation of a with change in current 
density recorded for electrolyses 1-3. Within the 
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limits of observation, the a-value is not dependent 
on the amount of back reaction as shown by the 
current efficiencies, on the concentration of the 
electrolyte, or on the fraction of chemical equiva- 
lents of K ion electrolyzed. The value, 


a= 1.0054+0.0005, 


can be stated as the electrolytic separation 
coefficient for the separation of potassium iso- 
topes at a mercury cathode over the temperature 
interval of 15°-50°C. 


B. The Relation Between a and the 
Atomic Weight 


It is of interest to note the relationship which 
becomes apparent when the a-value for the 
separation of potassium isotopes obtained in the 
present work is compared with the a-values given 
in the introduction of this paper for the elements, 
lithium, oxygen, and chlorine. Assuming an 
a-value for hydrogen of 5.000, the a’s for these 
five elements have been plotted against their 
respective atomic weights. There results a good 
smooth curve. The actual value assumed for 
hydrogen, which Topley and Eyring*® found to 
vary from 2.8 to 7.6 depending on the nature of 
the electrode, will not change the shape of the 
curve appreciably for atomic weights greater 
than 5 or 6. Figure 1 presents that part of the 
curve with atomic weights above 6. Since these 
a's have been determined at different electrodes 
(iron, mercury, and bright platinum) a cursory 
conclusion can be made that the mechanism of 
the separation is independent of the electrode 
material for elements with atomic weights 
greater than 5 or 6. The fact that the separation 
of oxygen and chlorine was effected at an anode 
and the separation of lithium and potassium was 
effected at a cathode would lead one to guess that 
the separation mechanism is similar at either 
anode or cathode, at least in its over-all effect. 

The data presented in this work do not serve to 
distinguish among the various mechanisms” that 


® B. Topley and H. Eyring, J. Chem. Phys. 2, 217 (1934). 

*H. Eyring and A. Sherman, J. Chem. Phys. 1, 345 
(1933) ; Horiuti and Polanyi, Nature 132, 931 (1933); R. P. 
Bell, J. Chem. Phys. 2, 164 (1934); B. Topley and H. 
Eyring, J. Chem. Phys. 2, 217 (1934); Fowler, Proc. Roy. 
Soc. Al44, 452 (1934); Bowden and Kenyon, Nature 135, 
105 (1935); O. Halpern and Ph. Gross, J. Chem. Phys. 3, 
452 (1935); Butler, Zeits. f. Elektrochemie 44, 55 (1938). 
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have been proposed to account for the elec- 
trolytic separation of isotopes. Efforts to explain 
the mechanism of separation have involved two 
kinds of factors: (1), kinetic factors determined 
by processes leading to different rates of elec- 
trolysis of the different isotopic species, and (2), 
an equilibrium factor due to isotopic exchange 
between the electrolyzed material and the resid- 
ual electrolyte. Sufficient precise auxiliary data 
are not available to test directly the various 
proposed mechanisms. However, the dependence 
of the observed a-value on temperature may be 
employed to make a rough conclusion as to 
whether a kinetic or an equilibrium factor plays 
the major role in the separation. Neither the 
proposed kinetic mechanisms, which are based on 
different rate determining steps, nor the isotopic 
exchange mechanism possesses sensible tempera- 
ture dependence and such small differences in the 
temperature coefficient of a as are predicted 
could not be detected by the relatively precise 
method of analysis employed in this work. On 
first consideration, it may appear that Bell’s* 
kinetic mechanism, which is a quantum-me- 
chanical application of Gurney’s*® theory of 
overvoltage, refutes the latter statement. Bell’s 
mechanism is based on the supposition that the 
rate limiting step in the separation is the pene- 
tration of a solvent sheath surrounding a 
hydrated ion by an electron. On this basis, he 
predicts that the temperature dependence of the 
In a should be equal to —a/T. This, however, 
appears to be incorrect. It is believed that the 
temperature coefficient of In a should be equal to 
—In a/T. When this correction is made, it is seen 
that the Bell-Gurney mechanism, also, does not 
possess sensible temperature dependence. 

If the observed a represents a partial approach 
to an exchange equilibrium in addition to a 
kinetic mechanism, and if the exchange a is 
assumed to differ markedly from a kinetic a, it 
would be expected that temperature would have 
a noticeable effect on the observed a, by bringing 
about a change in the extent to which equilibrium 
is approached. Since a is constant within the 
error limits of this determination over the tem- 
perature range of 15°-50°C, it is most probable 


4 R. P. Bell, J. Chem. Phys. 2, 164 (1934). 
% R. W. Gurney, Proc. Roy. Soc. (London) A134, 137 
(1932). 
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that a kinetic mechanism effects the major part 
of the separation. 

It is apparent that before any proper explana- 
tion of the mechanism of the electrolytic separa- 
tion process can be forthcoming, there must be a 
much greater extension of the experimental 
material. The writer has work in progress at 
present on the electrolytic separation of nitrogen 
isotopes which is based on the discharge of the 
ammonium ion at a mercury cathode. The experi- 
mental curve of Fig. 1 predicts the a-value for 


GRIFFITH 


nitrogen to be near that of oxygen. It will be of 
interest to check this prediction. Because of the 
much more rapid method of mass spectrometer 
analysis employed in this nitrogen work instead 
of the long and tedious crystal suspension method 
employed in previous work, much greater varia- 
tion in experimental conditions of the electrolyses 
are made available. It is hoped that the results of 
the nitrogen work will furnish experimental data 
that can serve as a more adequate basis for 
interpreting the mechanism of the separation. 
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An X-Ray Diffraction Investigation of the Thermal Decomposition of Silver Oxalate*'} 
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GRIFFITH 


Kodak Research Laboratories, Rochester, New York 


A single crystal of silver oxalate was studied by x-ray diffraction during the process of thermal 
decomposition. The diffraction patterns indicated that some fragmentation of the single crystal 
occurred early in the reaction. Except for distortion, the monoclinic structure was maintained 
until all trace of silver oxalate disappeared. The gradual appearance of polycrystalline silver was 
apparent. The diffraction patterns showed that a portion of the silver assumed preferred 
orientations, the ratio of oriented to unoriented silver remaining constant throughout the 
reaction. Four different preferred orientations were found and were defined relative to the 


silver oxalate lattice. 


ILVER oxalate, when heated, undergoes 

thermal decomposition, the products being 
metallic silver and carbon dioxide.'? The crystal 
structures of silver oxalate* and metallic silver‘ 
are known. It seemed quite probable, therefore, 
that information regarding. the behavior of the 
silver oxalate structure during decomposition, 
the growth of silver, and its crystalline form and 
orientation with respect to the silver oxalate 
structure could be obtained by means of x-ray 
diffraction studies during the process of thermal 
decomposition. Such information is useful in 
connection with the development of photo- 


“* Paper ea eg before Division of Physical and 


Inorganic Chemistry, American Chemical Society, Atlantic 
City, New Jersey, April 8-12, 1946. 

t Communication No. 1068 from the Kodak Research 
Laboratories. 

1 J. Y. Macdonald and C. N. Hinshelwood, J. Chem. Soc. 
127, 2764 (1925). 

2 J. Y. Macdonald, J. Chem. Soc. 832, 839 (1936). 

. L. Griffith, J. Chem. Phys. 11, 499 (1943). 

4R. W. G. Wyckoff, The Structure of Crystals (Chemical 

Catalog Company, New York, 1931), second edition, 
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thermographic processes® and also in connection 
with the theory of photographic development. 
A thick plate-like crystal of silver oxalate, 
prepared by slow evaporation of a concentrated 
ammoniacal solution,?, was mounted on_ the 
goniometer head of a Weissenberg x-ray diffrac- 
tion camera and so oriented as to rotate about 
its a axis. The crystal was then reduced to a 
cylindrical form having a diameter of approxi- 
mately 0.4 mm, by rotating it against a small 
piece of sponge soaked in ammonium hydroxide. 
The goniometer head could be removed from the 
camera without disturbing the attached crystal. 
The experimental procedure consisted in heat- 
ing the crystal, still attached to the goniometer 
head, for a period of time in an oven at a con- 
trolled temperature ranging from 100° to 140°C 
for the various periods of heating. After each 
interval of heating, the crystal was returned to 


’U. S. Patent, 1,976,302, S. E. Sheppard and W. 
Vanselow (assigned to Eastman Kodak Company). 
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THERMAL DECOMPOSITION OF SILVER OXALATE 


TABLE I. Diffraction patterns of silver oxalate during thermal decomposition. 








Photograph Heat treatment 
N Type Time (Hr.) Temp. °C 


Silver oxalate pattern 


Remarks 
Silver pattern 





. 1} 100 
R 3 100 


Normal single crystal. None. 
Normal pattern plus ten weak None. 


spots just off 0 and +1 layer 


lines. 
Normal pattern. Extra spots 


Three weak powder lines. 


slightly stronger. 


Normal pattern slightly diffuse 


Six silver lines. Maxima on lines. 


at large angles. About 50 extra 


: spots. . 
Normal pattern same as in No. 


Pattern considerably stronger. 


6. Extra spots stronger. 


Normal pattern. Spots at large 


One line visible. Maxima on 


angles slightly diffuse. line. 


Normal pattern becoming more 


Same as No. 7. 


diffuse. No additional extra 


spots. 
Same as No. 8. 


About same as No. 9. 
Spots diffuse. 


Normal pattern and extra spots 


Two silver lines with distinct 
maxima. 

Lines much stronger. 

Four lines present. All have 
maxima. 

Intensity greater. 


becoming weaker and more 
diffuse. 


Same as No. 12. 

Same as No. 13. 

Normal pattern and extra spots 
are weaker. 

Same as No. 14. 


Pattern considerably weaker. 


Much weaker. All extra spots 


Intensity greater than No. 12. 

Same as No. 13. 

Stronger. Eight lines. Maxima 
at same points on lines. 

Six lines visible. Same maxima 
present. 

— lines visible. Maxima on 
all. 

Pattern stronger. 


present. 


Weaker. Spots quite diffuse. 
Spots quite diffuse. 
No evidence of pattern. 


No evidence of pattern. 


Same as No. 19. 

Same as No. 18. 

Eight lines. Maxima not as 
sharp as on No. 18. 

Maxima not as sharp as on 














the camera where rotation and usually zero-layer 
Weissenberg photographs were obtained, using 
filtered copper K radiation. By this procedure the 
effect of thermal decomposition on the silver 
oxalate diffraction pattern and the appearance 
and growth of the diffraction pattern of the 
thermally produced silver could be observed. 
The photographs obtained and the observa- 
tions concerning these photographs are sum- 
marized in Table I. Column 1 gives the number 
of the photograph; Column 2 indicates the type 
of photograph, W, indicating Weissenberg and R, 
rotation photograph. Columns 3 and 4 give the 
thermal treatment received since the previous 
photograph. Notes concerning the silver oxalate 
and silver patterns are in Columns 5 and 6. 
Typical photographs are shown in Figs. 1-7, 
these figures representing the photographs listed 


in Table I as 3, 7, 20, 23, 8, 18, and 22, re- 
spectively. 

Figure 1, which is photograph No. 3 of Table I, 
is the diffraction pattern obtained by rotating a 
single crystal of silver oxalate about its @ axis. 
The pattern consists of sharp spots, all dis- 
tributed along straight lines across the pattern, 
the so-called layer-lines. There is no evidence of 
the silver diffraction pattern on this photograph. 
That silver was undoubtedly present was indi- 
cated by the darkening of the crystal, but the 
amount was too small to give a detectable dif- 
fraction pattern. Figure 2, photograph No. 7, 
contains, in addition to the pattern for rotation 
of a single crystal of silver oxalate, spots which 
are not on layer-lines. These spots can be 
accounted for by assuming that the crystal has 
undergone fragmentation in which small portions 
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Fic. 1. Rotation photograph. Heat treatment: 13 hr. at 100°C. 


45) 0 oa 
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Fic. 2. Rotation photograph. Heat treatment: 19 hr. at 100°C. 


of the original single crystal break away and 
assume orientations in which their @ axes are 
not parallel to that of the parent crystal. These 
extra spots first appear on photograph No. 4 
and increase in number and intensity on going 
from No. 4 to No. 7. All rotation photographs 
following No. 7 up to and including No. 20, 
Fig. 3, show the same extra spots at about the 
same intensity relative to that of the normal 


silver oxalate pattern. Beyond photograph No. 
20, as is shown in Fig. 4, there is no evidence of 
the silver oxalate pattern nor of the extra spots. 
Examination of the series of photographs shows 
that as the reaction progresses, the silver oxalate 
pattern and the extra spots become gradually 
weaker and more diffuse. This is caused by 
disappearance of silver oxalate and the distortion 
of the structure of the remaining silver oxalate. 
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THERMAL DECOMPOSITION 


OF SILVER OXALATE 


Fic. 3. Rotation photograph. Heat treatment: 


61 hr. at 100°C 
17} hr. at 110°C 


60 hr. at 125°C 
39 hr. at 130°C 





Fic. 4. Rotation photograph. Heat treatment: 


61 hr. at 100°C 


39 hr. at 130°C 


173 hr. at 110°C 42 hr. at 140°C 
60 hr. at 125°C 


Apparently, the time required for the decom- 
position to occur in the parent crystal and in 
the fragments which split off is about the same. 
The increase in number and intensity of the 
extra spots between photograph No. 4 and 
No. 7 may be caused by the growth and increased 
perfection of the fragments which split off from 


the parent silver oxalate crystal. The extra spots 
do not appear on the Weissenberg photographs 
because of the presence of the layer-line screen 
in the camera when the photographs were taken. 

The first photograph exhibiting any evidence 
of the diffraction pattern of silver is No. 5, on 
which three weak powder lines are visible. These 
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Fic. 5. Weissenberg photographs. Heat treatment: 19 hr. at 100°C. 


lines and those on succeeding photographs are 
caused by face-centered cubic, polycrystalline 
silver. On photograph No. 6, the silver pattern 
is considerably stronger and close examination 
shows that the powder lines are not of uniform 
density over their length but exhibit definite 
density maxima. This can be seen in Fig. 2, 
which is photograph No. 7. Fig. 5, photograph 
No. 8, a Weissenberg photograph, shows one 
powder line for silver. This line exhibits density 
maxima. As the decomposition progresses, the 
lines on the diffraction pattern increase in in- 
tensity and number. Density maxima are present 
on the lines at the same points in all photographs. 
Comparison of Figs. 3 and 4 shows that the 
density maxima are less pronounced on Fig. 4— 
the photograph taken after all the silver oxalate 
pattern had disappeared—than on Fig. 3, which 
was taken before the disappearance of the 


oxalate pattern. This may be caused by the 
orienting influence exerted on the silver crystal- 
lites by the residual silver oxalate lattice. Con- 
tinued heating at 130°C after the silver oxalate 
pattern disappeared did not diminish the maxima 
any further. 

The fact that the lines due to silver are of 
about the same breadth on all patterns, allow- 
ance being made for differences in exposure, is 
consistent with the view that when the growth 
of a small crystal is initiated it rapidly grows to 
final size, the over-all process taking place by 
the initiation and growth of new crystals rather 
than by the simultaneous growth to a large 
size of a small number of crystals. The mass of 
polycrystalline silver left at the completion of 
the reaction was a pseudomorph of the original 
crystal. 

All the maxima on the diffraction lines of the 
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Fic. 6. Weissenberg photograph. Heat treatment: 
61 hr. at 100°C 60 hr. at 125°C 
173 hr. at 110°C 


silver pattern on the rotation photographs can 
be accounted for both as to position and relative 
intensity by assuming that a portion of the 
polycrystalline silver which is formed assumes a 
preferred orientation with respect to the silver 
oxalate lattice such that the [011] direction of 
the si'ver structure is parallel to the a axis of 
silver oxalate. The fact that the intensity maxima 
are arcs rather than spots indicates a range of 
orientations. From the lengths of these arcsone 
may conclude that the orientation of nearly all 
of those crystallites which are oriented lies within 
+4° of that described here. 

Having established that the [011] direction of 
the criented silver structure is parallel to the a 


axis of the silver oxalate structure, we can: 


define completely the relative orientation in space 
of the two structures by determining the angle 


between some direction in the silver oxalate 
lattice which is normal to its a axis, for example, 
the 6 axis, and a direction in the lattice of the 
oriented silver which is normal to its [011] 
direction, for example, the [011] direction. The 
angle between these directions can be determined 
by well-established methods® from the zero-layer 
Weissenberg photographs shown in Figs. 6 and 7. 
From these patterns it was determined that there 
are at least four different preferred orientations 
of silver with respect to the silver oxalate struc- 
ture. In all four cases the [011] direction of 
silver is parallel to the @ axis of silver oxalate. 
The four orientations are then completely defined 
by the angles between the [011] directions and 
the } axis of the silver oxalate structure. These 


6M. J. Buerger, X-Ray Crystallography (John Wiley and 
Sons, New York, 1942). 
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Fic. 7. Weissenberg photograph. Heat treatment: 


61 hr. at 100°C 


39 hr. at 130°C 


173 hr. at 110°C 42 hr. at 140°C 
60 hr. at 125°C 


angles are —30°, 35°, 75°, and 100°, with a 
range of +5°. These four preferred orientations 
of silver account for all the density maxima which 
appear on the diffraction lines of silver in the 
photographs shown in Figs. 6 and 7. 

The preferred orientation of silver with respect 
to the silver oxalate structure probably arises 
from the attempt of the silver crystals to grow 
_ in such a fashion that at the interface between 
silver and silver oxalate there will be as nearly 
as possible a correspondence in direction and 
magnitude of the silver-silver distances. There 
are no close silver-silver distances in these two 
structures which agree in magnitude. The closest 
distances in metallic silver are 2.88A and 4.08A. 
Those in silver oxalate are 3.34A, 3.35A, and 
3.58A. The process of crystal growth must in- 
volve a considerable adjustment of silver-silver 


distances. The oriented silver probably grows 
in such a fashion as to produce a minimum of 
readjustment of distances. The probability that 
a silver crystal will be in a position of preferred 
orientation is about the same as the probability 
that it will have a random orientation, since, 
as is shown by the diffraction patterns, a con- 
siderable portion of the silver is randomly 
oriented. The ratio of oriented to unoriented 
silver remains about the same throughout the 
process of thermal decomposition. Come 
heating of the mass of silver after decomposition 
is complete does not change the ratio of oriented 
to unoriented silver. 


CONCLUSIONS f 


e ef 
In the crystal of silver oxalate which was 
investigated, thermal decomposition is accom- 
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panied by fragmentation of the crystal, the 
fragments produced being slightly disoriented 
with respect to the original crystal. The un- 
decomposed silver oxalate maintains its structure 
until the reaction is complete, except for distor- 
tion produced by removal of carbon dioxide from 
the lattice. The silver formed by the decomposi- 
tion exists in the form of face-centered cubic 
crystals having a final size of approximately 
10-° cm. Some of the crystals of silver assume 
preferred orientations with respect to the silver 
oxalate lattice, the ratio of oriented to un- 
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oriented silver remaining approximately constant 
during and after the reaction. The explanation 
of the particular preferred orientations assumed 
is probably to be found in the structures of the 
two crystals but it is not immediately evident. 
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Following up the ideas of W. Kuhn's skein theory, J. J. Hermans has recently elaborated a 
theory of the properties of dilute solutions of highly polymerized molecules. Hermans describes 
the behavior of linear molecules by the diffusion of their end points. He finds that the contribu- 
tion of a single molecule to the viscosity coefficient as well as to the double refraction of the 
solution is proportional to the square of the degree of polymerization, in agreement with the 
experiments. In the present article we investigate the statistical behavior of the individual links 
of the solved molecules. In this way we arrive at results equivalent with those of Hermans. 
Our method, however, can also be applied to molecules which possess branching points and 


rings. 


1, INTRODUCTION 


¥ connection with Kuhn’s skein theory, 
Hermans! developed a simplified method to 
describe the influence of macromolecules on vis- 
cosity and double refraction. In this method the 
behavior of the molecule is characterized by the 
diffusion of the end points. This can be described 
as if two mass points which exert a fictive elastic 
force upon each other are diffusing in the liquid. 
The equilibrium distribution of the vector 
joining the end points, showing spherical sym- 
metry in a liquid at rest, is deformed under the 
influence of the flow. It is assumed that this flow 
does not alter the elastic force. By means of a 
relatively simple calculation Hermans shows that 
the contributions of a molecule both to the vis- 
cosity and to the birefringence of flow are pro- 


‘J. J. Hermans, Physica 10, 777 (1943). 


portional to the square of the number of ele- 
mentary links, in conformity with experiment 
(Staudinger, Signer). 

In the following we attempt to develop a more 
detailed theory. Our method also enables us to 
give a treatment of macromolecules which possess 
branching points and rings. 


2. THE BEHAVIOR OF A SYSTEM OF MUTUALLY 
CONNECTED PARTICLES IN A FLOW DE- 
RIVED FROM A VELOCITY POTENTIAL 


Let the systematic force which the liquid 
exerts on the ith particle be —{,v;, where v; is 
the velocity of the particle with respect to the 
liquid. This velocity can be written —v,’+v;,”, 
if v,/ represents the velocity of the liquid and v,”’ 
that of the particle. If the streaming is irrota- 
tional and stationary and if its velocity potential 
is V(r), we have v,/=—V¥(r;). Consequently, 
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the force which the liquid exerts on the ith par- 
ticle can be written 


dU 
X= —-——- 5 ¢;4+-X/, (z), (1) 
Ox; 


U=): (ri). (2) 


Here X;’ is a random force of such a nature 
that for a free particle it would maintain the 
equipartition. The influence of the pressure- 
forces which has here been neglected will be con- 
sidered in Section 6. 

The equations of motion are the same as those 
which would apply to our system in a liquid at 
rest if it were subject to a field of force with 
potential energy U. Introducing generalized 
coordinates g', g’---g* and conjugated momenta 
pi; p, and writing 


(y), 


(Di = Axq') (3) 


T= 20x1G"¢' = 30" pxpi, 


for the kinetic energy 7, the molecules will be 
distributed in p, g-space according to Boltzmann: 


W Ik dpidg. = const. e~(7+¥)/6 TT, dp,dq:, (4) 


where U must be regarded as a function of the 
q's. Boltzmann’s constant times the absolute 
temperature has been called 8. 

Average values of functions of g and g can 
now be calculated in principle. A function which 
is linear in the q’s (or p’s) will on the average be 
zero. The average value of a function which is 
quadratic in g can simply be found by first 
reducing it to the average value of a function 
F(q::--qs) of the g’s alone by means of the rela- 
tion 

(Gd!) = a". (5) 


Averaging a function F(q::--qs) of the q’s 
alone amounts, as is well known,” to computing 
the ratio 


(F(q1> + + qo) )w 


= f Fe-U'9(A)+ TJ, dg / f e-Ul(A)* TIndg‘, (6) 


where A is the determinant of the coefficient a;,.: 


A=Det | ay. (7) 


*W. Gibbs, Principles of Statistical Mechanics, Chap. 6. 
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3. PEARL-NECKLACE MODEL OF A 
MACROMOLECULE 


As a model for a macromolecule we consider 
a number of particles with coordinates r1(x;y,2;) 
++:Ij++-Iy, in which every two successive par- 
ticles are connected by a weightless rod of 
length L. Let the mass m and the frictional 
constant ¢ be the same for all these particles. 
The direction cosines of the ith “link,” i.e., of 
the rod which connects the ith particle with the 
(t+1)st particle, will be called a;, Bi, yi (vector 
w;). If x, y, z (r) are the coordinates of the center 
of gravity of our model, one has 


L 
mat Nest ieet 7s -ta;—[N—-i- 1 Jai: 


—[N—i-—2]aiz2°--—an-1), (8) 


N L? 
> Z;? ~ WN } Lv yay + N&, (9) 
1 


_ /uN—») v2 wu 


4 (10) 
\v(N— up) 


Suv 
. vou 


mL? 


T = (11) 
2N 


Low Sur(@yer) +2Nm(e?+y? +2”). 

The simplest model is obtained if it is assumed 
that the two rods which meet in a particle are 
completely free to rotate with respect to each 
other ; a link would then be more or less analogous 
to W. Kuhn’s “statistical chain-element.”’ In 
that case we introduce, for instance, in addition 
to x, y, z, the polar angles 3,, 9, : 


a, =COS 0y=Cy, 
B,=sin 3, COS gy =SyCy, 
Yu.=sin 3, SIN gY,=SySy. 
The total number of coordinates is 2N+1. The 


coefficients a, in the kinetic energy are easily 
derived from (11) using the following formulas: 


(Gyery) = (SySy+CyCyC yr) 8 py +5 pSrC wv Ou Gr 
HAC yS Syrup Pr +S plrS pO pO 

- sont (gu- *’) (12) 
Su» =sin (Pp — ¢y) 


With 2, 3, 4,5 links, these formulas lead to the 
following expressions for the determinant A of the 
coefficients a;: (the angle between the uth and 


the (u+1)st link is denoted by y,): 





MACROMOLECULES IN 


INHOMOGENEOUS FLOW 


A (apart from a numerical factor). 


$1°s97(1 as 2 cos” v1); 


$1°S9"s3"(1 — }[[cos? yi+cos? y ]), 


$17S0°53754?(1 — 3 cos? Yi+cos? p2+cos? W3 ]+¢ cos? y cos? Ws), 
$17S975375 4?55?(1 — tLcos? Yi+cos? po+cos? ¥3+ cos? Ws | 


+ ,'5[cos? y1 cos” ¥3+ cos? ¥; cos? ~s+ cos? Pe cos? Ws }). 


For N—1=5 we have calculated A only for the 
special case that all links are lying in the same 
plane; for N—1=2, 3, 4, however, the formulas 
(13) are of general applicability. We have not 
yet succeeded in deriving a general expression 
for A holding for arbitrary values of N, but its 
form can be surmised. 

In a liquid at rest the probability that 3, and 
y, have values between #, and #,+d#,, and ¢, 
and g,+dg,, respectively, is 


const. A? J], dd,d¢,. (14) 


From the formulas (13) it is seen how the 
orientations of the various links are correlated: 
The correlation is strongest between successive 
links. The probability that they are perpen- 
dicular to each other is somewhat larger than 
that they are parallel or antiparallel; in the case 
of two links these probabilities are in the ratio 
(4/3)§21.15. 

A somewhat less artificial model would be one 
in which a given constant value yp» (for instance 
the supplement of the tetrahedron angle 109.5°) 
is assigned to the angle between two successive 
links. The number of degrees of freedom in that 
case amounts to 3+2+N—2=N-+3, but it is 
now difficult to introduce generalized coordinates 
in a simple manner. If desired, however, one can 
start from the former case, i.e., retain the com- 
plete number of 2N-+1 degrees of freedom, but 
complete this by introducing a potential energy 


U'=3K {(vi—o)? + (Y2—Yo)? 
++++(Ww-2—Wo)*}. 


We would thus introduce explicitly the bending 
vibrations of the valence bonds instead of 
replacing these by rigid bonds. It is clear that 
the incomplete freedom of rotation of a link with 
respect to another link from which it is separated 
by one neighbor could, if desired, also be ac- 
counted for by a suitable potential energy. 


(15) 





4. INTERNAL FRICTION IN A NON-UNIFORM 
FLOW 


In a volume element of sufficiently small size 
the Poiseuille flow which occurs in the experi- 
ments is of the well-known laminar type 


u=xy, v=0, w=0. (16) 


This flow, however, is not irrotational and we 
will rather replace it by 
u= —dW/dx = jKy, 
v= —OWV/dy = $Kx, 
nei (17) 
w= —dW/dz=0, 


W = —SxKxy. 


Superposing a uniform rotation u= xy, v= —}xx 
on this flow, one obtains the laminar flow (16). 
If x is not too large, the effect of this rotation can 
be neglected, as will be shown in Section 6. 
According to the considerations of Section 2, 
the effect of the flow (17) can be described by a 
potential energy U (compare Eq. (2)): 


(18) 


N 
U=—3o« D xii. 
1 


Omitting the coordinates of the center of gravity, 
this becomes 


KL? 
U= ~ ON. om LuvAyBy, 


since it is clear that the coefficients g,, which 
occur in the dependency of }>xy on the a’s and 
B’s are the same as those which occur in the . 
dependency of >>? on the a’s (see Eq. (9)). 
The influence of the flow is described by the 
Boltzmann factor exp (—U/@) which for suf- 


ficiently small x can be replaced by 
eVle=1—U/@. (19) 


If PzzP2y:+- represents the stress tensor in the 
liquid, the coefficient of internal friction 7 is 
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defined by 


Pp ae 
oe Nay | ax 


= Kn. 


We will now calculate the average contribution 
6P.,, of a single macromolecule per unit volume 
to the value of P.,. The contribution 6 per 
molecule to the viscosity will then be given by 


1 
bn =-6P ry. (20) 


K 
Two effects contribute to 6P.,, the motion of the 
particles and the tension in the links. The former 


amounts to 
r= —mM Viedidmy 


because the probability that in unit time a plane 
of unit area perpendicular to the y axis is 
traversed by the ith particle in the direction of 
the negative y axis is —y,;, and the momentum 
transferred in this process is ma;. 

The second contribution results from the 
tensions in the links. If S, is the tension in the 
uth link, this tension contributes an amount 
S,a, to 6P,, if the link intersects the plane men- 
tioned. The probability for the occurrence of 
such an intersection is L8, per molecule; hence 


67 =L Dd u( Sue uB uw 


If the total force on the ith mass point result- 
ing from the tensions is called X;’’, Y;’’, Z;’’, we 
have 

=>; S;%a;, 


Y/=D; SB, 


where the sum extends over all links 7 which 
meet in the ith particle. Usually there will be 
only two of such links, but in a ramified molecule 
—for which this calculation also holds good— 
there can be three or more. Considering now the 
virial-like expression 


Vedi Xi'y( =D: Y's), 


in which we sum over all the particles, and sub- 
stituting (21), each S;a,; will occur twice, 
viz., for those two particles which form the ends 
of the links concerned, but in these two cases the 
a,’s are of opposite sign. Consequently we may 
write 


V=Q. Spau(vi— ys’) = — Diu Suey LBy, 


where y; and y,’ characterize the two ends of the 


(21) 
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uth link, and where a, is the direction cosine of 
the link 7—7’. Now, if X; represents the force 
which the ith particle experiences from the 
liquid, we have according to Eq. (1): 


X = MX; +e 
Ox 


X ff =m;#;— 


S ididers Xi. 
i 


In the average value of >>; X;’’y; in Boltzmann's 
distribution the contributions of the last two 
terms in the third member vanish, and we finally 
obtain 


b= — UL AKXi' Vi) w= — mM Vd Ei)aw 


(25), 


6P ry = 6) +671 = —m DY i((EYidwt (EVs) 


Si (—) . (22) 


The sum of the first two terms in the third 
member represents the derivative with respect 
to time of the average value of —m >0; &:y;, and 
therefore contributes nothing. For the remaining 
expression one can write )>¢{Xi""yiw, if Xi” 
represents the systematic force resulting from 
the streaming liquid. That this expression is equal 
to «én also follows immediately from the con- 
siderations given by J. M. Burgers’ in connection 
with Einstein’s work. Burgers studied the changes 
which forces — on the liquid in a certain 
region (here —X;’"’) will bring about in the 
flow of the liquid outside that region, and 
arrived at a formula which is equivalent to ours. 

With the aid of (18) and (20) it follows from 
(22) that 


= 36 Lily? w=6F Did (eP +92 +2:") wv 


ai Dandie 
~ 6N uv Zuv\ LV) we 


(23) 
To a sufficient approximation we may here take 
the average over the Boltzmann distribution 7 
the liquid at rest (U=0 in Eq. (4)). In the third 
member we have introduced the abbreviated 
notation 


(@,,) = (uv) 


3 J. M. Burgers, “Viscosity report. II,” Proc. Roy. Acad. 


Amsterdam [1] 16, 128 (1938). 
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for the cosine of the angle between the wth and 
the vth link. 

An accurate calculation of the expression (23) 
would require a precise knowledge of the corre- 
lation between the orientations of two different 
links. Now, in the free pearl necklace these cor- 
relations are small for neighboring links and 
negligible for more remote links. Moreover, in 
those cases with a small number of links where 
the calculation was carried out, the correlation 
only depends on the second power of the cosine 
of the angle between the two links, which means 
that ((uv))s4 becomes zero. It would seem that 
this also applies to arbitrary N-values. Since 
((uu))a = 1, we may therefore write 

tL? N-1 cL? N 
” “a 2 Bun =, 2 ol h) 
cL? (N-1)N(N+1) 


~ 6N 6 





or, if N is assumed to be large: 
bn (1/36) ¢L?N?. (24) 


As in the analogous formula derived by Hermans, 
this equation shows 6y to be proportional to the 
square of the number of links, in conformity with 
Staudinger’s viscosity rule. 

If the free rotation of the links is of a nature 
such as that known from organic chemistry, i.e., 
if two successive links always enclose a given 
angle a, we get 


((u))w=1, ((u, w+1))y=Cos a, 
((u, w2))y=cos? a, ((uv))y=cos!4—"la. 


To a sufficient approximation we may then write 
5 : ¢L?N (25) 

=—¢L?N? ctg? -. 

™=36 ae 


If the links are connected by ordinary C—C 
bonds, a is the supplement of the tetrahedron 
angle, cosa amounts to 3, ctg?a@/2 equals 2, 
and the influence on the viscosity becomes, 
therefore, twice as large as in the case of the free 
pearl necklace. 


5. DOUBLE REFRACTION IN NON-UNIFORM FLOW 


The birefringence is described by a dielectric 
tensor €22, €zy=€yz, etc., which is related to the 
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polarization tensor pzz, pz, resulting from the 
influence of a single molecule per unit volume 
according to the formulas 


€zr=€ot4aGprz, €2y=40Gpry, etc. (26) 


Here é') is the refractive index of the liquid 
and G represents the number of macromolecules 
in unit volume. To simplify matters, we assume 
that the contribution of each link to the polariza- 
tion tensor shows cylindrical symmetry round 
the direction of the link and can therefore be 
written 


1 0 
pPr=a 0 0 
0 1 


a?— 4 By QXYr 
+b] By Br—-} Brym J. (27) 


Yaar yrBx Yn sie 3 


The total polarization tensor of the molecule 
becomes 


N-1 


p= Pr. 
1 


To a first approximation the average contribu- 
tion resulting from a@ is independent of x and 
simply amounts to 

0 

1 


The average contribution which results from the 
anisotropy constant 6 will toa first approximation 
be proportional to «x and originate from the term 
— U/T in the expression (19). 

From considerations of rotational invariance 
it is at once obvious that only the xy component 
of p in (27) gives a contribution which differs 
from zero. Its average value per molecule per 
unit volume will be 

N-1 
Dry=b QD) ((arBx)av)av 


1 


(i a 
= DON » &B uv Zur yy) Dav. 


(28) 


The angular brackets enclosing the product 
in the third member of (28) means that we are 
to integrate over a Boltzmann distribution 
exp (—T/6)IIdpdg. Since no q occurs, we can at 
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once average over the q’s, with A! as weight 
factor : 


(Fu= f FANT da / [AT aa, (29) 


It is suitable first to use the following formula 
which follows immediately from rotational in- 
variance and in which the choice of the coordi- 
nate system no longer plays a part: 


3 (a8) (a yB,+ a8») dav 
= (1/15)( {3 (uw) (Av) — 3 (ur) } Dw. 


We thus obtain for the polarization tensor per 
molecule : 


(30) 


Pez _ Puy saa puZNa, Pez = Pyz a 0, Pry mm, 


31) 
C= Fe, gor 3((MH) 00) u— (Cu) 

~ 608 Auv Spr is M Av MV) ws. 

This shows that the polarization tensor is 
deformed as a result of the flow; it has become 
“tri-axial.’’ For electrical oscillations parallel to 
z the polarization has remained Na, as in the 
liquid at rest. In the directions perpendicular to 
z which make angles of 45° and —45° with the 
x and y axis, however, the polarization has re- 
spectively increased and decreased by an amount 
Cb. 

The constant C is proportional to ¢, and there- 
fore to the viscosity of the liquid, and inversely 
proportional to the absolute temperature. Above 
all, however, we are interested in its dependence 
on the number of links in the molecule. To discuss 
this dependence, we note that according to (31) 
everything depends on the correlation between 
the orientations of different links, in as far as 
\, uw, and » are different. For links which are far 
apart, this correlation is practically negligible. 
In the model of completely free rotation of neigh- 
boring links with respect to each other only the 
correlation between two adjacent links differs 
perceptibly from zero. The main terms are then 
doubtless those for which yp and » in (31) are 
exactly or almost equal. In this case the average 
value of 3(Au)(Av) practically cancels that of (uv) 
if X differs perceptibly from yu, and the sum in 
formula (31) becomes approximately 
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while for the constant C we obtain the approxi- 
mate expression 


(32) 


If in (28) we had only retained the term with 
\=u=v, we would have arrived at the same 
result. Thus the double refraction by a single 
molecule becomes proportional to the square of 
the number of links, as was also found by 
Hermans. It is to be expected that in a more 
accurate calculation applied to the pearl necklace 
with completely free rotation, the result will have 
to be corrected by a factor which does not differ 
much from unity. With the partially free rotation 
of actual molecules the analogous correction 
factor will be larger than 1; for C—C bonds a 
preliminary estimate yielded a value between 
2 and 3. 


6. REMARKS ON THE APPROXIMATIONS 
INVOLVED 


I. Influence of the Pressure Forces 


In formula (1) the force which results from the 
non-uniform pressure in the flow of liquid has 
been neglected. If az; represents the x component 
of the acceleration a at the center of the ith 
particle and if m,’ is the mass of liquid replaced 
by the particle, the x component of the force 
concerned is m,’az;. Since in a stationary flow 
we have 


10 
a,=—- —(VW)?, 


(y), (2 
2 Ax » @) 


the influence of the pressure would result in an 
additional term 


—3 Lim) (V¥;)? 


in U. If V is considered small of the first order, 
our approximation amounts to omitting a term 
of the second order. With the potential YW = — }xxy 
used by us, the ratio between the pressure force 
ma and the frictional force ¢v’’ is of the order 
m’x/, which for a sphere of radius 10-7 in water 
(n=0.01) is about 10~'*x, i.e., negligible. 


Il. Influence of Rotation 


So far we have not accounted for the rotation 
of the liquid. If the flow of the liquid consists 








(32) 
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exclusively of a uniform rotation, the average 
behavior of our system can still be described on 
the basis of an assembly in temperature equi- 
librium if we insert the distribution 


W II dpdg=exp [—(T—(QD)/@]]] dpdg. (33) 


Here Q is the angular velocity of the liquid and 
I the angular momentum of the system. On 
account of the fact that in the absence of forces 
I would be an additive integral of motion, such 
a distribution would be stationary in an ideal 
gas whose particles are systems of the type con- 
sidered here. Maxwell already pointed out, that 
with free particles this distribution corresponds 
to a uniform rotation. of a gas with angular 
velocity 2. Similarly, in the present case, where 
each particle possesses in addition internal de- 
grees of freedom, the distribution mentioned in- 
cludes a superposed uniform rotation of each of 
the mass-points constituting the systems. This is 
at once evident from the identity 


T- (QI) = x. 3m (“?2+y?+2,) a > m(Q(rt;]) 
= Li 5m (t:— {Qr;]})?- 3m,((Qr; ])? ’ 
since it shows that the distribution (33) is 
equivalent to the distribution 
W =e-(7' +0018. 
if T’ is considered as the kinetic energy referred 


to the coordinate system rotating with angular 
velocity 2, while 


U'= —4} Di m,(LQr; })? 


represents the potential of a centrifugal force 
acting on all the particles. The influence of the 
pressure forces can be accounted for by the 
potential 


ad =} >: m;'((Qr; })?. 


If the system consists of spheres with the same 
specific weight as that of the liquid, m;=m,’, and 
U’ and U” exactly cancel. 

If now the system is placed in a liquid which 
itself rotates with angular velocity Q, the sys- 
tematic and the random forces exerted by the 
liquid will be the same in the rotating coordinate 
system as in a liquid at rest: the distribution 
(33) will, therefore, be left undisturbed. 

If the flow is stationary but of an arbitrary 
nature and not irrotational, the equilibrium dis- 


tribution can no longer be described by one of 
the simple assemblies known in the kinetic 
theory of gases. Using a method known from the 
theory of the Brownian movement,‘ it is then 
still possible to construct an equation of dif- 
fusion which describes the change with time of 
an arbitrary distribution in phase space. We 
would then have to determine the stationary 
solution of this equation. In the special case of 
the flow (16), representing the superposition of 
a pure rotation and an irrotational flow, both of 
which are proportional to x, the mathematical 
problem becomes more simple. Yet it would not 
be permissible to assume that the stationary dis- 
tribution is now described by a simple com- 
bination of (33) and (4), i.e., by 


W=exp [—(T+3«I.—3«f 2 xy)/0]. 


This is seen, for instance, from the fact that with 
increasing x the system will gradually be oriented 
parallel to the direction of the flow,’ whereas 
according to this formula there can only be 
question of a preference orientation at 45°. Each 
of the two terms following 7 in the exponent 
applies only in the absence of the other one. 
However, if their influence on the distribution 
is small, errors will only occur if one would be 
interested in effects of the order x’. Effects of the 
order « will be described by the distribution 


W =e-7T!9{1—3«I,/O0+3K¢ > xy/6}. 


As the term with J, cannot have any influence 
on the double refraction (since it corresponds to 
a uniform rotation of the liquid) we are led to 
conclude that it was actually permissible to 
neglect the rotation in the preceding section. 
The rotation can also be neglected in the vis- 
cosity. The terms which must be added to the 
sum >}. X”’y in Section 4 in order to account for 
the rotation in the determination of 6Pz,, 
exactly cancel in as far as quantities ~« are 
concerned. This follows at once from the fact, 
that 6P,,=0 in pure rotation, but it can, of 
course, also be verified by direct calculation (the 
term with.—3«J,/6 in the last mentioned formula 
will then also play a part). 

A simple estimate shows that even with 
molecules of a length 10-* cm, it will not be 
simple to realize « values of such a magnitude 


4H. A. Kramers, Physica 7, 287 (1940). 
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that «? effects would be perceptible in dilute 
solutions. 


7. RAMIFIED MOLECULES 


In virtue of the derivations of the formulas 
(23) and (31) for viscosity and birefringence, 
these formulas also hold good for branched 
molecules; we will only have to insert other 
expressions for g,,. These g,,’s occur, for instance, 
in the expression ; 


mL? 7 
die 3 Mw Luv(hr), (0) — (dyudy + BuBv +7 u¥r), 


representing the kinetic energy of the molecule 
if the movement of the center of gravity is 
omitted. The value of g,, can now be determined 
in the following manner. 

Cut the molecule at the wth link and at the 
vth link; it will then fall apart in three pieces, to 
wit two end pieces which contain, respectively, 
half the uth and half the vth link, and a middle 
piece. The number of particles in these three 
parts is called Ni, N2, and N3, respectively ; their 
sum is N. With p=» we get N3=0. 

If uv, the expression 

m 2 

T wv =— { Zune) +28 y( ar) +2,,( 0%) }, 

2N 
will be equal to the kinetic energy of a macro- 
molecule, in which only the links u and v possess 
non-zero angular velocities @, and @, while all 
other links perform translations only, in such a 
way that the center of gravity remains fixed in 
space. 7,, therefore also represents the kinetic 
energy of a system of three particles with masses 
mN,, mN2, and mN;, of which both 1 and 3 and 
3 and 2 are connected by a rod of length L. Since 
the g’s are numerical coefficients, it will suffice 
to consider the case where the particles momen- 
tarily are lying on a straight line. Introducing 
the momentary velocities v1, v2, and v3: 


Ni + Nov2+ Nyv3=0, 


v3—0,=L|a,| P Ve—v3=L|a,|, 


we get 
T = 3m(Nyw?+ Nove? + N33") 


mL? 


= {Ni(N3+Noe) (44) +2NiN2(i0) 


2N 
+(Nit+N3)N2(vr)}. (34) 
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This shows that the value of g,, is equal to the 
product of the numbers of particles in the end 
pieces which are formed when cutting the th 
and the vth link: 


Luv=NiNo. (35) 


Specializing 4 =», or also directly from (34), we 
obtain 


Zunp=MiNe (36) 


where NV; and N2 are the numbers of mass-points 
in the two pieces in which the molecule falls 
apart when cutting the uth link. 

We will now restrict ourselves to the viscosity. 
In order to compute the friction from (23), we 
may assume that with completely free rotation 
of two adjacent links the average value ((uv)),, 
for two different links (uv) is almost always 
zero, as it was in the unbranched pearl necklace. 
For the effect of a single molecule on the friction 
this gives 

cL? tL? 


cL? 
6n =—— =— (37) 
n 6N i Sup 6 


(NiN2)m, 


where the angular brackets with subscript »4 now 
indicate the average value over all possibilities 
of cutting. 

For an unbranched molecule this average value 
becomes 


N-1 


(NiN2)w=— N\(N—WN)) 


4 


1 
auf x(1—x)dx = N’/6 
0 


in conformity with the formula (24). For a 
ramified molecule (NiN2) will always be smaller. 
For instance, if the molecule consists of s 
branches of N/s links each, all starting from a 
single point : 


(NiN2)u = {1/(2s) — 1/(3s?) } N°. 


For s=1 and for s=2 this is equal to N?/6, as 
it should be, but for s>2 it is always smaller 
than N?/6. 


8. VISCOSITY CONTRIBUTION OF A MACRO- 
MOLECULE CLOSED TO A RING 


Let a linear molecule consisting of N particles 
be closed to a ring by joining the Nth mass point 
to the 1st particle by means of an Nth link of 
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length L. The kinetic energy will still be given by 


ma? = 


gs AG a Lur( hv), 


38 
ON ms (38) 


but as a result of the new kinematic bond the 
possible values of 7 will be independent of the 
link which is assumed to have closed the ring. 
The coefficients g,, may, therefore, be replaced 
by the average values of all g,,’s for which the 
uth and the vth link are separated by the same 
number of links in the ring: 
1 N-1 


— N Lp Su+p, v-+p 


=— || 


er 


|u—v|)*}. 
The link between N and 1 will be equivalent 
to the others; the numbers y and » will vary from 
1 to N. In the sum the numbers n+p (and v+ p) 
must be considered as identical with 4+p—N 
(and v-+p—J) in all cases where u+p (or v+p) 
would be larger than N. The third member in 
the equation is obtained by realizing that the 
sum must be extended over all possible positions 
of the closing link, so that it consists of two sums 
of the type }\u(N—yw)=N*/6 (compare Eq. 
(24)). If the distance between yu and » is de- 
scribed by an angle 2nr=22|y—v|/N, which 
would represent the angular distance between 
these links if the ring were a circle, one can also 
write 
Gw25N?2(1—37+377), 


O<r<1. (39) 


The identity of 
_ mL? 


3 DL Gurl av) (40) 


with (38) results, of course, from the relations 


N 
DX (a%) =0 
v=1 
which in their turn are a result of the ring- 
condition 


s wo, = 0. (41) 


1 
From (41) it also follows that >°,,, (4%) =0, 
and we may therefore subtract a constant from 
§.». In other words, instead of (39) one might 
also write 


Ju = —3N?r(1—7). (42) 
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The influence of a ring-shaped molecule on the 
friction will be given by 


cL LL? N-1 


cL? 
eal GE 2 Suv (uv) ww =—— x Guv( (uv) wv. 


43 
6N “7 6N ad 


To calculate the average values ((uv))w, one 
would have to know the distribution of the o’s. 
This distribution is determined in principle by 
the square root of the determinant A of the 
quadratic expression 7(q:---qg;) after trans- 
formation into suitable generalized coordinates. 
It is easy to see that even our former A}, in 
which #1¢1,:--3v-igw—-1 served as generalized 
coordinates, will represent the distribution con- 
cerned, if we restrict ourselves to those values of 
0 and ¢ for which the distance between the Nth 
and the ist particle is LZ. In fact, it would have 
been a sufficient approximation if the kinematic 
relation which expresses this condition is omitted 
and replaced by the introduction of a potential 
energy U’=}K(rw:°—L’)*. Here ry: represents 
the distance between the 1st and the Nth par- 
ticle, while K is a constant of sufficient mag- 
nitude. The distribution of the w’s would then 
be described accurately by A} exp (— U’/6), and 
this amounts exactly to the prescription given. 

However, without going into detailed calcu- 
lations, the sum in (43) can be determined for a 
pearl necklace with completely free rotation if 
we take into account that the correlation be- 
tween two different links u and v which leads to 
a non-zero ((uv))w-value, will almost always be 
the same, independent of whether u and » are 
now lying close together or not. It thus follows 
from the ring-condition (41) that 


0=D.( (ur) dav = ((uH) wt Dov (Cur) dv 


pv 


=1+(N—1)((uv))w ux, (44) 


(uv) dav eS 
Consequently, if we use the expression (42) for 
Juv: 


1 
Le Gur((ur) w= Dow m= ze Juv 


N ux 


1 1 
=0+— wf 3N?r(1—7)dNr 
N 0 


=3N*(5—3) = iy. my 
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Since the analogous sum was N*/6 for the free | 


linear pearl necklace, it appears that the ring 
formation has reduced the influence on the friction 
to 4 of the original amount. 

For a linear pearl necklace with fixed angle a 
between two successive links >>, ((uv))w was equal 
to ctg? 3a, and the sum needed only to be taken 
over a small number of neighboring links (unless 
a would be too close to zero). It is clear that in 
the case of rings, the formula (44) must then be 
replaced by 


1 
(uv) daw anv—— cts? 2, 


and that the sum (45) must also be multiplied by 
a factor ctg? $a. Thus, here again, as was to be 
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expected, it still holds good that ring closing 
reduces the effect to 3 of the original figure. 

There will now be no further difficulty in the 
problem of arbitrarily ramified molecules con- 
taining one or more separate rings. As soon, 
however, as the different rings are no longer 
separate, i.e., if they have one or more chain 
segments in common, the ((uv))\’s can no longer 
be determined in this simple manner. 


Editor’s Note: 

This article is a translation by J. J. Hermans and H. A. 
Kramers of an original appearing in Physica in January, 
1944. Because of wartime restrictions imposed by the Germans, 
the original appeared in Dutch. The innovation of publishing 
a translation was thought to be justified by the inaccessibility 
of the original to many readers because of the language barrier. 
The editors and publishers of Physica gave their kind per- 
mission for the publication of this translation. 
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Conductometric investigations were made of electrolytic 
solutions in xylene and of a hydrocarbon oil both of which 
have conductivities of from 10-” to 10-" mho/cm at room 
temperature. The electrolytes introduced were composed 
of an organic basic and an organic acidic material, the 
former being an amine. A three-component system (amine, 
aliphatic acid, and a phenol) was generally required in 
order to obtain high conductivities at relatively low concen- 
trations of electrolytes. It is shown that the effect of the 
phenol is only partly caused by the increase of dielectric 
constant; the chief process is probably a specific one. On the 


1. INTRODUCTION 


¢ certain previous publications'~* the author 
has introduced e.m.f. measurements in the 
study of insulating liquids for the purpose of 
obtaining additional information on the nature 
of electrolytic ions in such systems. Up to that 
time, the measurement of the electrical conduc- 
tivity was essentially the only tool used for this 
particular purpose and it seemed that the addi- 
1 Andrew Gemant, J. Chem. Phys. 10, 723 (1942). 


* Andrew Gemant, J. Chem. Phys. 12, 79 (1944). 
* Andrew Gemant, J. Chem, Phys. 13, 146 (1945). 





theory that the mechanism consists in the formation of an 
addition-compound that subsequently dissociates into ions, 
the experimental data were subjected to an analysis. It 
appears that generally five or more single molecules com- 
bine to form a compound and that the concentration of the 
compound is a small percentage of the total solute present. 
The results have a bearing on the problem of deterioration 
of insulating oils. The high conductivities of the latter, 
even in the presence of certain colloidal constituents, are 
ultimately caused by electrolytic ions of the type in- 
vestigated. 


tion of the e.m.f. method, which is most useful 
in the study of aqueous and alcoholic soiutions, 
would be of value in this field, too. Results ob- 
tained in this series of researches to date seem 
to justify this enterprise. 

The author used the glass electrode for the 
measurement of hydrogen ion concentrations in 
hydrocarbon solvents and hydrocarbon oils. 
E.m.f. data obtained from oils to which a strong 
acid (organic sulfonic acid) was added could be 
interpreted quantitatively on the basis of a 
reversible electrode action of the glass. When, 
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however, the technique was used for oils that 
were oxidized in the laboratory, or deteriorated 
in service, the hydrogen ion concentration ap- 
peared much too low to account for the meas- 
ured electrical conductivity. The conclusion was 
drawn, therefore, that the cation in such oils was 
not the hydrogen ion. It should be pointed out 
that this conclusion was based on direct proof 
obtained from specific e.m.f. data. This result is 
in keeping with those obtained by J. D. Piper 
and co-workers‘; these authors showed that 
solutions of purified acids in paraffin oil do not 
produce the high conductivity present in oxi- 
dized oils. 

Since the cations in oxidized or deteriorated 
oils, which exhibit increased conductivity, are 
not hydrogen ions, the further conclusion was 
drawn that they must in all probability be 
organic in nature, just as are the anions. What 
their exact chemical nature is could probably be 
found from e.m.f. measurements, if appropriate 
electrodes were available. In attempting to do 
further work in this direction, the next step in 
the research seemed to be the establishment and 
conductometric investigation of solutions in hy- 
drocarbon oils of electrolytes of which both ions 
are organic. Obviously, organic amines are best 
suited for such a purpose, although there are 
other compounds (certain triphenylcarbinol de- 
rivatives, for instance) that under certain circum- 
stances, particularly in combination with halo- 
gens, show cation-like behavior. The conclusion 
that the cations in hydrocarbon solvents are 
amines is supported in the literature. Fuoss 
and Kraus® investigated the conductivities of 
tetraisoamylammonium nitrate in mixtures of 
dioxane and water. J. B. Dickey® has obtained a 
series of patents on treating yarns with certain 
compounds (diethylcyclohexylamine salts of 
phosphinated ketones and the like) in order that 
static electricity be avoided during subsequent 
weaving. Since the yarns during weaving are in 
oil, it is likely that the organic salts are dissolved 
in the oil and form semi-conducting solutions. 
J. D. Piper and collaborators! also have ex- 


‘J.D. Pi r, A. G. Fleiger, C. C. Smith, and N. A. 
Kerstein, Ind. Eng. Chem. 32, 1510 (1940) ; 34, 1505 (1942). 
*R. M. Fuoss and C. Kraus, J. Am. Chem. Soc. 55, 2387 


rig 
°j. Dickey, U. S. Pat. 2,256,112; 2,286,791, and 
2,286 194, 1941 and 1942. 


amined, among a great variety of added ma- 
terials, the effect of salts of organic amines upon 
the electrical conductivity of paraffin oil; they 
obtained high conductivities at 80°C. 


2. SOLUTIONS IN XYLENE 


Although the majority of data in this paper 
refer to a hydrocarbon oil, several sets of meas- 
urements were carried out with xylene and 
benzene as solvents. As will be shown in the 
section dealing with hydrocarbon oil, appreciable 
conductivities at relatively low concentrations of 
the electrolytes are obtainable only in a three- 
component system, consisting of an amine, an 
aliphatic acid, and a phenolic compound. In 
solvents such as xylene or benzene, however, 
probably because of their higher dielectric con- 
stant, two-component systems (amine and ali- 
phatic acid or amine and phenolic) are sufficient 
to cause reasonably high conductivities, although 
even in these solvents the three compounds to- 
gether are far more effective. In this section 
solutions in xylene of a-naphthylamine and 
phenol are dealt with in detail. All chemicals 
used in this research were high grade commercial 
products, and were used without additional 
purification. 

The tests consisted in preparing a series of 
solutions of various normalities with regard to 
the solutes added, and measuring the conductivi- 
ties at room temperature in a cell designed by 
J. C. Balsbaugh.? A variable-frequency power 
factor bridge was used, and the range of fre- 
quencies extended from 50 to 20,000 c.p.s. As 
explained later in this paper in the part dealing 
with solutions in oil, no variation of the con- 
ductivity with frequency was observed during the 
many checks at several frequencies. In addition 
to the conductivity, the dielectric constant, too, 
was determined for each solution of the series. 

The obtained conductivities were corrected for 
the conductivities observed when either of the 
two components was added separately. This was 
necessary because the ion formation from elec- 
trolytes consisting of organic cations and anions 
was the only point of interest in this connection. 
Thus, the conductivities caused by either com- 
ponent alone were deducted from the total. If the 


7J. C. Balsbaugh and A. G., Assaf, Ind. Eng. Chem. 13, 
515 (1941). 
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Fic. 1. Conductivity at room temperature vs. concentra- 
tion of solutions in xylene. The conductivities in this as in 
all other figures are corrected for initial conductivities of 
the starting solutions, as explained in the text. 

Curve 1 =Concentration of a-naphthylamine: varying; that of phenol: 
0.23 normal. 

Curve 2 =Concentration of a-naphthylamine: 0.23 normal; that of 
phenol: varying. 

Curve 3 =Concentration of a-naphth'ylamine: varying; that of phenol: 
0.16 normal. 

Curve 4=Concentration of a-naphthylamine: 0.16 normal; that of 
phenol: varying. 

Curve 5 =Concentration of a-naphthylamine: varying; that of phenol: 
0.08 normal. 

Curve 6 =Concentration of a-naphthylamine: 0.08 normal; that of 
phenol: varying. 


hypothesis that ion formation from such types of 
electrolytes plays an important role in liquid 
dielectrics is correct, then the corrections applied 
should be relatively small parts of the total. This 
expectation was very strikingly fulfilled in the 
case of oil as solvent, as will be shown later. For 
xylene the deducted amounts were larger, ranging 
from about 20 to 35 percent of the total. 

In each series the concentration of one com- 
ponent was kept constant while that of the other 
was increased by steps. In three sets the concen- 
tration of a-naphthylamine was kept at, re- 
spectively, 0.08, 0.16, and 0.23 normal; in three 
other sets the concentration of phenol was kept 
at the same percentages. The results obtained 
from these six series is shown in Fig. 1 in which 
the conductivity in mho/cm is plotted against 
the normality of the component that was varied 
by steps. The conductivities are all corrected as 
described above; they are of the order of 10-™ to 
10 mho/cm. It was not deemed necessary to 
present the uncorrected conductivities, as ob- 
served, for all series of measurements presented 
in this paper; as an illustration, however, the 


GEMANT 














CONDUCTIVITY, MHO PER CM 




















002 0.04 
CONCENTRATION, MOLE PER LITER 


Fic. 2. Conductivity vs. concentration of solutions in 
xylene; conductivity data reduced to a dielectric constant 
of 2.45. 


Curves 1 to 6: the same as in Fig. 1. 


uncorrected data will be given for three particular 
series of solutions dealt with in Section 3a. 

One of the reasons for the increase in con- 
ductivity with increasing concentration is the 
increase of the dielectric constant. In order to 
find out in what manner the conductivity de- 
pends on the concentration of dissociable material, 
the effect of the varying dielectric constant must 
be eliminated. Although the variation of the 
latter was rather small in this series (from 2.39 to 
2.48), such a procedure was considered necessary. 
It was carried out in the following manner. 

One of the six series (with 0.08” a-naphthy!- 
amine) was repeated with 3.5 percent by volume 
chloroform added to the solutions. It was as- 
sumed that the essential role of the chloroform 
was to raise the dielectric constant without 
acting specifically upon the various solutes. By 
this means the dielectric constant was increased 
throughout by an increment close to 0.06. The 
conductivities at corresponding concentrations of 
phenol were found to differ by a factor close to 
1.60. In the narrow range of dielectric constants 
involved it is safe to take a simple exponential 
variation between the conductivity and the 
dielectric constant. In other sections of this paper 
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a different procedure was employed, in which no 
use of the special assumption of an exponential 
variation was made. In order to eliminate the 
role of the dielectric constant, all conductivities 
were reduced to an arbitrarily selected value of 
dielectric constant. The value 2.45, lying in the 
middle of the range, was selected. The conduc- 
tivity o at a dielectric constant D is then reduced 
to o2.45, its value at a dielectric constant 2.45, by 
the relation: 


02.45 = 01.60@-45—P)/0.06, (1) 


The corrected conductivity data are presented 
in Fig. 2; the curves, as can be seen, are some- 
what different from those in Fig. 1. The variation 
of « shown in Fig. 2 is caused solely by the 
increase in concentration of the two solutes. 

A possible mechanism of ion generation as 
observed in these systems is the formation of an 
addition compound with subsequent ionic dis- 
sociation (this subject is further discussed at the 
end of this paper). While it is not maintained 
that the experimental data furnish compelling 
proof for this explanation, it is thought to be at 
least a useful guiding hypothesis for the interpre- 
tation of the data. It is possible that a continua- 
tion of these studies will indicate the desirability 
of a different approach to the problem; one of the 
various other possibilities is, for instance, the 
theory of multiple ion formation, as suggested by 
Fuoss and Kraus.® 

Let us then denote the basic constituent by B, 
the acidic one by A, and assume that ” molecules 
of B and m of A form a compound B,A»: 


nB+mA2B,An. 


If concentrations in mole/liter are denoted by ¢ 
with respective indices, we have 


Kate "CA "= CB Ams (2) 


TaB_e I. Conductivities of solutions in hydrocarbon oil. 








Concentration 
in mole/liter 

Conductivity 
in mho/cm 


1.7 10-4 
— J 2.9 

— 0.05 ; 2.1 
0.05 0.05 F 90 
0.205 — ‘ 36 
0.205 0.038 , 1700 


Tributylamine Oleic acid o-Cresol 


0.05 0.05 
0.05 











where K,=association constant. If B+ denotes a 
positive and A~ a negative ion, formed by 
dissociation of the addition compound (this is the 
simplest, although by no means the only possible 
kind of dissociation), the reaction is 


B,Am@Bt+A- 
and, with K,z=dissociation constant: 


Kae, m — CB+CA-. (3) 
Since 
CBp+>=Ca-—=Ci, (4) 


where ¢c;=ionic concentration: 
KaKacp"ca" =C?. (5) 
Since the conductivity o is given by 
=10-*Agc;, (6) 
with Ag=equivalent conductivity, Eq. (5) yields 


o=kcep"?c4””, . (7) 
with 
k=10-*Ao(KKa)}. (8) 


In applying Eq. (7) and (8) to the present case, 
an inspection of Fig. 2 reveals that curves 1, 3, 
and 5 (variable amine content) follow approxi- 
mately a square root function, while curves 2, 4, 
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Fic. 3. Conductivity vs. concentration of solutions in 
xylene; calculated from Eq. (7). 


Curves 1 to 6: the same as in Figs. 1 and 2. 
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TABLE II. Dielectric constants and conductivities (uncorrected) for three series of solutions in hydrocarbon oil. 








Starting solution: 
0.05n oleic acid 
0.05 o-cresol 


Conduc- 
ity of tivity in 

tributyl- Dielectric 10-4 
amine constant mho/cm 


Normal- 


Starting solution: 
0.05” tributylamine 
0.05” o-cresol 


Dielectric 
constant 


Starting solution: 
0.05” tributylamine 
0.05 oleic acid 


Conduc- 
tivity in 
10-14 


Conduc- 
tivity in 
107-14 
mho/cm 


Normal- 
ity of 
o-cresol 


Dielectric 


constant mho/cm 





0.0000 2.210 
0.0096 2.245 
0.019 2.281 
0.028 2.301 
0.036 2.310 
0.044 2.317 
0.052 2.319 
0.060 2.324 
0.067 2.326 





2.232 
2.266 
2.289 
2.303 
2.312 
2.323 
2.331 


0.0000 
0.011 
0.022 
0.032 
0.042 
0.053 
0.062 
0.072 


2.234 1 
2.253 7 
2.275 17 
2.290 32.5 
2.310 60 
2.325 93 
2.341 135 
2.355 187 


75 
1 








and 6 (variable phenol content) are approxi- 
mately linear. It appears, then, that the main 
addition compound consists of one molecule of 
amine (m=1) and two molecules of phenol 
(m= 2). The obvious disagreement of some of the 
curves (clearly of curve 2) with this simple be- 
havior indicates that probably several addition 
compounds are formed, but that the one just 
mentioned prevails. 

The conductivities at concentrations corre- 
sponding to the set of curves in Fig. 2 were 
calculated from Eq. (7) with n=1 and m=2. The 
value of k was adjusted so as to obtain an order of 
magnitude of 10-" for o at c=0.06. The value 
used was 2.210-'°. The plot of these calculated 
data is shown in Fig. 3 which is reasonably 
similar to Fig. 2. 

In further verification of the mechanism here 
depicted, the measurements of Lucasse* may be 
quoted. The melting points of mixtures of p- 
toluidine and acetic acid have a maximum of 
49°C at 66 percent of acetic acid, corresponding 
to a compound: 1 amine, 2 acid molecules— 
exactly the same ratio as the conductivity curves 
in Fig. 2 indicate. Thus Eq. (7) reproduces the 
measured data with one arbitrary constant k in a 
satisfactory manner. 

The value of this constant, as required to 
satisfy the experimental data, may be used for an 
estimate of K,, the association constant. In 
Eq. (8) the order of magnitude of Ao is 60 
(mobility of organic ions in a solvent of viscosity 
of 1 centipoise is around 30) and that of Ka is 


8 W. W. Lucasse, R. P. Koob, and J. G. Miller, J. Phys. 
Chem. 48, 85 (1944). 








10-8; this order of magnitude was obtained? for a 
number of strong electrolytes (acids) in xylene 
and dioxane. With these figures a value of 10 is 
obtained for K, from Eq. (8). This order of 
magnitude is significant: It shows that the com- 
pound formation is not complete (which is the 
case if K, is much larger than 1), nor is the 
amount of compound formed only a very small 
fraction of the components (which is the case if 
K, is much smaller than 1). From Eq. (2) the 
concentration of the compound can be computed, 
using K,=10; if, for instance, both compenents 
are present in a normality of 0.05, the concentra- 
tion of the compound is 0.001 which is about 2 
percent of the amine. 


3. SOLUTIONS IN HYDROCARBON OIL 


The hydrocarbon oil used was a typical trans- 
former oil of a viscosity of 0.14 poise at 25°C, a 
dielectric constant 2.18, and an electrical con- 
ductivity 110-4 mho/cm at 25°C. Three 
different systems of solutes were investigated 
which will be treated separately in the following. 


(a). Tributylamine—Oleic Acid— 
o-Cresol Solutions 


In testing hydrocarbon oils, it was observed 
that addition of an acidic and a basic component 
did not as a rule produce a noticeably increased 
conductivity at room temperature. The author 
was first inclined to believe that, concerning 
hydrocarbon oils, his conclusions with regard to 
the importance of organic cations were incorrect. 
Neither aliphatic acids, nor phenolic compounds 
—both weak acids in aqueous solutions—pro- 
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duced an effect, when added separately to solu- 
tions of the amine. When, however, acids of both 
types were added simultaneously to an oil con- 
taining an amine, the effect became pronounced. 
Table I shows that none of the three possible 
combinations in pairs of the three materials 
added produced a conductivity noticeably higher 
than that of the oil alone, while the three ma- 
terials together had a rather pronounced effect. 
The last two lines show that even if the amine 
and cresol were present in excessively high 
concentrations, producing a conductivity of 
medium level, addition of relatively small 
amounts of acid increased the conductivity about 
fifty times. The correction to be applied to values 
of ¢, in order to obtain the conductivity resulting 
from the combination alone, is only 2 to 3 percent. 

Three series of measurements were carried out, 
each starting with a combination of two of the 
three solutes in concentrations of 0.05 mole/liter 
each. To these solutions the third solute was 
added in increasing concentrations. The dielectric 
constants and the total conductivities, as meas- 
ured, for these three series are shown in Table II. 

The dielectric constant D of the starting solu- 
tions was about 2.23, increasing with increasing 
concentration to about 2.34. This increase, while 
rather small, is significant because of the generally 
low level of D. The process of reducing the 
g-values to one selected value of D was rather 
important. It was done by repeating all three 
series with an oil containing 10 percent by volume 
of chloroform. By this means the above quoted 
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DIELECTRIC CONSTANT 


Fic, 4, Logarithm of conductivity ratio of tributylamine 
—oleic acid—o-cresol solutions in hydrocarbon oil vs. 
difference in dielectric constants. 

Corresponding figures on curve refer to identical solutions except 10 


percent by volume of chloroform added to increase the dielectric 
constant. X refers to curve 1 of Fig. 5. O refers to curve 3 of Fig. 5. 
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range of dielectric constants was raised to a level 
of about 2.40 to 2.50. By examining correspond- 
ing pairs of data (with and without CHCls;), 
having conductivities « and o’ and dielectric 
constants D and D’, various values of ¢/o’ were 
correlated with values of (D—D?’). Or, since a 
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Fic. 5. Conductivity vs. concentration of solutions in 
hydrocarbon oil. 


Curve 1: Tributylamine varying; oleic acid 0.05”; o-cresol 0.057. 
Curve 2: Tributylamine 0.05”; oleic acid varying; o-cresol 0.05”. 
Curve 3: Tributylamine 0.057; oleic acid 0.05; o-cresol varying. 
(Conductivity data reduced to a dielectric constant of 2.35.) 
xX © O and solid lines: experimental data. 
Dotted lines: calculated from Eq. (10). 
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Fic. 6. The same as the solid lines of Fig. 5, plotted on a 
bilogarithmic scale. 
Slope of curve 1 =0.38. 


Slope of curve 2 =0.38. 
Slope of curve 3 =1.35. 
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Fic. 7. Logarithm of conductivity ratio of tributylamine- 
myristic acid—a-naphthol solutions in oil vs. difference in 
dielectric constants. 

X refers to curve 1 of Fig. 8. 


+ refers to curve 2 of Fig. 8. 
O refers to curve 3 of Fig. 8. 


plot of log (¢/o’) vs. D was attempted, values of 
(log «—log o’) were correlated with values of 
(D-—D’). 

These correlated data were then plotted in the 
manner shown in Fig. 4: differences in ordinates 
were plotted against differences in abscissas, 
pairs of data being shown by identical figures on 
the curve. As can be seen, the absolute value of 
the conductivity does not enter into the plot, and 
the curve is assumed to hold independently of the 
latter. With the aid of this curve all data were 
reduced to D=2.35, an average value for the 
range covered. Curve 1 of Fig. 5 refers to varying 
amine concentration, curves 2 and 3 to varying 
acid and o-cresol concentrations, respectively. 
The conductivities increase to about 2X10-” 
mho/cm, an order of magnitude often found in 
deteriorated oils. The values are lower than for 
corresponding concentrations in xylene partly 
because of the higher viscosity of the solvent. 

The same guiding hypothesis that was used in 
the previous case might be used here, too. Apart 
from the 2 molecules of base B, and m molecules 
of acid A, / molecules of a phenolic compound P 
will combine into a compound B,A»,P1: 


nB+mA+/P2B,AmP1 


that subsequently supplies a positive and a 
negative ion. Instead of Eq. (5) we then have: 


KiKacr"ca"cp' =c;, (9) 
and instead of Eq. (7): 


o=kep”!*c4""cp"”. (10) 
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Inspecting Fig. 5, curves 1 and 2 are found to 
follow approximately a square root function, 
while curve 3 is a 14 power function. The same 
curves plotted bilogarithmically in Fig. 6 reveal 
straight lines with slopes of 0.38, 0.38, and 1.35, 
respectively. The dotted lines in Fig. 5 were 
calculated from Eq. (10) according to these 
powers and made to cross at c=0.05 as they 
should if Eq. (10) were strictly correct. Since 
curves 1 and 2 do not coincide, Eq. (10) can be 
only approximately valid. 

With n=1, m=1, and /=3, it appears that one 
molecule of amine, one molecule of acid, and 
three molecules of phenol combine into one: thus, 
a fairly large molecule results. The possible 
structure of such an addition compound will be 
discussed in the last section (discussion) of this 
paper. 

The constant k corresponding to the calculated 
dotted lines in Fig. 5 has a value of 1.910-°. 
This permits K, to be estimated from Eq. (8). In 
the oil of viscosity 0.14 poise, it can be estimated? 
that Ao is around 4 and Kg is of the order of 
magnitude of 210-'®. This.gives K,=1 X10". 
Hence, with all three components present in a 
concentration of 0.05 mole/liter each, the concen- 
tration of the compound is 4X10~‘, about 1 
percent of the amine present. 


(b). Tributylamine—Mpyristic Acid— 
e-Naphthol Solutions 


As mentioned in Section 2, the chemicals used 
in this research were not subjected to laboratory 
purification. This implies that a repetition of 
these tests with purified materials would yield 
numerically somewhat different results. In par- 
ticular, the objection could be raised against the 
data presented under (a) that oleic acid and 
o-cresol always contain admixtures and that the 
effect was, at least in part, caused by them. The 
low conductivity found in any pair of the three 


TABLE III. Conductivities of solutions in a hydrocarbon oil. 








Concentration in mole /liter ee 
Conductivity in 





Tributylamine Myristic acid a-Naphthol mho/cm 
0.05 0.05 — 1.4x 10" 
0.05 — 0.05 7.9 

— 0.05 0.05 5.5 
0.05 0.05 0.05 400 
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materials, however, can be brought against such 
objection. 

Toclarify this point further, a combination was 
chosen in which oleic acid was replaced by 
myristic acid and o-cresol by a-naphthol, both 
crystalline materials. The program for this part 
of the work was the same as that for the previous 
part. Three sets of solutions were investigated ; 
each set contained two of the materials each 0.05 
normal and the third material in increasing 
amounts. The ‘zero values’ are shown in 
Table II]. The conductivity of the solution con- 
taining all three components is 4 X 10-”, more than 
four times that of the previous combination. This 
finding is reassuring in that impurities have proba- 
bly only a secondary role in the effect observed. 

The effect of the dielectric constant was 
evaluated in the same manner as before; the 
result of this ‘‘plot of differences” is shown in 
Fig. 7. Both this plot and Fig. 4 show that the 
effect of the dielectric constant increases with 
decreasing value of the latter, in keeping with the 
theory of Bjerrum.* 
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Fic. 8. Conductivity vs. concentration of solutions in 
hydrocarbon oil. 














Curve 1: Tributylamine varying; myristic acid: 0.05”; a-naphthol: 
.OSn. 
Curve 2: Tributylamine 0.05”; myristic acid: varying; a-naphthol: 


.05n. 
Curve 3: Tributylamine 0.05”; myristic acid: 0.052; a-naphthol: 
_ varying. 
(Conductivity data reduced to a dielectric constant of 2.35.) 
+ O and soiid lines: experimental data. 
Dotted lines: calculated from Eq. (10). 
————— 


*D. A. MacInnes, The Principles of Electrochemistry 


peinhold Publishing Corporation, New York, 1939), p. 


The conductivities, reduced to D=2.35, are 
shown in Fig. 8. As the dotted curves (3 and 

+ power) show, the former is not in agreement 
with the experimental data below 0.01 normality. 
It appears that the type of the reaction changes 
with increasing concentration of the components. 
The experimental data are, however, not con- 
sidered sufficiently complete to permit further 
interpretation of this upward trend of the curves 
at the origin. Using the location of the dotted 
lines, an evaluation analogous to the previous one 
gives k=7.5X10-*; K,=2 X10‘, and the concen- 
tration of the compound for all three components 
present in 0.05 molar concentration is 6X10-%, 
i.e., twelve percent of the amine present. The 
large conductivity effect obtainable with this 
combination is caused, according to the mecha- 
nism assumed, by an increased tendency for 
addition compound formation. The compound 
molecule containing three molecules of a-naphthol 
must have appreciable size. 


(c). Triethanolamine—Oleic Acid Solutions 


In addition to the series discussed, some tests 
were carried out with an amine containing 
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Fic. 9. Conductivity vs. concentration of solutions in 
hydrocarbon oil, 


Curve 1: Triethanolamine: varying; oleic acid: 0.06n 
Curve 2: Triethanolamine: varying; oleic acid: 0.032. 
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hydroxyl groups: triethanolamine. Because of the 
low solubility of this compound only low concen- 
trations were investigated, giving clear solutions 
—at least to the naked eye. The concentration 
limit at which inhomogeneity appeared was 
rather sharp, and measurements were carried out 
only well below that limit. 

The dielectric constant was remarkably con- 
stant in these measurements (from 2.20 to 2.22), 
so that no correction in this respect was neces- 
sary. The results from two series are presented in 
Fig. 9. Curve 1 refers to an oleic acid content of 
0.06 normal; curve 2, to one 0.03 normal. The 
concentrations of the amine are shown on the 
abscissa; the range covered was below 0.0015 
normal. Solutions up to an amine concentration 
of about 0.01 normal were clear. 

Because of the low solubility of this amine (in 
fact, it is insoluble in oil alone and dissolves only 
in the presence of the acid) the existence of a 
second phase, contributing by a Maxwell-Wagner 
effect to the high conductivity, might be sus- 
pected. In this case o should be a function of 
frequency, and in particular, if measured with 
d.c. voltages, distinctly less. While a check of o 
using several frequencies was occasionally carried 
out in all the series reported, a few data relating 
to this particular combination are given in 
Table IV. Even the second of the listed solutions, 
having the high conductivity of 10-!° at room 
temperature, shows no variation in the range 
between a frequency of 20,000 and d.c. voltages; 
this solution, by the way, remains clear to the 
naked eye even after prolonged standing. Ac- 
cordingly, a Maxwell-Wagner effect on inter- 
phases can be ruled out. 

The data obtained in this combination are not 
deemed sufficiently complete for the drawing of 


TABLE IV. Conductivities at various frequencies of 
solutions in hydrocarbon oil. 














Concentration in mole/liter Frequency, Conductivity in 
Triethanolamine Oleic acid c.p.s. mho/cm 

0.00082 0.06 10,000 4.8xX10-" 
0.00082 0.06 1,000 4.6 
0.00082 0.06 100 4.4 
0.00082 0.06 50 4.6 
0.0068 0.06 20,000 96 X10-” 
0.0068 0.06 10,000 96 
0.0068 0.06 1,000 96 
0.0068 0.06 d.c. 100 
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quantitative conclusions. In attempting to apply 
the mechanism suggested (compound formation 
with subsequent dissociation), the steepness of 
the curves suggests that a relatively large number 
of molecules are involved in the compound forma- 
tion. An estimate of about eight molecules of 
amine and four molecules of acid is in keeping 
with the trend of the curves; the relatively large 
size of the compound is in agreement with the 
observed low solubility of the amine. 


4. DISCUSSION OF THE RESULTS 


While, for the sake of a systematic presenta- 
tion, part of the discussion, particularly a quanti- 
tative interpretation of the data, has been in- 
cluded in the corresponding sections dealing with 
experimental results, two subjects deserve special 
consideration. One of these is the mechanism of 
conductivity in the three-compound systems in 
oil. The other subject is the bearing of the 
presented results on the problem of deterioration 
of insulating oils. 


(a). Mechanism of Conductivity in the Three- 
Component Systems in Oil 


If the conductivity is produced in a two- 
component, base-acid system, there can be little 
doubt that the mechanism is compound forma- 
tion with subsequent ionic dissociation. Conse- 
quently the equations presented in Section 2 of 
this paper are probably safely founded. The 
extension of this mechanism to the three com- 
ponent base-acid-phenolic system, however, is 
not quite as obvious. One might first think of 
another possibility, namely, that the phenol 
simply plays the role of an addition to the 
solvent, increasing the dielectric constant of the 
latter. Hence, its action would be merely electro- 
static. This possibility cannot be strictly dis- 
proved, but tests carried out with chloroform to 
show the effect of the dielectric constant upon the 
conductivity rule against it. Table V gives a few 
figures that are instructive in this respect. Solu- 
tions containing the phenolic component have a ¢ 
over 10-”; those without the phenolic have, in 
spite of a considerable amount of CHCl;, a ¢ 
below 10-”. The first solution has a o larger than 
the second, although its dielectric constant is 
lower than that of the second; the same is true if 
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DROCARBON SOLUTIONS 


f solutions in hydrocarbon oil. 











Concentration, mole/liter 


Tributyl- ; 
amine Oleic acid Myristic acid o-cresol 


: Dielectric 
Chloroform constant 


Conductivity 
mho/cm 


a-naphthol 





0.072 


0,05 0.05 
0.05 0.05 
0.05 0.05 
0.05 

0.05 - ~ 
0.05 - - 


0.072 
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the fourth and fifth solutions are compared. 
It can-be concluded that, although the purely 
electrostatic effect of the components is not 
negligible, the phenols have an effect on o 
markedly greater than the electrostatic one. 

Assuming, therefore, a compound formation 
between the three components, a question arises 
regarding the chemical structure of the com- 
pound. The amount of information contained in 
this paper is clearly not sufficient to answer this 
question in a definite manner. However, some 
indication as to the correct answer can be found 
in using the concept of resonance.’° 

The unshared pair of electrons of the amino 
nitrogen attracts a proton; this is the reason for 
the formation of an addition compound between 
an amine and an acid. The same holds for an 
addition between an amine and a phenol or a 
naphthol. If, now, an addition between an acid 
and a phenol should take place that leaves both 
hydroxyls free (hence, not an ester formation), 
then such a compound would be considerably 
stabilized through resonance in the presence of an 
amine, even if such a compound would be 
unstable in the absence of an amine. As an 
illustration of such resonating structures, formu- 
las | and II are given; the position of the electron 
pair of the nitrogen is different in the two cases. 
If more than one molecule of the phenol or 
naphthol is attached to the acid, the stability of 
the compound is further increased. The experi- 
mental data presented indicate this to be the 
case. The addition of the phenol to the acid is 
probably due to secondary valencies; this struc- 
ture obviously needs further elucidation. It 
should be pointed out that there are theoretical 
considerations that might be raised against the 


G. W. Wheland, The Theory of Resonance (John Wiley 
and Sons, Inc., New York, 1944), p. 167 et seq. 


likelihood of such resonating structures. 
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Once a compound of the type suggested is 
formed—and, if it is formed, its concentration 
should be reasonably large, as shown in the 
previous sections—its dissociation to a small 
extent into ions is understandable. Even if the 
aliphatic acid itself in the oil is extremely weak, 
the presence of the phenol groups considerably 
strengthens its acidity..° The reason for this 
effect is, partly at least, the —CO* dipole-moment 
of the phenol. More than one phenol group in the 
compound would further enhance this effect. The 
result of the dissociation is then the positive 
HN*(R); ion with the rest of the compound as 
the negative ion. 


(b). Bearing of the Results Presented on the 
Problem of Deterioration of Oils 


In investigating ionic solutions in oils, the 
question arises, to what extent are the ionic 
components responsible for the high conduc- 
tivities of deteriorated oils. It is certain that such 
oils quite often reveal the presence of colloidal 
materials that appear to be related to their high 
conductivity. This standpoint is particularly 
emphasized by J. D. Piper and his co-workers." 

In the following it is shown that, even if 
colloidal particles add to the conductivity of oils, 


4 Joseph Sticher and J. D. Piper, Ind. Eng. Chem. 33, 
1567 (1941), 
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the ultimate carriers of the electricity are in all 
probability molecularly dispersed ions. Three 
cases are now considered. 

(1) If the inhomogeneity takes the form of a 
microscopic second phase dispersed in the oil, a 
Maxwell-Wagner effect might contribute to the 
conductivity. This, as the theory of the effect 
reveals, is ultimately caused by the bulk con- 
ductivity of the two phases, hence, by electrolytic 
ions. 

(2) If the inhomogeneity is a colloidally dis- 
persed second phase, a so-called lyophobic colloid, 
the particles might, owing to their charge, exhibit 
cataphoresis, thus contributing to the total con- 
ductivity.” Judging from aquecus solutions, the 
charges upon such particles are caused by prefer- 
ential adsorption of one kind of the electrolytic 
ions present in the bulk of the liquid. 

(3) If the inhomogeneity is a lyophilic colloid, 
in the majority of cases a so-called macromo- 
lecular solution, then an electrolytic dissociation 
of such macromolecules contributes to the con- 
ductivity. Examples are proteins, soaps, and the 
particles in Oden’s sulphur sol. In these cases 
the concepts of the collodial particle and of 
the macromolecular electrolytic ion merge into 
one. 

It is not maintained that the solutions investi- 
gated are identical with deteriorated oils. The 


2 Andrew Gemant, J. Phys. Chem. 43, 743 (1939). 
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exact nature of the ionic components can certainly 
vary within considerable limits. 

Concerning the presence of acids and aromatic 
alcohols, these are probably formed during the 
oxidation of the oil. Concerning the basic com- 
ponent, the average total nitrogen content of oils 
appears to be about 0.03 percent, i.e., 0.02 molar. 
Assuming this nitrogen to be largely basic in 
character, this concentration is exactly in the 
range that was found to produce conductivities 
of the order of magnitude 10-” to 10—" mho/cm 
(the same as in deteriorated oils). In particular, 
if traces of hydroxyamines are formed during 
oxidation, the effectiveness of the nitrogen is 
considerably increased, as has been shown. It is 
generally agreed that the nitrogen in oils is 
present in the form of cyclic compounds, and it is 
of interest to note in this connection that substi- 
tution of the cyclic piperidine for tributylamine 
produced conductivities comparable to those 
presented in Section 3. Nevertheless, the presence 
of other types of organic cations, such as oxonium 
ions and triarylcarbinols, has also to be considered, 
although they are much weaker bases than the 
amines. 

A more accurate answer as to the chemical 
nature of organic ions present in oxidized oils 
could be given, if specifically acting electrodes for 
e.m.f. measurements could be developed. The 
author is at present engaged in a research to find 
out whether this is possible. 
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Accommodation coefficients toward bright platinum for 
air, oxygen, nitric oxide, carbon monoxide, carbon dioxide, 
nitrous oxide, methane, hydrogen, and helium are com- 
puted from the (1/K.) versus (1/P) slopes of the thermal 
conductivity measurements of Taylor and Johnston and of 
Johnston and Grilly. The results—which cover the temper- 
ature range 80°K to 380°K—are tabulated for the several 
gases, and are collected in a graph where comparison with 
results obtained by other investigators are also included. 
Our results with hydrogen rise toward a coefficient of unity 
with lowering temperature, in confirmation of measure- 


LTHOUGH the literature contains compara- 

tively little data on the accommodation 
coefficients of gases below room temperature, it 
is generally assumed that these tend to approach 
unity as T approaches zero. In connection with 
measurements of the thermal conductivities of 
gases, recently completed in this laboratory” * 
for the temperature range of 80° to 380°K, by 
the hot wire method, it was possible to evaluate 
accommodation coefficients on the cell wire of 
bright platinum. 


RELATION BETWEEN ACCOMMODATION COEFFI- 
CIENT AND RECIPROCAL PLOT 


The calculation of accommodation coefficients, 
from the thermal conductivity data, is based on 
the experimental slopes of the reciprocal plots of 
“apparent” conductivity versus pressure, for 
various pressures used in the cell. The “‘apparent”’ 
conductivity, Ka, of the gas is obtained, as out- 
lined in the initial paper? of this series, by apply- 
ing to the experimental (uncorrected) conduc- 
tivity, K,, all the necessary corrections except 
that for the “temperature-jump” effect. The 
latter is the effect of the temperature discon- 

* Present address: American Petroleum Institute Re- 
search Project 44, National Bureau of Standards, Wash- 
ington, D. C. 

‘Fellow of the Charles A. Coffin Foundation 1939-40; 
E. I. du Pont Fellow 1940-41. 


*W. J. Tayl % a ' 
219 (1946). or and H. L. Johnston, J. Chem. Phys. 14, 


sean L. Jolnston and E. R. Grilly, J. Chem. Phys. 14, 233 


ments by others, but those for other gases drop steeply at 
temperatures between room temperature and that of liquid 
air. Our data are in generally good agreement with the 
recent accurate work of Amdur, Jones, and Pearlman, at 
room temperature. The accommodation coefficients which 
are derived from the reciprocal K versus reciprocal P plots 
differ slightly from Knudsen accommodation coefficients, 
and are referred to as ‘‘Maxwell-Knudsen”’ coefficients. 
Equations are derived which give the relationship between 
the reciprocal slopes, the temperature jump constants g’, 
and the M-K accommodation coefficients. 


tinuity which exists at a gas-solid interface when 
heat flows across the interface.‘ The temperature 
in the body of a gas will vary in a way deter- 
mined by the boundary conditions of the system. 
Thus the temperature varies linearly with dis- 
tance for heat flow between parallel plates, and 
linearly with the logarithm of the radial distance 
for heat flow between concentric cylinders (as in 
a hot wire cell). However the particular way in 
which the gas molecules incident upon the surface 
are reflected from it (or re-emitted after adsorp- 
tion), and, further, the way in which the reflected 
and incident molecules interact through colli- 
sions, affects the condition of the gas within a 
few mean-free-path lengths of the surface. As a 
result the temperature in the gas near the surface 
will deviate from the functional relation which 
obtains in the body of the gas, and the extrapola- 
tion of the latter relation to the surface will yield 
a temperature different than the temperature of 
the surface. This temperature discontinuity will 
be denoted by 67.5 

Experiment has established, and the theory, 
in so far as it has been developed, supports the 
fact that 67 is proportional to the temperature 
gradient in the gas at the surface (as extrapolated 
from the body of the gas), and is inversely pro- 


4E. H. Kennard, Kinetic Theory of Gases (McGraw-Hill 
Book Company, Inc., New York, 1938), Sections 176-177. 

5 The actual temperature of the gas at the surface will 
also differ from the temperature of the surface, although by 
an amount less than 67°. 
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portional to the pressure of the gas (or directly 


oT oT 
6 =g—=(¢'/p)—, (1) 
or or 
where g’ is independent of the pressure, p. The 
constant g has been called the temperature jump 
distance.‘ It is suggested that g’ be called the 
temperature jump coefficient. 

In the hot wire cell a temperature discon- 
tinuity is found at the surface of the cell wire, 
and to a lesser extent at the surface of the cell 
wall. The discontinuity which enters into the 
derivation of the equations for the present cell 
is 67, as defined here. As a result of these dis- 
continuities, the apparent conductivity of the 
gas is a function of the pressure of the gas (as 
well as of other constants of the gas and of the 
apparatus). It has been shown previously? that 
a plot of the reciprocal of the apparent con- 
ductivity, Aq, versus the reciprocal of the pres- 
sure, yields a straight line, and that the intercept 
at (1/p)=0 is the reciprocal of the true con- 
ductivity, K, or 


1/Ka=1/K+A/P. (2) 


The slope A is related to the coefficient g’ of Eq. 
(1) by the expression? 


d(1/Ka) 1 
~ d(1/p) In (r2/r1) 
XL (g1’/ Kiri) +(g2"/Kore) |, (3) 


where the subscripts 1 and 2 refer to the surfaces, 
or temperatures, of the cell wire and the cell wall, 
respectively. The term (ge’/Kore) in Eq. (3) 
makes only a small contribution, since (g2’/Ke) 
~(gi'/K,), and (1/re) =0.032(1/r1). It is there- 
fore permissible to set (go’/Kz2) equal to (g1'/K1), 
whereupon Eq. (3) may be solved for gi’, 


gi’ =In (r2/r1)-L(1/r1) +(1/re) + KA. (A) 


If in Eqs. (1) to (4), r is expressed in centimeters, 
and p in centimeters of mercury, then g’ has the 
units centimeters times centimeters of mercury. 
On introducing the dimensions of the present 
cell, Eq. (4) becomes, 


gi’ =0.0426-K,A. (4a) 





The constant g,’ is dependent upon the nature of 
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proportional to the mean-free-path length). Thus 
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the gas and of the surface, but is independent of 
the dimensions of the apparatus® or of the 
pressure of the gas. It should be emphasized that 
the validity of Eq. (1) is confirmed by the fact 
that the observed data are precisely in accord 
with Eq. (2), and that the values of gi’ calculated 
from Eq. (4) are not dependent upon any detailed 
theory of the interaction of the molecules of the 
gas with the surface. The usefulness of g’ lics 
in the fact that it may be used with Eq. (1) to 
calculate 57 for the same gas and surface in any 
other situation in which they occur. In such a 
case the derivative in Eq. (1) is merely the 
derivative normal to the surface, whether the 
surface is curved or flat. 

The accommodation coefficient of the gas 
molecules on the surface has been calculated from 
the value of g’ and a theoretical relation con- 
necting these two constants. 

Knudsen’s accommodation coefficient, for heat 
conduction across a gas-solid interface, is defined 
as the ratio‘ 


a=(E;—E,)/(E:—E,’), (5) 


when &; is the energy carried (per unit time) by 
the stream of molecules incident upon the sur- 
face, E, is the energy carried by the molecules 
reflected (or re-emitted) from the surface, and 
E,’ is the energy that would be carried by the 
reflected molecules if they were issuing as a 
stream out of a body of gas at equilibrium at 
the temperature of the surface.? This accom- 
modation coefficient is therefore the ratio of the 
net energy actually flowing across the surface to 
the value the latter would have if the incident 
stream of molecules were unchanged but the 
reflected stream were replaced by a stream in 
“thermal equilibrium’’ with the solid. A theo- 
retical relation between g’ and Knudsen’s accom- 
modation coefficient has not as yet been obtained 
and probably no simple relation between them 
exists, because of the complex nature of the 
actual stream of incident molecules. However, 
it is possible to relate g’ to a slightly different 


6 It is assumed that, as in the case in the present measure- 
ments, the length of the mean free path in the gas is small 
relative to the diameter of the cell wire. 

7 The “pressure” of this gas is to have a value such that 
the issuing stream contains the same number of molecules 
as the actual incident and reflected streams. 
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ACCOMMODATION COEFFICIENTS 


constant, defined as 
a’ =(E;—E,)/(Ei —E,’). (6) 


E;, E,, and E,’ have been defined following Eq. 
(5),2 and E,’ is the energy the incident stream of 
molecules would carry if the temperature gradient 
in the body of the gas continued uniformly to 
the surface (in which case the temperature dis- 
continuity at the surface would be 67 of Eq. (1)). 
It is suggested that a’ be called the Maxwell- 
Knudsen accommodation coefficient, since the 
definition is based upon the work of both men.‘ 
The expression relating g’ and a’ is* 


K 2-—a’ 
¢'= @eRMT)(— ~—)/( ), (7) 
2C,+R a’ 


where R is the molar gas constant, M is the 
molecular weight of the gas, and K and C, are 
the thermal conductivity and heat capacity (at 
constant volume) of the gas at the temperature 
T(°K) of the surface. On introducing the value 
of R and appropriate conversion factors, and 
solving for a’, Eq. (7) becomes: 





0.583 C,+0.994 7 | 
mT 
(MT)! K 


a’ =| (1/2)+ 


where the units are: C,, calories per degree per 
mole; K, calories per centimeter per degree per 
second ; and g’, centimeters times centimeters of 
mercury. 

The only approximation of any significance in 
the derivation of Eq. (7) is the assumption that 
the heat capacity for internal energy (excluding 
translation) of a stream of molecules issuing 
from a body of gas at equilibrium is the same 
(per mole) as for the entire gas. This assumption 
is not likely to be seriously in error unless there 


TABLE I. Accommodation coefficients of air on 
bright platinum. 





T,*K d(1/K) /d(1/P) q a’ T,°K d(1/K)/d(1/P) 8’ 








94 378 ; 0.62 224 405 83 
104 145 1.11 248 434 97 
139 254 0.90 270 448 107 

299 0.85 280 491 121 
505 0.62 318 547 149 
383 52 0.75 344 576 167 
313 : 0.88 382 574 181 


* However in the present case the stream for E,’ is defined 
‘o contain the same number of molecules as the stream for 
FE,’ (rather than for E;, as in footnote 7). 
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Fic. 1. Accommodation coefficient of air on bright platinum. 


is a strong correlation between the average 
internal energies of the molecules and their 
translational velocities. Account is taken of the 
fact that the heat capacity for translational 
energy is 2R for the issuing stream as compared 
with 3R/2 for the entire gas. 

It seems desirable to recognize that g’ yields 
the accommodation coefficient a’, rather than a. 
Measurements at lower pressures, such that free- 
molecule conduction exists, are generally con- 
sidered to yield the Knudsen coefficient, a. 
However, it should be pointed out that approxi- 
mations have been necessary in the theoretical 
treatment of free-molecule conduction.® A com- 
parison of values of a and a’ would indicate the 
difference between the values of E; and E,’ in 
Eq. (6). 

The accommodation coefficients for air, ob- 
tained by Taylor and Johnston® are given in 
Table I, and are graphed in Fig. 1. Those for 
the eight gases investigated by Grilly and 
Johnston’ are given in Tables IJ-IX, inclusive, 
and are graphed in Fig. 2. Units employed in 
the several tables are as follows: 7, wire tem- 
perature in degrees Kelvin; P, centimeters of 
mercury; K, (apparent), calories per centimeter 
per second per degree X 10°; g’, centimeters X cen- 
timeters of mercury X 10°. 

The pressures, at which apparent K’s for air 
were obtained by Taylor and Johnston are given 
in Table IV of their paper, and are not quoted 
here. Certain of Taylor and Johnston’s results— 
namely those for wire temperatures of 104° and 
139°K—are based on plots of three points that 
were not quite linear. The remaining points 
below 175° were based on runs at only two 
pressures. This probably accounts for the spread 


°E. H. Kennard, reference 4 (Sections 178-179). 
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Fic. 2. Accommodation coefficient of several gases on bright platinum. Experimental points, O @ This research, 
+0 Other investigators (references indicated). 


TABLE IT. Accommodation coefficients of oxygen on of their data below 175°K. For purposes of 
bright platinum. , : : 

_ extrapolation—i.e., to obtain thermal conduc- 

d(1/K) tivities—the slight non-linearity is less serious 

T,°K PrP: Ps Ki Ks _—Ks_d(t/P) «’ a’ ~—than in the determination of accommodation 




















oo ee SS eS eS oe |S lf coefficients, where slopes themselves are used. - 
isg Leas 421 920 2708 240? 2771 403 §0 94? — Comparisons with results obtained by other ff 
i ay 1 Gos 3a Sa 35 4g) Gt 43 investigators” are included in Fig. 2. The ff _ 
He Liss 416 1109 44uz 44@ 448i i473 943 measurements in references 13 to 17 were made 
= 7 ee ee te St tem fe is 87; §6at low pressures, where free molecule conduction rT 
383 1.556 4.08 8.61 7.676 7.868 7.920 750 258 0.63 exists, and presumably yield values of Knudsen’s : 
accommodation coefficient a, while those in : 
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TABLE III. Accommodation coefficients of nitric oxide on references 10 to 12 were made at higher pressures 














2 
bright platinum. (p>1 cm), and should yield essentially the 2 
———— — =e ; 
d(1/K) 1C, T. Archer, Proc. Roy. Soc. (London) A165, 474 3 
a a : . _ ao 1938). ; 
T,°K Pi Pz P K K2 Ks d(i/P) ¢’ _,. = - , , 3 
; : : : ilk. “- 11H. S. Gregory (Spencer-Gregory), Proc. Roy. Soc. 3 
139 1.851 5.72 10.72 2.763 2.781 2.795 940 121 0.33 (London) A149, 35 (1935). 3 
56 1.646 4.28 10.86 3.139 3.153 3.161 440 63 0.63 12 ¢ , : de Phil. Mac ine 
182 1.387 6.46 10.01 3.689 3.717 3.720 364 61 0.76 H. Spencer-Gregory and E. H. Dock, Phil. Mag. 25, = 
199 1.312 5.27 9.97 4.020 4.064 4.070 450 80 0.70 129 (1938). : 
216 2.530 5.89 10.65 4.381 4.398 4.404 396 78 0.77 3 1.. B. Thomas and F. Olmer, J. Am. Chem. Soc. 65, 
235 1.531 6.75 10.48 4.733 4.780 4.785 410 87 0.78 1036 (1943) refe 
253 1.511 5.07 10.53 5.064 5.114 5.129 444 101 0.76 Ms : 776 use 
273 1.697 5.20 10.40 5.461 5.515 5.530 466 113 0.75 W. B. Mann, Proc. Roy. Soc. (London) A146, 77 thi 
290 1.871 4.73 10.44 5.776 5.827 5.842 446 114 0.79 (1934). ; 
335 1S SH) Wee Gay GUS GI SE GIB HM. Knudsen, Ann d. Physi 34,993 (1911), 
337 L711 4.86 10.53 6.621 6.704 6.730 500 147 0.77 ‘6 H. H. Rowley and K. F. Bonhoeffer, Zeits. f. physik. ret 
353 1.756 4.70 10.44 6.912 7.001 7.037 544 165 0.76 Chemie B21, 84 (1933). Recalculated by Thomas and Dj 
369 1.712 4.36 10.56 7.188 7.292 7.341 592 188 0.73  QOlImer, reference 13. ° IC 
383 1.761 4.78 10.58 7.437 7.557 7.597 596 197 0.72 177. Amdur, M. C. Jones, and H. Pearlman, J. Chem and 
4 ’ - s J wy € . « « > Je , . 





ee ~—~Phys, 12, 159 (1944). 
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accommodation coefficient a’.18 Both types of rated with gas, with the exception of some of 
measurements refer to a platinum surface satu- the measurements of Mann!‘ which were made 
on a gas free surface.!® 


TaBLE IV. Accommodation coefficients of hydrogen on 


bright platinum. TABLE VII. Accommodation coefficients of carbon 


monoxide on bright platinum. 











d(1/K) 
Kz d(i/P) 2’ 


13.87 
15.38 
19.83 
22.41 
24.77 
27.02 
29.20 
30.61 
31.37 
32.94 
35.40 
37.48 
39.52 
39.55 
41.14 
43.95 
45.11 
46.70 
48.31 





a, 


d(1/K) 
Ks; d(i/P) 


790 
550 
350 
350 
560 
410 
510 
440 
440 
430 
440 
445 
450 
450 
520 
460 
470 
550 
650 
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TABLE V. Accommodation coefficients of carbon dioxide on 


bright platinum. , : 
~ Accommodation coefficients of methane on 


bright platinum. 






























































d(1/K) 
Ke K3 d(i/P , 
3 d(1/P) a a(1/K) 
4.37 2.069 875 0.46 Ki d(1/P) 
4.24 2.273 2.280 710 0.53 
4.95 2.516 2.522 620 0.59 4. 2.438 y 
4.17 2.739 2.747 900 0.48 4, 3.174 Y 
5.17 2.990 2.996 590 0.61 4. 3.593 . 
4.66 3.229 3.238 590 0.62 4.2: 3.979 . 
4.05 3.487 3.496 570 0.65 4. 4.394 ‘ 
| 4.08 3.730 3.743 580 0.62 3.98 4.775 0.73 
es 4.08 4.007 4.021 610 0.64 4.96 5.203 0.62 
:duc- 3.30 4.169 4.186 580 0.64 4.41 5.640 0.73 
4.31 4.616 4.637 670 0.57 5.06 6.061 0.68 
rious 4.71 4.891 4.912 700 0.56 5.95 6.492 0.69 
: 4.63 5.229 5.257 710 0.56 5.06 6.948 0.67 
ation 4.21 5.567 5.594 620 0.61 4.34 7.406 0.74 
c 5.17 7.930 0.78 
od 4.95 8.853 0.74 
. : ; : ; ; 4.64 9.374 0.67 
other TaBLE VI. Accommodation coefficients of nitrous oxide on 5.53 9.974 0.63 
bright platinum 5.30 10.535 0.68 
The gntp ‘ 5.07 11.112 0.58 
made d(1/K) 
ction T,“K Pi Pz € Ks d(1/P) a’ TABLE IX. Accommodation coefficients of helium on 
isen’s 2415 4.19 2.171 0.50 bright platinum. 
f 1.772 4.77 2.434 0.53 
se in 1.906 4.57 2.756 0.71 
1.962 4.79 2.864 0.54 d(1/K) 
sures , 4.92 3.167 0.65 Tx -P P2 E , 
th 1847 4.52 3.415 0.61 . a 
e . . 3.661 0.61 1.885 4.46 15.65 220 0.44 
1.448 4.31 : 3. 9: 3.949 0.60 2.054 4.88 16.94 200 0.51 
: 1.748 4.78 i . 193 4.207 0.61 2.219 5.46 19.78 256 0.49 
5, 474 “ 1.611 6.04 <. te e ms 4.498 0.59 2.188 5.42 21.27 382 0.44 
1.591 5.15 y .733 4.814 0.58 2.139 5.32 22.95 3 359 0.51 
. 2.367 5.12 11.20 5: 5.132 0.54 2068 5.11 24.52 427 0.48 
. Soc. 373 2.049 5.30 11.13 5. ‘ 5.484 0.53 2078 5.06 463 0:50 
2.203 5.67 10.96 5.731 5. 5.814 0.57 2139 5.27 475 0.54 
ig. 25, 2.228 5.49 521 0.52 
2.178 5.30 0.52 
- * However, the accommodation coefficient calculated in 2.019 4.98 0.52 
. 65, le : : : ge 2.113 5.07 0.53 
rences 10 to 12 is not identical with a’ because the theory 2073 5.13 0.53 
6, 776 used by these workers differs from that presented here. In 2.188 5.38 0.53 
this connection, see the comment of Kennard (reference 4, yo ry 0.54 
Pages 323-324), on the formula used by Dickins, reference 1655 4.31 ‘ : on 
—" 20 (Proc. Roy. Soc. (London) A143, 517 (1934)). Gregory 380 1.775 4.60 36.04 38.85 40.03 0.48 
's and (reference 11) presumably used the same formula as 
ickins, while Archer (reference 10) and Spencer-Gregory §=————— 
Chem. and Dock (reference 12) used a different theory developed 19 However, the work of Amdur, Jones, and Pearlman 


by one of them. (reference 17) indicates that in some earlier work, including 








440) GRILLY, 


In the figure the ordinate scale for each gas 
has been shifted for better distribution of the 
points, but qualitatively all the gases are in 
order of the magnitudes of their accommodation 
coefficients. 


DISCUSSION OF RESULTS 


Our own values of a’ agree reasonably well 
with the values of a obtained by Amdur, Jones, 
and Pearlman,’ which are probably the most 
reliable values of a (on a gas-saturated platinum 
surface) in the literature. This comparison is 
perhaps the best evidence now available that 
these two accommodation coefficients do not 
differ very greatly. 

As shown by Fig. 2, the literature is generally 
lacking in measurements of accommodation 
coefficients below room temperature—hydrogen 
(and some data with helium, on glass) being 
the only gas for which data are published. The 
accommodation coefficients for hydrogen rise 
with descending temperature—as confirmed in 
our own work with this gas. It is generally 
assumed that the accommodation coefficients of 
all gases—on all surfaces—tend toward unity as 
the temperature is lowered. This is because the 
surface becomes covered with an adsorbed layer 
of the gas which is thought to be conducive to 
more efficient energy interchange. However, our 
own results with gases other than hydrogen 
consistently indicate sharp decreases in accom- 
modation coefficients at temperatures somewhat 
below that of the room. This effect is, in fact, 
most marked for those. gases that are most 
strongly adsorbed. We plan, sometime, to meas- 
ure the accommodation coefficient of hydrogen 





that of Thomas and Olmer (reference 13), the pressure of 
the gas may have been too low to attain complete satura- 
tion of the surface. Accommodation coefficients measured 
in the region of partial saturation would be expected to bea 
function of the pressure of the gas. 


TAYLOR, 


AND JOHNSTON 


in the neighborhood of its boiling point to dis- 
cover whether hydrogen also will behave like 
the other gases in its region of high adsorption. 

Some differences among results by different 
investigators is probably to be expected since 
we are dealing with a surface effect and surface 
conditions are probably not exactly comparable. 

The observed accommodation coefficient is 
dependent on the average efficiency of transfer 
to the surface of the several kinds of molecular 
energy (translational, rotational, and vibrational) 
which may be excited in the gas. Some experi- 
ments of Knudsen” indicate that the accom- 
modation coefficients for translational and rota- 
tional energies in hydrogen are essentially the 
same. 

However, it is probable that in general the 
efficiency of energy transfer is somewhat different 
for internal degrees of freedom (particularly for 
vibration) than it is for translation. This has 
been considered by Shafer, Rating, and Eucken,” 
who assign much lower coefficients to vibrational 
and rotational than to translational energy. For 
example, for carbon dioxide at 3°C they assign 
an accommodation coefficient of unity for trans- 
lation and rotation but of only 0.13 for vibration. 

A relatively low accommodation coefficient 
for vibrational energy may partially explain the 
decrease of the accommodation coefficient for the 
total energy which is generally observed at 
higher temperatures. However, it cannot explain 
the decrease at low temperatures observed for 
several gases in this research. 

We see no obvious explanation of the shapes 
of the curves at low temperatures, and regard it 
as desirable to investigate the effect further with 
other gases, and with hydrogen and helium at 
liquid hydrogen temperatures. 


20 W. Knudsen, Ann. d. Physik 6, 129 (1930). 
21K. Shafer, W. Rating, and A. Eucken, Ann. d. Physik 
42, 176 (1942). 
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Statistical Mechanics of Multimolecular Adsorption 


II. Localized and Mobile Adsorption and Absorption 
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It is shown that the transition from localized to mobile adsorption (in the first adsorbed layer) 
takes place at rather low temperatures for potential barriers of around 1000 cal./mole or less, so 
that localized physical adsorption should be a very rare phenomenon at the temperatures 
usually employed in adsorption experiments. Since the BET (and Langmuir) isotherm actually 
assumes localized adsorption, a new approximate isotherm equation is derived here on the 
basis of a mobile first layer obeying a two-dimensional van der Waals’ equation. This isotherm 
is in semiquantitative agreement with the two-dimensional phase changes observed by Jura and 
co-workers, and by earlier workers, at very low pressures. The localized-mobile transition in the 
absorption of hydrogen by metals is mentioned briefly. The relationships between (1) the two- 
dimensional equation of state of a monolayer, (2) the adsorption isotherm of the monolayer, and 
(3) the ordinary three-dimensional equation of state of the gas being adsorbed are discussed in 


some detail. 


I. INTRODUCTION 


* most cases at ordinary temperatures, a 
molecule which is held to the surface of a solid 
by van der Waals forces is able to migrate rather 
freely over the surface.'* In order to do this the 
molecule must, in general, be able to pass over 
potential barriers (of the order of 400 cal. /mole™ 
in physical adsorption) which arise because of the 
periodic structure of the surface of the solid and 
the resulting periodic variation (in two dimen- 
sions) of the energy of adsorption. However, at 
sufficiently low temperatures, such a molecule 
will be restricted to a localized!” site at a po- 
tential minimum, with vibrations about the 
equilibrium point replacing translational motion 
in two dimensions. At intermediate temperatures, 
the situation is rather more complicated. 

We have been referring above to molecules in 
the first adsorbed layer. Whether these molecules 


_ “For a recent discussion of this point, see S. Brunauer, 
The Adsorption of Gases and Vapors (Princeton University 
Press, Princeton, 1943), pp. 448-454. 

The term “localized” will be used here to refer to 
molecules whose motions in the two surface dimensions 
can be described as primarily vibrational, with only rela- 
tively infrequent shifts from one potential minimum to 
another. These shifts can of course occur in any case by 
the process of evaporation and condensation. By ‘“‘mobile”’ 
molecules will be meant molecules whose motions in the 
{wo surface dimensions are essentially those of a two- 
mensional gas or liquid, with no effective constraint 
ue to the periodic variation of the energy of adsorption 
over the surface. The energy of adsorption itself of course 
constrains the motion to two dimensions. 


are localized or mobile, the behavior of molecules 
(if any) in second and higher adsorbed layers is 
generally assumed to resemble the behavior of 
molecules in the liquid state, and hence these 
latter molecules are not localized. However, at 
temperatures below the melting point of the 
adsorbed phase, this would naturally not be the 
case. 

As is well known, mobile monolayers exhibit a 
phase change which is analogous to ordinary 
condensation. This phase change is usually 
studied experimentally on the surface of a liquid, 
but there have been a few studies of two- 
dimensional condensation on the surface of a 
solid.? Presumably these experiments involving 
solid surfaces have always dealt with mobile 
rather than localized monolayers. The theoretical 
problem involved in studying the condensation 
of a mobile surface layer is essentially the same as 
for the usual type of condensation. A completely 
rigorous treatment has not been accomplished 
because of mathematical complications, but cer- 
tain semi-empirical methods are quite useful. 

It is clear from theoretical considerations’ that 
“condensation” of a localized monolayer is also 
possible. The theory of the phase change in this 
case is a good deal simpler (though still requiring 

2 See reference la, pp. 437-440. Also, see reference 10. 

3 R. H. Fowler and E. A. Guggenheim, Statistical Thermo- 


dynamics (University Press, Cambridge, England, 1939), 
pp. 429-443, 
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approximations) than for a mobile monolayer, 
essentially because of the restriction to discrete 
sites which makes it somewhat easier to enumer- 
ate possible states and energies of the system. A 
statistical treatment leads to Langmuir’s isotherm 
if interactions between adsorbed molecules are 
neglected, and to a sudden phase change, at 
temperatures below a critical temperature, if 
interactions are not neglected. As far as the 
writer is aware, the condensation of a localized 
monolayer has never been studied experimentally. 

The general problem is thus rather complicated. 
Not only are different phases possible for an 
adsorbed layer (for example, gas, liquid ex- 
panded, liquid condensed, solid), but one must 
also consider (1) the low temperature localized- 
mobile transition, (2) the possibility of multi- 
molecular adsorption at higher pressures, and, in 
this case, (3) the properties of second and higher 
layers. 

In Section II we shall consider the transition 
from localized to mobile adsorption when the 
adsorbed phase is quite dilute. A few remarks will 
also be made about the absorption of hydrogen 
by metals, since this is a closely related problem 
in some respects. In Section III we derive an 
approximate adsorption isotherm for multimolec- 
ular adsorption when the first layer is mobile. 
Relations between the two-dimensional equation 
of state and the adsorption isotherm and between 
the two- and three-dimensional equations of 
state are discussed in Sections IV and V, re- 
spectively. 


II. THE LOCALIZED-MOBILE TRANSITION 
Adsorption 


We consider in this paper only the simplest 
case. Several possible generalizations will be 
obvious. We assume that: (1) The potential 
energy function restricting free migration of an 
adsorbed molecule over the surface is of the 
form 


2ry 
U= V(1-cos =) +4 V(1=cos : 


2ax 2ry 
3V(2-c0s cos , (1) 
d d 


49 If the dimensions of the molecule are large compared 
to d, the situation is more complicated. But there will still 
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where the maximum potential barrier is Vo=2V, 
and d is the distance between potential minima 
along both x and y axes (the surface lattice is 
thus assumed to. be simple cubic); and (2) the 
number of molecules, N, adsorbed in the area 
%=/? is sufficiently small so that (a) interactions 
between adsorbed molecules may be neglected 
(the mobile case may then be treated as a perfect 
gas and the two-dimensional gas pressure will be 
less than the two-dimensional vapor pressure of 
the condensed phase—solid or liquid—so con- 
densation will not take place when mobility is 
acquired, as will ordinarily be the case), and 
(b) the number of sites available to each of the NV 
adsorbed molecules in the localized case is es- 
sentially n? =/? /d?. That is, the probability of any 
two molecules being in the same site is assumed 
negligible, and it is not necessary to invoke a 
special “exclusion principle” restricting the num- 
ber of molecules per site to one. Although the 
following discussion will be based on the above 
simplifications, the general conclusions arrived at 
should not be affected by further refinements. 

It is easy to see that under assumption (2), 
above, a partition function for the system of the 
form 


1 
a (2) 


nr’ = 


where f is the partition function for a single 
molecule over the entire area YW, is correct not 
only for the mobile case but also for the localized 
case. In the latter case we would ordinarily write 


n*! 


=—___—_—_—-f,, (3) 
N!'(n?—N)! 


where f, is the partition function for a single 
molecule in a single site. Since we are assuming 
n*>> N, Eq. (3) becomes 


N 1 
gu ah v= (ntf (4 


which is the same as Eq. (2) with f=n°f,, thus 
taking into account correctly the degeneracy (n’) 
in the localized case. 

be a periodic potential function. Other things being equal, 


the potential barrier should be smaller the larger the 
molecule. 
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KX 


We proceed now to formulate f for the two 
degrees of freedom of interest here. In order to do 
this, a study must first be made of the possible 
energy levels of the system, the system being a 
single molecule of mass m in the area Y=? with 
the potential energy given by Eq. (1). The two- 
dimensional quantum-mechanical problem here 
is essentially equivalent to the one-dimensional 
problem encountered in studying restricted rota- 
tion in molecules such as ethane. This makes it 
possible to use some of the results and tables 
obtained by Pitzer and Gwinn," after a suitable 
correlation of parameters. Also, it appears that a 
rather accurate method of approximation sug- 
gested by Pitzer and Gwinn should be quite 
adequate and very convenient for adsorption 
studies, so we shall consider it in some detail as 
applied to the present problem rather than dis- 
cuss further the more exact results obtained by 
these authors. The approximation consists of 
using the following partition function: 


Jucomente osc— quantal 
S =ferassical x . (8) 


harmonic osc—classical 





The classical partition function here is 


i. 
mh USS 


Xexp (—H/kT)dxdydpdp,, (9) 


where 


pb py? 2ax 2ry 
W= + +4V( 2-cos ——cos — . (10) 
2m 2m d d 


On performing the integrations over the mo- 
menta, we have 


2rmkT 
Sotess = : exp (— V/kT) 


12 


' 1 V 2ax . 
x| f exp (= —c0s ax, (11) 
0 2kT d 


2amkT ; (—V/eTVI ( V (12) 
oa 2 exp (— 2] —— }, 12 
h " — 


9 
« 
lee 


“’K. S. Pitzer, J. Chem. Phys. 5, 469 (1937); K. S. 
Pitzer and W. D. Gwinn, J. Chem. Phys. 10, 428 (1942). 
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Fic. 1. Heat capacity of adsorbed argon. 


where Jo(u) is a modified Bessel function of the 
first kind. Also, 


Sues ose—quant — n*/[1 — exp ( —h V, ‘Rk T) }*, (13) 
Suse ose—class = n*(kT/hyv)?, (14) 


with 
v= (V/2md?)!. (15) 
Substituting Eqs. (12)—(14) into Eq. (8) gives 
2an*ue—*"Iy?(u) 
— sm (16) 


K=(2h?/md?V). (17) 





u=V/2kT, 
The two limiting expressions are, as expected: 
f=an?/(1—e**), (18) 


2amkT 
usmall: f= i?, (19) 


9 
~ 


u large: 


Equation (18) follows from 


lim Io(u) = . (20) 
wo" (Qe)! 


and we have used *d?=/ to obtain Eq. (19). 
Finally, the complete partition function for the 
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two degrees of freedom is 


1 
ae (21) 


The contribution of these degrees of freedom to 
the thermodynamic functions may be obtained in 


the usual way. For example, for the heat capacity, 
one finds 


T,(u) T;*(u) 
c= Nel —1 +2uf 1 —_——-— | 
Io(u) = I?(u) 
2K*u*e *) 


+—————;. (22) 
(1—e-K»)2| 


As an illustration we have plotted C against T, in 
Fig. 1, for different values of Vo=2V, for argon 
adsorbed on potassium chloride. The true value 
of Vo is probably about 400 cal./mole.> Actually, 
to save labor, Table IV of Pitzer’s*® paper was 
used instead of Eq. (22). The connection between 
our notation and his is found by putting 


I (moment of inertia) = mr? = ml? /4r°, 


since /=2zar. This becomes clear if one considers 
the problem of a particle constrained to move in a 
circular path. The motion can be thought of 
either as rotation or translation (even if a 
perfectly reflecting barrier is introduced, blocking 
the path). One finds that K is related to the 
variable n?/ITV used by Pitzer by 


10-%6n2 /T V = 31.25K°, (23) 


where V is in cal./mole. © 

The theoretical form of the heat capacity 
curve, as given in Fig. 1, has apparently not been 
observed experimentally. In fact the only work of 
this type seems to be that of Simon and Swain'® 
in which the heat capacity of adsorbed argon on 
charcoal was measured. This experimental heat 
capacity curve resembles that of a one-dimen- 
sional harmonic oscillator. There have been 
several conflicting interpretations of this result. 
In any case, this example is unsatisfactory from 
our point of view because of the high porosity of 
the adsorbent and the large amount of argon 


5W. J. C. Orr, Trans. Faraday Soc: 35, 1247 (1939). 
The curves given are not exact, of course, since the theo- 
retical potential function used is more symmetrical than 
the actual (potassium chloride) one. 





adsorbed, leading to complications not taken into 
account in the present theory. 

The transition from vibration to translation is 
continuous and takes place over an appreciable 
temperature range. However, one may indicate 
the general temperature region in which the 
transition occurs by noting that the heat capacity 
for the two degrees of freedom passes through a 
maximum at about T= Vo/10R, this temperature 
being fairly insensitive to K (values of K of 
interest here will correspond in general to 
10-**n?/TV <8). It is therefore clear, since V4 is 
generally of the order of 500 cal./mole or less, 
that the transition is in progress at quite low 
temperatures and that true localized physical 
adsorption should be a very rare phenomenon 
indeed at temperatures ordinarily used in adsorp- 
tion work, although intermediate cases should be 
somewhat more common (e.g., around 90°K). 
Certainly the mobile case is the rule and not the 
exception. Thus the above analysis tends to con- 
firm the similar view expressed by Brunauer' 
concerning the mobility of physically adsorbed 
molecules. In chemical adsorption the situation 
is, of course, radically different because of much 
higher potential barriers (see the discussion of 
absorption below). 

At the low pressures being considered here the 
adsorption isotherm is linear, the amount 
adsorbed being proportional to the pressure. The 
constant of proportionality is in turn pro- 
portional to the partition function f, Eq. (16). It 
is through the partition function f that the 
amount of adsorption at low pressures and tem- 
peratures will depend on the localized-mobile 
transition. 


Absorption 


As is made clear by Fowler and Guggenheim," 
the absorption of hydrogen by certain metals 
bears some resemblance to the adsorption of 
molecules on the surface of a solid. Both localized 
and mobile absorption are to be expected and 
both are found experimentally. We should merely 
like to point out here that there is an analogous 
localized-mobile transition problem in absorp- 
tion, but in three dimensions instead of two. 
Since localized absorption is very common, the 


6 See reference 3, pp. 554-563. 
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potential barriers must be much greater in this 
case. This is caused, at least in some cases, by 
chemical absorption—the formation of metal 
hydrides, MH,, at the absorption sites (metal 
atoms). 


Ill. MULTIMOLECULAR ADSORPTION ON A 
MOBILE FIRST LAYER 


The Brunauer-Emmett-Teller’? theory of multi- 
molecular adsorption is a_ generalization of 
Langmuir’s theory of unimolecular adsorption. 
As is pointed out in the statistical versions,°° 
these theories are based on the assumption of 
localized sites. Since localized physical adsorption 
is really rather uncommon, it would seem worth- 
while to derive a multimolecular adsorption 
isotherm on the assumption of a mobile first 
layer. We shall confine ourselves in this paper to a 
treatment which is only slightly more general 
than the corresponding statistical treatment? for 
a localized first layer, which leads to the BET 


“equation. The BET model neglects interactions 


between molecules in the same layer. We shall 
also neglect these interactions, except in the first 


"See reference 1a, Chapter VI. 
oR. H. Fowler, Proc. Camb. Phil. Soc. 31, 260 (1935). 
T. L. Hill, J. Chem. Phys. 14, 263 (1946); A. B. D. 
Cassie, Trans. Faraday Soc. 41, 450 (1945). 
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layer where they will be taken into account in an 
approximate manner (by the use of Eq. (24)). 

We assume for simplicity that the X molecules 
in the first layer obey van der Waals’ equation in 
two dimensions. It should be recognized at the 
outset that the choice (there seems to be no 
satisfactory alternative at the present time) of 
such a crude equation of state will give the 
resulting theory only semi-quantitative signifi- 
cance. However, the theory should still be useful 
in contributing to an understanding of the 
problem. We also neglect the existence of more 
than one liquid phase in using this equation of 
state. Thus we shall write 


(g+a’X?/%*)(A— Xb’) =XkT, (24) 


where ¢ is the two-dimensional pressure, % is the 
area, and a’ and b’ are constants. The associated 
partition function (for these two degrees of 
freedom) which leads to this equation of state is 


1 ¢24amkT x 

|= ——a-x0 exp (a’/X wn) | . (25) 
We therefore write for the separate partition 
functions of the X molecules in the first adsorbed 
layer and the N—X molecules in higher layers, 
respectively, 


1 
Qs= =| ii exp (€1/kT) 


a 


2amkT 4 
x — exp (a’/X ae) » (26) 
12 
(N—1)! 


Co XK 1)! 





Lirexp (e,/kT) ’-*, (27) 


in which 7 is the internal partition function, j, the 
partition function for vibrations normal to the 
surface, — €, the potential energy of a molecule in 
the first layer due to interaction between the 
molecule and the surface, and jz, and —ey, the 
partition function and potential energy of mole- 
cules in higher layers (assumed to be liquid-like). 
Equation (27) is the same as Eq. (4) in the first 
paper of reference 9, while Eq. (26) is the analog 
of Eq. (3) in that paper. We are thus assuming 
that molecules in the first layer behave like a two- 
dimensional van der Waals’ gas and that addi- 
















0.1 x 0.5 
Fic. 3. 
tional molecules may build up higher layers on 


top of the X sites presented by the first layer. 
The equilibrium value of X is found from 


0 log QsQx 0 
ax , 


(28) 


We neglect unity compared to N and X in 
Eq. (27) and find, from Eq. (28), 


(N—X)(%—X0’) Xb’ 2a X 
oo (-gtet) 

XxX? A—Xb’ WAkRT 

ju exp [(er—«1)/kT ] 

ij, QemkT/h? 








(29) 





The chemical potential of the adsorbed phase is 


A 8 log OsQz 
kT aN 

N-X 
(30) 


| ts 
= log —log 71 -——. 
RT 





For the gas phase, 


Me po(T) 
pe +log DP, 
kT kT 





(31) 









where p is the gas pressure and yo(7°) need not be 
specified further for our purposes. At equilibrium, 
HA=hG: 

N-X €L 


—log j -—— = 
"BT 


wo(T) 
kT 








log +log p. (32) 
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For the pure liquid 
_  €4 po(T) 

i oe aa 

kT 

where /o is the vapor pressure of the liquid at 7, 


Subtracting Eq. (33) from Eq. (32) gives 
(N-—X)/N=p/po=x 








+log po, (33) 













(34) 





X=N(1—x). (35) 


If we now substitute Eq. (35) into Eq. (29) to 
eliminate X, there results 


6(1—x)? 
~ x(1 — 6+ 6x) 










6(1—x) 
xXexp F 
1—6+ 6x 











—a(1 -»)| (36) 












b' jin 2rmk T 
C= 





exp [(e1:—€z)/RT ], (37) 
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UL i” 


6= Nb’ /A, 





a=2a’/b’RT. 





The quantity @ corresponds to v/v» in the BET 
theory. It is not possible to solve Eq. (36) for 6 
or x explicitly. In order to calculate an isotherm, 
for given values of C and a, one must find by trial 
the value of x which satisfies Eq. (36) for each 
value of @ chosen. 

There are several special cases of interest: 

(1) At low pressures (x1), Eq. (36) simplifies 
to 











: & 0 
x =— —— exp (—-00). (38) 
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(34) 
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It should be emphasized that this equation can 
also be derived directly from Eq. (26) alone, if 
only unimolecular adsorption is allowed. That is, 
it follows from Eq. (24) and does not depend on 
any of the assumptions concerning the behavior 
of layers other than the first. In most cases of 
interest, Eq. (38) is an excellent approximation 
to Eq. (36) up to values of @ of about 0.7 or 0.8. 
Then second and higher layers begin to build up 
and @ becomes increasingly larger than Eq. (38) 
would predict. We have plotted, in Fig. 2, 6 
against Cx for several values of a, using Eq. 
(38). The critical value of a is a,-=27/4, or 
T.= 8a’ /27b’k. This is, of course, also the critical 
temperature of the equation of state, Eq. (24). The 
other critical values are 6.=4% and x,=0.0869/C. 
These results are of interest in connection with 
the recent communication of Jura, Loeser, 
Basford, and Harkins'® reporting a first-order 
phase change for n-heptane on the surface of 
silver at 15°C and p=0.019 mm (x=0.00070). 
Clearly, the present theoretical treatment is in 
semiquantitative agreement with their results 
(see also reference 2). As we have arbitrarily 
chosen van der Waals’ equation as the equation 
of state, one should not expect more than this. 
Since C will generally be of the order of 100 or 
more, x, is in the neighborhood of 10-*-10~+. 
Hence this phase change is to be expected at 
quite low pressures ordinarily, as is emphasized 
by Jura and co-workers on the basis of their 
experimental work. 

(2) If we set a’ =0 and b’ =0, the first layer is a 
perfect gas. Equation (36) then reduces to 


WC x 


N=— ——_, (39) 
b’ (i-—x)? 
a result obtained by Cassie" by an incorrect 
argument.® 
(3) If the exponential factor in Eq. (36) is 
ignored, one obtains, on solving for @, 


6=Cx/(1—x)(1—x+Crx), (40) 


which is the BET equation. 
Returning now to the general isotherm, Eq. 
(36), we have plotted in Fig. 3 as an illustration 


© G. Jura, et al., J. Chem. Phys. 13, 535 (1945). See also 
reference 20. 
" A. B. D. Cassie, Trans. Faraday Soc. 41, 450 (1945). 
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the isotherm for C=490 and a=8. The dotted 
curves are ordinary BET curves for comparison, 
with c=490, 81, and 27. The solid points above 
the a=8 curve are for a=12 and those below are 
for a=4. It will be observed that Eq. (36) is not 
very sensitive to a (except at low pressures). It 
is also not very sensitive to C, when x>0.1. For 
example, in this region the points for a=4 and 
C=490 are also approximately the points for 
a=8 and C equal to about 10. However, the 
sensitivity to C at low pressures is considerable, 
as is indicated in Fig. 2. 

Perhaps the most noticeable difference be- 
tween the mobile isotherm, Eq. (36), and the 
BET isotherm is that the former associates much 
lower pressures with values of @ less than, say, 
0.8. In general this is in the direction of a 
correction of a well-known fault of the BET 
equation, but unfortunately Eq. (36) seems to 
over correct considerably. It should be remarked 
that this property (low pressures) of Eq. (36) is 
caused primarily by the fact that we are taking 
into account interactions between molecules in 
the first layer and not to the fact that the 
molecules are mobile. As will be seen in a later 
paper of this series, if one takes these interactions 
into account in a localized first layer a con- 
densation can take place® at pressures which are 
also very low. But in this localized condensation 6 
jumps from almost zero to almost unity, and 
6.=}. 

While the equations of this section should not 
be expected to provide a really satisfactory 
agreement with experiment, it is of some interest 
to know to what extent they are actually appli- 
cable and useful. In a paper to be published at a 
later date, Dr. George Jura will discuss this 
question, using experimental data obtained at the 
University of Chicago. 

Finally, the question arises as to whether a 
more satisfactory equation of state might be 
used. A possibility that suggests itself is the 
equation of Devonshire.” Actually, the general 
form (Fig. 4) of the adsorption isotherm thus 
obtained at low pressures (corresponding to 
Eq. (38)) is not nearly as satisfactory as Eq. (38) 
itself, obtained from van der Waals’ equation, 
judging from the experimental curves of Jura and 


22 A. F. Devonshire, Proc. Roy. Soc. 163A, 132 (1937). 





448 


co-workers (which are the only ones available up 
to the present time). For example, 6, (Devonshire) 
=0.71, which appears to be much too high. The 
adsorption isotherm is derived in Appendix I. 


IV. RELATIONS BETWEEN THE TWO-DIMEN- 
SIONAL EQUATION OF STATE AND THE 
ADSORPTION ISOTHERM 


The equation of state of an adsorbed phase can 
be related thermodynamically to the adsorption 
isotherm through the Gibbs equations (Appendix 
II). However, thermodynamics can tell us only 
very little about the equation of state and 
adsorption isotherm separately. The following 
equations represent the best that can be done 
(the notation is summarized in Appendix IT): 


0A 
p(t, T) =exp | (—) —ary| /er], (41) 
ON/ 4,7 


g(', T) = —(0A/d%)y, r. (42) 
Thermodynamics does not provide A (2, N, T) or 
wT). In order to proceed we must turn to 
statistical mechanics. 

We shall confine the discussion 
monolayers as this is the case in which ¢ is of 
most interest. A corresponding treatment can be 
given for other systems, using the appropriate 
partition function. We have, 


A=-—hkT log Q, (43) 


where Q is the partition function. The problem 
thus becomes one of finding-Q. In analogy with 
Eq. (26), we employ here 


1 
o- J ii, exp (€:/kT) 





2amkT 1 
Je 


9 


2 
l N 
=—O.0,, (44 
7" Q;, (44) 
where Q, is the usual configuration integral 
on= f- :  femnrdsady,: --dxndyn, (45) 


and W(x1,-:--yy) is the potential energy of 


interaction between the N adsorbed molecules. 





to mobile 
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In Eq. (44), Q1 is independent of % so that 





0A d log Q 
gs —-——=eI- , (46) 
a ay 
a) log QO; 
= kT--———_. (47) 
art 


We now express ? in terms of Q,. In Eq. (41), 





0A d log Q 
— = —kT——_, (48) 
aN aN 
d log Q, 
=#7 (tog N—log -——~) (49) 
Ol 


(assuming that Q; is independent of N) and, as 
is well known,’ 


h 3 


——= —log 


p(T) (2amkT)ikT jg 
| (50) 
kT 


where jg is the internal partition function of a 
molecule in the gas phase. Substituting Eqs. (49) 
and (50) in Eq. (41), there results 

j Njn exp (e:/RT) | 

jg (2xamkT)*kT 


d log Q,; 
=Nexp (-— ), (51) 
ON 





p'=p 





which defines p’. When the internal degrees of 
freedom of adsorbed molecules are not appreci- 
ably influenced by interactions with the surface, 
j=jca. It is often convenient to rewrite Eq. (51) 
in the form (compare Eqs. (36) and (37)): 


p’ =—{ — ——— exp [(a—ez)/kT] 


Polit h? 
d log Q, _ 
= N exp (- — =), (52) 
aN 


in which we have made use of Eq. (33). 
Equations (47) and (51) give ¢ and # in terms 
of Q,. Q, depends on N and Y and on the inter- 
molecular forces between the N adsorbed mole- 
cules. If the nature of these forces is known, it is 
possible in principle to derive Q, from Eq. (45). 
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However, because of mathematical difficulties, 
approximations must be resorted to in practice. 

When the surface does not effectively influence 
the intermolecular forces between adsorbed mole- 
cules, the right-hand members of Eqs. (47) and 
(51) are independent of the properties of the 
surface and depend only on I, 7, and the prop- 
erties of the gas being adsorbed. If the adsorption 
of a given gas is studied on several different 
surfaces, all of which satisfy the above condition, 
it is clear that ¢(f, T) and p’(T, T) should be the 
same for all the surfaces. That is, for surfaces of 
this type, the two-dimensional equation of state 
of a given gas is independent of the surface, and 
the adsorption isotherm I'(p, JT) depends on the 
surface only to the extent that in plotting IT 
against p, p must be multiplied by an appropriate 
constant (which depends on the surface) for each 
surface in order to make the different adsorption 
isotherms coincide. This might be the basis of an 
approximate method of determining the surface 
area of a solid from the known surface area of 
another solid“ (or of determining the relative 
surface areas of several solids). For, if both 
surfaces are in the class being discussed, any 
unusual features of the adsorption isotherms 
(e.g., phase changes) will occur at the same value 
of f= N/%. Hence if Ni /%l1= Ne/%s and if Ni, Ne 
and %{, are known, %, may be calculated. 

The question naturally arises as to what com- 
binations of adsorbents and adsorbates are likely 
to exhibit the phenomenon, discussed above, of 
the independence of ¢(I', T) and I'(p’, T) on the 
adsorbent. There will obviously be a wide range 
of examples, depending primarily on the polar 
properties of adsorbent and adsorbate.'5 The 
most unfavorable case, from the present point of 
view, would be the adsorption of a substance 
with a strong permanent dipole moment on a 
plane of an ionic crystal having all ions of the 
same type (e.g., the 100 plane in caesium iodide). 
The permanent dipoles of the adsorbed molecules 
will all be oriented approximately in the same 
direction thus introducing a repulsive orientation 
* See pp. 423-426 of reference 3. For a van der Waals gas 

Q,=(&%— Nb’) exp (a’N2/AkT). 
For a perfect gas, 0, =. 


sae also D. H. Bangham, Proc. Roy. Soc. A147, 175 


" See Chapters VII and VIII of reference 1a. 


effect between them. There will be an additional 
and similar effect because of a further dipole in- 
duced by the ions of the crystal. The attractive 
dispersion forces and the short range repulsive 
forces between adsorbed molecules should be 
relatively unaffected by the surface. But the 
over-all effect of the adsorbent on the forces be- 
tween the adsorbed molecules would be con- 
siderable in this case. On the other hand, the 
most favorable combination would be the adsorp- 
tion of a gas without permanent moments 
(dipole, quadrupole, etc.) on a non-polar surface 
(e.g., argon on graphite, charcoal, metals, mer- 
cury). In such a case the dispersion and repulsive 
forces between the adsorbed molecules should be 
essentially the same as in the gas phase. Without 
an experimental or a detailed theoretical study of 
the question, it is rather difficult to say very 
much about intermediate cases. However, it is 
clear that the most important criterion for the 
independence of ¢g(f,7) and I(p’,T) on the 
adsorbent is that the adsorbent be non-polar 
(non-ionic, no oriented dipoles). This is perhaps 
more crucial than the presence or absence of 
permanent moments in the adsorbed molecules. 

The above argument implies that the most 
suitable gases would be those that are small and 
spherically symmetrical. There is another reason 
for this. Regardless of whether or not the surface 
is polar, a large or non-spherical molecule will 
tend to be oriented on the surface in such a way 
as to give a large ‘“‘contact”’ between the surface 
and the adsorbed molecule (e.g., heptane mole- 
cules would lie flat on the surface rather than 
assume all possible orientations with the same 
probabilities as in the gas phase). This will 
clearly lead to a different intermolecular potential 
function for adsorbed molecules as compared to 
molecules in the gas phase. Also, the extent to 
which a molecule is flattened by a surface will be 
a function of the nature of the surface (even if 
the surfaces considered are non-polar). Hence, 
not only will the intermolecular potential func- 
tion be different (in this case), for adsorbed 
molecules and molecules in the gas phase but the 
potential function will also be different for 
adsorbed molecules on different surfaces, because 
of varying degrees of flattening. 

A further complication is that the most de- 
sirable molecules, such as argon, may not be 
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TABLE I. Two-dimensional critical constants for n-heptane. 








Experimental 


Theoretical 
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Lennard-Jones 























Silver Ferric oxide Graphite and Devonshire Van der Waals j 
Tx 14.9°C <T2. <25°C 29°C 31°C 13°C —3°C | 
we 31042 << 220042 900A? 400A? 82A? 110.42 | 
7 3104? <> ; . ’ y At | 
Gc ¢e>0.18 erg/cm? 0.45 erg/cm? 1.05 erg/cm? 1.9 erg/cm? 0.87 erg/cm? 7 
completely mobile at temperatures sufficiently Experimental evidence on the ratio of the two- 
low to study two-dimensional condensation. dimensional critical temperature to the three- i 
In view of the above restrictions there is dimensional critical temperature, 72-/73-, is not \ 
naturally some question as to just how many, if at all plentiful. There is a certain amount of | 
any, actual systems will satisfy the various re- indirect evidence'® which is contradictory. The | 
quirements (to some reasonable degree of accu- best (until recently) direct evidence”®:'*.' indi- \ 
racy). This question must be settled by experi- cates that two-dimensional condensation does é 
ment. At the present time, experimental data do not occur at temperatures well below 73., but the ( 
not exist for a suitable test of this point. Practi- ratio T,-/T3. was not actually measured. 
cally all of the work in this field (phase changes On the theoretical side, using the method (for S 
on solids) has been done at the University of spherical molecules) which Lennard-Jones and i 
Chicago by Harkins, Jura, and co-workers, but’ Devonshire!® applied successfully to three-dimen- n 
the systems studied so far are not suitable, for the sional condensation, Devonshire” found that, g 
reasons discussed above, to test the arguments of according to the theory, T2./73-=0.53. As pointed t 
this or the following section. out by Devonshire, this value agrees qualitatively t 
with the experimental work mentioned above!”."® t 
V. RELATIONS BETWEEN THE TWO- AND THREE- (as far as that work goes). Of considerably more J 
DIMENSIONAL EQUATIONS OF STATE interest is the fact that, in the condensation of a 
In this section also we shall only consider %-heptane on silver recently reported by Jura and . 
mobile monomolecular adsorption. There has ©0-workers,"® the critical temperature was found 
been some interest in the past in the relation be- to lie between 14.9°C and 25°C. As will be clear . 
tween the two- and three-dimensional critical from the discussion at the end of the preceding . 
temperatures but there has not been as much section, heptane is not a very satisfactory mole- " 
discussion of the more general problem of the cule for this type of comparison; but it is still 3 
relations between the two- and three-dimensional ‘mteresting that, using the known value of 7s. * 
equations of state. The two equations are obvi- =266.8°C, the temperatures 14.9°C and 25°C 7 
ously closely connected. In the ideal case (dis- ead to a value of T2./T3. between 0.53 and 0.55, le 
cussed at the end of Section IV), the law of force in excellent agreement (somewhat fortuitous, no Vi 
between molecules is essentially the same, but the doubt) with Devonshires prediction.”° d 
two equations of state differ in that: (1) in the See po, 440-444 of reference 1a. 
adsorbed phase, the configuration integral Q, “a gee and K. Neugeberger, J. Phys. Chem. 40, 
(Eq. (45)) is ae oe compared to “fF, Simon and R. L. Swain, Zeits. f. physik. Chemie ~ 
#1,**-8v; and (2) in the adsorbed phase, TLE. Seen and A. F. Devonshire, Proc. Roy. : 
Soc.'163A, 53 (1937). ag 
d log Q, 20 After the present manuscript was submitted, a paper cu 
g = kT———_ (53) appeared by Jura, Loeser, Basford, and Harkins (J. Chem. th 
Dy | Phys. 14, 117 (1946)) in which this agreement is also dis- I. 
as compared to on hh ghee pond ghia 
d log Q, ferric oxide and on graphite. These values are given in As 
Paki-———. (54) Table I. They are somewhat more exact than the values ob 
OV published in the above reference. The T2-/T3e ratio for mi 
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Although the theory of Lennard-Jones and 
Devonshire’® is quite successful in predicting 
critical temperatures of gases, it is not very 
successful in predicting critical volumes and 
pressures. Hence, it is perhaps not surprising to 
find that the theoretical critical two-dimensional 
pressure g, and area Y., calculated (see below) 
from the theory and the experimental three 
dimensional values for n-heptane, do not seem to 
agree very well with the two-dimensional experi- 
mental values. This may be seen in Table I?° and 
in the curves given by Jura and co-workers.'®.?° 
We also include in Table 1, for purposes of com- 
parison, the corresponding theoretical values 
(also based on the three-dimensional experimental 
values) obtained if one uses the appropriate two- 
and three-dimensional van der Waals’ equations 
(see below). In the critical region one should not, 
of course, take van der Waals’ equation very 
seriously, but these results are still of some 
interest. Although neither Devonshire’s theory 
nor van der Waals’ equation should lead to very 
good values of the critical area, it is probably 
true that the large discrepancy between the 
theories and experiment is owing at least in part 
to some special cause. That is, as mentioned by 
Jura and co-workers,”® the experimental critical 
areas are really much higher than one would 
expect. 

The van der Waals values in Table I were 
obtained by carrying through the two-dimen- 
sional statistical treatment of the hard sphere 
model with an attractive energy proportional to 
r-* for r>ro. The approximations made were the 
same as in the corresponding three-dimensional 
treatment given by Mayer and Mayer.” This 
leads to the evaluation of the two-dimensional 
van der Waals constants. If the two- and three- 
dimensional equations of state are 


a’ N? 
(+) tao = wer, (55) 


both ferric oxide and graphite is 0.56. The critical areas 
and pressures for heptane on the three surfaces are not in 
agreement. This is certainly not surprising (see the dis- 
cussion at the end of Section IV). One gets the impression 
that 7. is not nearly so sensitive to varying conditions as 
Y%- and ¢y,. 

1 J. E, Mayer and M. G. Mayer, Statistical Mechanics 
(John Wiley and Sons, Inc., New York, 1940), pp. 262-269. 
‘As Mayer and Mayer point out, the correct results are 
obtained here by a convenient but slightly incorrect argu- 
ment (cancellation of errors). 
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aN? 
(P+—)v- Nb) = NkT, (56) 





one finds 
Yo 2 T 3b b 
y= x(-) |=200- ( -) ’ (57) 
2 2\2n 
TuUory? Taf 3b: 
nme. - 
4 4b\ 2x 


where ap is the projected area of a hard sphere of 
radius ro/2 and —wp is taken as the potential 
energy of interaction between two molecules 
practically in contact (r—ro). Using the expres- 
sions for the critical constants in terms of the van 
der Waals constants, one obtains the following 
relations between the two- and three-dimensional 
critical constants: 


T 2. =0.50T 32, (59) 


2 


A 3rfV. 1\3 Ve\3 
—= —.—) -1.38(—) (60) 
N 2\N 2e N 





It will be noted in particular that 72,/T3.=0.50, 
according to this model. 

The corresponding equations for the theory 
of Lennard-Jones and Devonshire are ap- 
proximately : 


T2-=0.53T x¢, (62) 
%./N =1.03(V./N)}, (63) 
¢-=0.79P.(V./N)}. (64) 


These are based on the following approximate 
values estimated from the data and curves given 
by Lennard-Jones and Devonshire (in their 
notation except for T2, and T3,) : 


pvo/ NkT3-=0.36, PscSo/ NRT 2-=0.71 ? 
A/kT3,.=9, A,/RT 2 =8.6; (65) 
v-/Vo= 1.8, $-/So= 1.4. 
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It should be of considerable interest to find 
experimental values of the numerical constants” 
appearing in Eqs. (62)—(64) and in Eggs. (59)- 
(61). The question is, will these experimental 
constants, in suitable cases, be independent of the 
adsorbent and adsorbate? If, as there is some 
reason to expect, constant ‘‘constants”’ are indeed 
found, say, for small spherical non-polar ad- 
sorbates on non-porous non-polar adsorbents, the 
values of these constants should be useful in 
providing an additional check on theories of the 
liquid state and condensation. That is, if both 
two- and three-dimensional versions of a theory 
are worked out, as was done by Lennard-Jones 
and Devonshire, the theoretical constants of 
Eqs. (62)—(64) may be calculated and compared 
with experiment. 

It is planned to carry out in this laboratory 
some experimental work in connection with the 
theoretical problems discussed in Sections III, 
IV, and V. 

The author is indebted to Professor K. S. 
Pitzer, Dr. George Jura, and Dr. P. H. Emmett 
for their helpful comments. 


APPENDIX I 


We use the notation of Devonshire in this Appendix. 
For Devonshire’s model the configuration integral Q; is 


Ps i ee li $0 § | ‘dl 9) 
0;= 1 74 Ng exp py !.068 0.508) ; (69) 

with 
0=S)/s, (70) 


g ={" exp EB — (9) +20%m(y)} dy, (71) 


The @ defined above has approximately the same physical 
significance as the @ in Section III. Using a?=2s/v3N, 


*% The fact that constants are expected which are in- 
dependent of the gas (if the intermolecular potential func- 
tion is of the same form for different gases, and is identical 
for a given gas in two and three dimensions) follows from 
the law of corresponding states, especially as discussed by 
Fowler and Guggenheim (reference 3, pp. 342-345). 
Partly in their notation we have 


kT3-=const. X A*, 
V./N=const. Xr*®, 
P.=const. XA*/r®; (66) 
and also . 
kT2.=const. X A*, 
%./N =const. Xr*, 
¢c =const. X A*/r*?, (67) 


Combining these equations, 


T/T 3¢=const., 
%-/N =const. X (V-/N)3, 
¢e=const. X P.(V-/N)!. (68) 


L. HILL 


Eq. (69) becomes 


2a 


y y N 
0,={ 758 exp [ g71-060"—0.50%) | ; (72) 
This gives Q; as a function of s and N, After obtaining 


A lo : ‘ i . 
= meee, as required by Eq. (51), there results finally 
COL 


, WS 6. As f 4 Zn 
p ~ Daso* ¢ exp | — pp (424+ “) 


—0°(3.5 of) | (73) 
7 f / 


or, in a form which is more convenient for calculations, 


(= yea FNP | Fag t— ppt.060"—0.50%)] (74) 
The tables given by Devonshire for A,/RT =9 and 10 were 
used to calculate 27s9*p’/v3, as plotted in Fig. 4. 

The critical values of A,/kT and s/so are estimated by 
Devonshire to be approximately 8.6 and 1.4, respectively. 
Hence, @.=1/1.4 =0.71. 


APPENDIX II 


In studying an adsorbed insoluble surface phase of a 
single substance which is in equilibrium with the same 
substance in gaseous form, there are certain thermo- 
dynamic variables which are most useful: the surface 
area Y%, the equilibrium number of adsorbed molecules .V 
(in the area 1), the absolute temperature 7, the pressure p 
of the gas phase in equilibrium with the adsorbed phase, 
the two-dimensional (spreading) pressure ¢ of the adsorbed 
phase, and the Helmholtz free energy A of the adsorbed 
phase. %, NV, and 7’ may most conveniently be thought of 
as independent variables and ~, ¢, and A as dependent 
variables. Also, {, N, and A(&, N, T) are extensive proper- 
ties while T, p(X, N, T), and g(A, N, T) are intensive 
properties. Hence p and ¢ may be written p(I, 7) and 
¢(I', T) and the same is true of A/A, where f= NV /%, the 
surface concentration. 

We wish to consider the relations between the adsorption 
isotherm, which is obtained experimentally in the form 
p(l, T) or '(p, T), and the equation of state ¢(I, T). We 
assume for convenience that the gas phase is a perfect gas 
so that uo, the chemical potential of this phase, may be 
written 


ug=p(T)+kT log p. (75) 


In an exact treatment, we would replace p by the fugacity. 
The chemical potential wa of the adsorbed phase is, by 
definition, 


wa=(O0A/ON)M%,T. (76) 

At equilibrium, 
; MA=UG, (77) 
(0A /ON)4,7T =p°(T)+kT log p. (78) 


This equation relates A and p. The corresponding relation 
between A and ¢ is contained in the thermodynamic 
definition of ¢: 

g=—(0A/dA)y, 7. (79) 


The next step is to eliminate A from Eqs. (78) and (79). 
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Integrating Eq. (78) for fixed values of & and 7, In differential form, Eq. (85) is 


AQ, N, T)= f° [u(T)+ET log pILN, — (80) dg=VRT¢ log p. (86) 


since A(U, 0, T)=0. Therefore, using Eqs. (79) and (80), Equations (85) and (86) are the familiar Gibbs equations 
‘ ‘ : as applied to this type of adsorption. The derivation 

and keeping T constant in the following, < : : ; 
given here proves helpful in relating the thermodynamic 





-—y= [2 (wT) +kT log pin, (81) and statistical aspects of the problem. It will be noted 

vo amt from the derivation that Eqs. (85) and (86) are thermo- 

-erf, (° log ) AN. (32) | dynamically exact, except for the assumption that the gas 

. au Jy phase behaves like a perfect gas. These equations are also 

Now, since p= p(T), general, applying to monomolecular, multimolecular, mo- 
(2198 #) -(3) dlogp_ _V'dlogp (83) bile, and localized adsorption on a solid or liquid. Of 

Ou Jy Wn dT % ar ' course ¢ is not of much physical interest except for the case 


of mobile monomolecular adsorption. 
p The usual procedure in obtaining the equation of state 
<dT (84) from the experimentally determined adsorption isotherm 


Equation (82) becomes 


aT y 
g= k rf, pl log. 





ar Z - pai 
- I'(p) is to evaluate ¢(p) from Eq. (85), and then, by 
”? bining T _ lr) at the giv - 
e=bT | Pd log p. (85) combining I'(p) and ¢(p), find ¢g(f) at the given tem 
Jo perature. 
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Time-Dependence of Boundary Tensions of Solutions 


I. The Role of Diffusion in Time-Effects 


A. F. H. Warp anv L. Torpar 
Faculty of Technology, University of Manchester, Manchester, England 


(Received April 19, 1946) 


The variation of surface tension of solutions with time has in some cases, where the change is 
over within a few seconds or less, been explained on the basis of diffusion. This paper attempts a 
rigorous mathematical analysis of the role that diffusion might play in such time-effects. The 
limitations of diffusion theories which have been proposed previously are discussed. A general 
theory of diffusion to the surface is derived, which allows for back-diffusion and which makes no 
special assumptions of a physical nature. It is possible to use Fick’s equation to calculate the 
total amount of solute which diffuses from a semi-infinite bulk solution into the surface if the 
concentration immediately under the surface is known at various times throughout the process. 
It is shown how the latter information may be deduced from the variation of surface tension 
with time together with final equilibrium values of surface tension. The methods of this theory 
are applied to analyze recent data on time-effects of short duration. It is concluded that even in 
cases where the variation of the surface tension is over in less than a second the rate-determining 
process is not diffusion. Even for these very rapid changes one is therefore led to assume the 
existence of an activation barrier which determines the rate of adsorption. 


INTRODUCTION period of several days.2* On the other hand, 


ANY instances have been reported of the Variations of an apparently similar kind may be 
surface tension of solutions being a func-. complete in a few hours or within an hour.*~® In 


tion of time. In some cases this has been ascribed some cases the time-effect may last only for a 


to slow chemical changes in the bulk of the 2.N. K. Adam and H. L. Shute, Trans. Faraday Soc. 31, 


solution.! But many observed time-variations 205 (1935). 


. . . . aN « . . T 5. c ay “+ , 
cannot be explained in this way. The change in 38188 and Ht. £.. Thats, Tran Faany See, 28 


surface tension may occur very slowly, over a *K. S. G. Doss, Kolloid Zeits. 84, 138 (1938). 
enemies : ; 5K. S. G. Doss, Kolloid Zeits. 86, 205 (1939). 
‘A. Lottermoser and B. Baumgiirtel, Trans. Faraday °H. V. Tartar, V. Sivertz, and R. E,. Reitmeier, J. Am. 


Soc. 31, 200 (1935). Chem. Soc. 62, 2375 (1940). 
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few seconds to a few minutes.’ Finally, time- 
variations of surface tension have been observed 
which are complete within a fraction of a 
second.*-"' Quite recently, the authors have 
shown that interfacial tensions may also be 
time-dependent.” The same theory may not 
explain all these effects; in fact many different 
explanations have been proposed. A review of 
these explanations is not contemplated here, but 
merely an examination of the role that diffusion 
to the surface may play. 

The suggestion that diffusion of the solute to 
the surface was the cause of the alteration of 
surface tension with time was first made by 
Milner,'* who did not, however, treat the theory 
quantitatively. Equations to describe the phe- 
nomenon in terms of diffusion have been deduced 
by Bond and Puls,* Doss,® Ross,'* and Langmuir 
and Schaefer.'® Their treatments involve various 
approximations and simplifying assumptions. 
Although no generally exact quantitative treat- 
ment of diffusion to the surface has been put 
forward, it seems clear that many of the ob- 
served cases of alteration of boundary tension 
with time are very much too slow to be explained 
by diffusion. In general terms, there is some 
kind of activation barrier, which hinders the 
establishment of the adsorbed film.*!? For those 
cases where the alteration of surface tension with 
time is very rapid, diffusion theories so far put 
forward are insufficiently exact to determine 
critically whether the variation with time can be 
entirely attributed to diffusion. The object of 
the present paper is to develop a general theory 
of diffusion to the surface which will allow a 
definite decision to be made in any case as to 
whether diffusion is a complete explanation of 
the time-variation or whether some additional 
effect must be sought. 


7J. W. McBain and D. A. Wilson, J. Am. Chem. Soc. 
58, 379 (1936). 
8W. N. Bond and H. O. Puls, Phil. Mag. [7] 24, 864 
1937). 
9 o C. Addison, J. Chem. Soc., p. 252 (1944). 
10 C, C, Addison, J. Chem. Soc., p. 477 (1944). 
1 C, C. Addison, J. Chem. Soc., p. 98 (1945). 
12 A, F. H. Ward and L. Tordai, Nature 154, 146 
1944). 
; 13 $ R. Milner, Phil. Mag. [6] 13, 96 (1907). 
4S. Ross, J. Am. Chem. Soc. 67, 990 (1945). 
% J, Langmuir, and V. J. Schaefer, J. Am. Chem. Soc. 
59, 2400 (1937). 
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Fic. 1. Successive stages in concentration gradient from 
surface into the bulk solution during diffusion to the 
surface. (1) initial; (2), (3), and (4) intermediate; and 
(5) final stages. 


NATURE OF DIFFUSION PROCESS 


Before the diffusion theory is developed, it is 
convenient to give a physical description of the 
process. We may suppose that initially there is a 
solution of uniform concentration and at a 
certain instant a plane surface (solution /air) is 
created. At this moment the concentration of 
solute in the surface is the same as the bulk 
concentration. If the solute is a substance which 
is positively adsorbed, its molecules are in a 
position of lower potential energy in the surface 
than in the bulk. Although individual molecules 
move in all directions, the resultant flow of 
solute molecules is in a direction from the bulk 
to the surface. The region of the bulk solution, 
with a thickness of a few molecular diameters 
only, immediately next to the surface, is referred 
to here as the ‘“‘sub-surface.”’ At the initial stage 
in the formation of the adsorbed layer, solute 
molecules will diffuse from the bulk to the sub- 
surface more slowly than they pass from the 
sub-surface to the surface, on the assumption 
that equilibrium is established instantaneously 
between the latter two regions. This is so, be- 
cause initially the surface is practically empty 
and every molecule arriving at the surface is 
likely to arrive at an empty “‘site’’ and become 
adsorbed. For all practical purposes the sub- 
surface concentration may be regarded as falling 
from its initial value to zero immediately. 

Eventually, as the surface becomes fuller, there 
is an increased probability that a solute molecule 
in the sub-surface will arrive at a portion of 
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the surface already occupied by the adsorbed 
solute. Such a molecule will, of course, remain 
in the sub-surface, which will accumulate solute 
molecules more rapidly than these are passed on 
to the surface. As soon as the concentration in 
the sub-surface is sensibly different from zero we 
must consider the possibility of back-diffusion 
from the sub-surface to the bulk. The sub-surface 
concentration will continue to rise until eventu- 
ally, at equilibrium, it will have a concentration 
equal to that in the bulk (on the assumption 
that the bulk region can be regarded as an 
infinite reservoir of solute molecules). Successive 
stages in the concentration gradient are shown 
in Fig. 1. 

If the sub-surface concentration at each in- 
stant is known throughout this process, the 
diffusion of molecules from the bulk to the sub- 
surface might be treated by application of the 
ordinary diffusion equations. But in fact the sub- 
surface concentrations cannot be obtained di- 
rectly. The experimental measurements are the 
values of the surface tension at various times. 
Assumptions have to be made relating the surface 
tension to the sub-surface concentration. It is 
partly on account of the approximate nature of 
these assumptions that previous treatments of 
this problem have been inadequate. 

Doss**® avoided this difficulty in his study of 
the slow adsorption of benzopurpurin at an air/ 
water surface by measuring directly the surface 
concentration at various times by compressing 
the surface film and measuring the area under a 
definite pressure. However, his theoretical treat- 
ment was inaccurate. He neglected back-diffusion 
from the surface to the bulk. His calculations 
were based on the Einstein diffusion equation for 
the mean displacement of a particle, but in the 
application of this he assumed that the molecules 
reaching the surface all started from a region 
where the concentration equalled the bulk con- 
centration. This would only be true at the very 
beginning and at the end of the diffusion process, 
and the assumption is equivalent to neglecting 
the concentration gradient near the sub-surface. 
This same fundamental inaccuracy is also to be 
found in the treatment by Ross,’ who follows 
Doss’s method very closely. Ross compensates 
for certain deficiencies in the diffusion theory 
by supposing that the adsorption layer is multi- 
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molecular. The hypothesis of thick layers on the 
surface of solutions containing such simple 
solutes is, however, not acceptable. It has been 
shown recently'® that the assumption of mono- 
layers of the solute on the surface of these solu- 
tions satisfies the thermodynamical requirement 
that the activities of the solute in the bulk and 
surface phases, respectively, shall be propor- 
tional to each other. It would therefore appear 
unjustifiable to postulate the existence of multi- 
molecular layers to explain this new effect. 

A diffusion theory was proposed by Bond and 
Puls® based on Fick’s diffusion equation. Since 
this is a statistical equation based on the behavior 
of all the molecules, it should allow a general 
solution to be derived, and if appropriate bound- 
ary conditions are employed, it should cover 
both forward- and _ back-diffusion. However, 
Bond and Puls made the limiting assumptions 
that the adsorption layer could be supposed to 
be almost unoccupied and that the lowering of 
surface tension was proportional to the concen- 
tration. This restricts their treatment to the 
initial stages only of the diffusion process. In 
addition, their final integrated equation 


(¢—9)/q=exp [—(2/n')(t/7)*], 


(their Eq. (9)), which has been quoted and used 
by other workers, is in fact incorrect, presumably 
because of an erroneous integration procedure. 


SIMPLE DIFFUSION THEORY 


A simple theory may be developed first for the 
early stages in the diffusion when the backward 
movement of solute from the sub-surface to the 
bulk may be neglected. The surface is then 
sufficiently empty to take up solute molecules as 
fast as they arrive and the sub-surface concentra- 
tion is practically zero. The problem is analogous 
to that of the diffusion of heat from an infinite 
conductor, initially at uniform temperature, into 
a conductor of infinite conductivity placed in 
contact. 

Let n=concentration in the solution at a 
distance x from the surface at a time ¢t. From the 
equations for heat flow in the analogous case!’ 


16 A. F. H. Ward and L. Tordai, Trans. Faraday Soc. 
(in press). 

17H. S. Carslaw, Mathematical Theory of Heat Conduction 
(MacMillan and Company, London, 1921), p. 34. 
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we obtain the equation 


on No x? 
S-Sao(-i) 
dx (Dt)? 4Dt 


where mp is the constant bulk concentration and 
D is the diffusion coefficient defined by Fick’s 
equation 


(1) 





dn/dt= Dd?n/dx?. (2) 


At the extreme edge of the bulk solution, i.e., the 
sub-surface, where x =0, Eq. (1) is simplified to 


(0n/0X) 2-9 = No/(rDt)', (3) 


and the rate of diffusion across 1 cm? of the plane, 
x=0, is given by 


(9M/dt)z-0= —Dno/(xDt)'. (4) 


Hence, the amount of solute which has diffused 
across this plane into the surface up to the time, 
t, i.e., the surface concentration at that time, J/, 
is given by 


OM Dix? 
M= [GC 7% -di=an(—). 
T 


sThis equation is identical with one used by 
Langmuir and Schaefer! to explain diffusion of 
ions into monolayers. 


(5) 


GENERAL DIFFUSION THEORY 


It is also possible to develop a more general 
theory which will make allowances for back- 
diffusion. Such a theory may be derived from 
Eq. (2) and the treatment is applied to a single 
phase extending from the bulk (assumed infinite 
in extent) to the sub-surface, but not including 
the surface. No assumptions regarding the rela- 
tion between surface and sub-surface concentra- 
tions are introduced at this stage and the 
solution is therefore quite general. 

Consider linear diffusion of solute molecules 
towards the boundary plane x=0 of a semi- 
infinite liquid solution. The bulk concentration 
remains effectively constant (m9) during the 
process. The concentration in the sub-surface 
varies with time and may be expressed as ¢(t). 
To obtain a solution of the problem it is possible 
to utilize a mathematically analogous theorem 
in heat conduction, based on Eq. (2) and the 
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use of Green’s functions!? (p. 173). The tem- 
perature v at a distance x from the surface of 
the semi-infinite solid of thermal conductivity x, 
initially at temperature vp with its surface at 
temperature ¢(/) at time ¢ is given by 

















Vo S| r (ee y% 
y= —— exp | -———— 
2(axt)} | : 4xt J 
(s+2" x)" | 
—exp| — -| pdx 
L 4xt J 
re 
exp | —————— 
" 4x(t—r) 
—f o(r) aT, (6) 
Peny (t—7)} 


where 7(7<?#) is the time at which an instan- 
taneous point source of unit strength, at a 
distance x’ from the surface, generates heat. 
It is important to note that Eq. (6) is made up 
of two terms. The first term corresponds to a 
flow of heat from the bulk region (which is 
initially at temp. vo) towards the surface (initially 
at zero temp.). The second term represents a 
flow of heat in the opposite direction, from the 
surface at temp. ¢(¢) towards the bulk (initially 
at zero temp.). In the diffusion nomenclature, 
we are interested in obtaining 


(dn /OxX)2-+0, 


where , the concentration, replaces v; m) and D 
replace v and x. D is the diffusion coefficient. 
To simplify the problem we employ two sets of 


transformations. In the first term we may write 


x’=x+2(Dt)'é and x’=—x+2(D0)}é, 


for the first and second integrals, respectively. 
In the second term 


p=x/2(D(t—r))', thatis, t—r=x?/4Dyz’. 


With such a change of variables Eq. (6) re- 


duces to 
Ino z/2(Dt)? 

n=—[ exp (— edt 
Tv 0 


x? 


2 x 
+5 fo o(t- 
mw J2/2(p1)) 4D? 





(7) 


) exp (—y?)du 
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It will be convenient to treat the two terms in 
Eq. (7) separately. It will be assumed that the 
first term contributes 0”,/dx and the second 
dn2/dx to give the total dn/dx. Differentiation 
of the first term with respect to x gives 


On, No a 
et xp (- — J. (8) 
dx (xDt)} 4Dt 


If x tends to zero in Eq. (8), we obtain Eq. (3) 
and for the amount of solute, 14;, contributed 
by this direction of flow, an expression is ob- 
tained which is identical with Eq. (5). 

Turning our attention to the second term in 
Eq. (7) and writing a=x/2(Dt)', we obtain 


1 9 
2 


1? ' a? 
n= | (-—1) exp (—y?)du. (9) 
2 ue 


a 


The integrand in Eq. (9) may be expanded by 
Maclaurin’s theorem to give 


Ri — 


=> are) { exp (—u*)u- "dz, 
0 0 


and by repeated integration by parts we obtain 
the series 


7 F 
n= 40) [ exp (—p”)du 
2 a 


ra (— 1)" 
+X a——1"9(t) 
1 


r: 


1 2a” 
et a 
(2r—1) ' 2r—1)(2r—3) 





2 
| 
(2r—1)(2r —3)(2r—5) 
Differentiating with respect to a, we have 


T? ON» 

——=-— exp (—a’)¢(t) 

2 da 
D (—1)’ 
> ——"'o (t) + O(a’), 
1 ri(2r—1) 


and letting a—0 as x0 the final result is 


On, 1 2 (—1)'t" 
( ) a p g(t), 
Ox F <0 (rDt)} 0 (2r—1)r! 





the first few terms of this series being 


1 
- ~) b(t) +i¢’ (t) 
(xDt)} 


Zz. 
“4 ‘o"(O+- 16") 


It is more convenient to express the series (12) 
as an integral. It is easily shown that 


On» 1 * (2) 
( ) =— -f dz, (13) 
Ox J 2-0 2(xD)} 0 (t otis z)! 


where z is a variable; for putting (t—z)=@ we 
obtain an integral which can be expanded by 
Maclaurin’s theorem to give Eq. (12). 

The amount of solute, 1/2, which flows from 
the surface in the direction of the bulk, must 
next be determined. Starting from Fick’s equa- 


tion 
(0M2/dt)z-0= —D(dn2/dx) 2-0, 


we are led to consider the integral 


(2) . 
{,.=- (- SS dzdt. (14) 
(t—z)3 


It can be shown that Eq. (14) reduces to 


Dr»? ¢! o(z) 
Me= —) f —d 
0 (t—z)} 


To prove the equivalence of Eqs. (14) and (15) 
use will be made of Cauchy’s theorem on n-times 
repeated integrals. This theorem, together with a 
process involving differentiation and integration 
to fractional orders enables us to complete the 


proof. 
—-~ — fu-s )"-16(s)dz, 
“Pp (n) 

then 


If we put 
rs n) 
‘(t—z)"" 
“y ———-9(z)dz, (17) 
o (n—1)! 


and Cauchy’s theorem states 


t (t—z)"- -1 
J ——— $(z)dz= D~"9(t). 
0 (n—1 )! ; 
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The operator, D-, stands for fo’: --dz, whereas 
D' has its usual meaning d/dz. For our particular 
case we have 
* $(2) 
F(t)= | ———dz. 


0 (t—2)? 

From Eqs. (16)—(18) we obtain 
F(t)=1(2)D“9(0). (19) 
We are interested in finding the derivative, F’(¢), 


and to show that it gives an integral; such as 
occurs in Eq. (13). We can write 


 @ ¢* of) 
F'(t)=— — a 
dt 0 (t—2z)} 


=1'(3)DD~“‘$(t) =1'(2)D'o(t). (20) 


° 





From Eggs. (17), (18), and (20) we have 





1 
PW) =13)— 
r(—4) 


J (t—2)-19(2)dz 
0 


f o) dz, (21) 
0 (¢—z)! P 





Nin 


which is the result required to prove the equiva- 
lence of Eqs. (14) and (15). 
For the two opposing processes we now have 


Dv} ‘ 
u-2(—) {not — f #(e,dL(—2)}, (22) 


where the expression for M2 is written in a 
slightly different form which has advantages 
from the point of view of numerical computation. 

It is possible to derive the expression for M, in 
a very much simpler but less rigorous way. 
Consider a semi-infinite liquid solution with zero 
initial concentration, the surface concentration 
being kept constant at a value c. After the lapse 
of time ¢, the amount which has diffused inward 
from the surface (taken as the positive direction) 


is given by 
Div} . 
u’=2¢(—) ’ (23) 
us 


From this system we can reproduce another with 
a varying surface concentration ¢(t). It isassumed 
that the surface concentration is changed in a 


stepwise manner so that Eq. (23) still applies 
for a short time, during which the surface 
concentration is constant and yet the stepwise 
change is taken to approximate ¢(t), the latter 
being given by a smooth curve. Such an argu- 
ment assumes, justifiably, that each portion of 
the solute which diffuses from the surface during 
one of the periods with a constant surface con- 
centration, Cn, continues its diffusion independ- 
ently of any other portions of the solute. This 
is equivalent to the assumption that the diffusion 
is an ideal process, that is to say that the diffusion 
coefficient is independent of concentration. 
We have then 


Aco on the surface during diffusion from t=0 to t=T, 
Ac; on the surface during diffusion from t=1 to t=T, 


Ac, on the surface during diffusion from t=n to t=T, 


where 7 = time needed to reach equilibrium con- 
ditions between bulk and surface. The total 
diffusion inward from the surface in the interval 
between ¢=0 and t= T is obtained by summation 


Dvir 
E u'=2(=) > Ac,(T—2)}. (24) 
T 0 


We can make the approximation as good as we 
please by decreasing the duration of each step. 
In the limit, the expression (24) can best be 
represented by plotting ¢(t) against (J —?)! and 
multiplying the area under the curve by 2(D/rz)?. 
This is the integral 


Dv} p™ 
2(—) f o(t)dl (T—t)5], (25) 
Tv 0 


which is identical in form with the corresponding 
part in Eq. (22) except for the sign. In Eq. (22) 
Mz is naturally negative as it represents diffusion 
away from the surface taken in a negative 
direction. 


CALCULATION OF DIFFUSION COEFFICIENTS 


The Eq. (22) enables us in theory to obtain 
either M or D if one of the two is known. Un- 
fortunately, explicit integration in Eq. (22) is 
not possible as ¢(¢) is not known explicitly. We 
can, however, employ a physical artifice which 
makes graphical integration possible and gives M 
the surface concentration at time ¢. 
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If it is assumed that there is no activation 
barrier between the sub-surface and the surface 
and that diffusion is the sole reason for the delay 
in establishing adsorption equilibrium, then the 
time required for the solute to diffuse from the 
sub-surface to the surface is very much smaller 
than the time required for diffusion equilibrium 
to be established between bulk and sub-surface. 
At any instant, there is a state of equilibrium 
between the surface and the sub-surface with 
respect to the solute. 

The observed boundary tension is determined 
by the concentration of solute in the surface and 
this is in turn determined by the concentration 
in the sub-surface, ¢(t). If a certain value of the 
boundary tension (y) corresponds to a certain 
value of #(t), say c, at any moment during the 
diffusion, that particular value of y would also 
be obtained under equilibrium conditions when 
the concentrations in the sub-surface and the 
bulk were both equal to c. The relationship 
between y and ¢(¢) can therefore be obtained 
from observations made under equilibrium con- 
ditions when diffusion has ceased. Values of the 
equilibrium boundary tension are plotted against 
the corresponding concentrations. In this way, 
from observed values of y during the diffusion 
process, the variation of ¢(¢) with time may be 
found without the necessity of making arbitrary 
assumptions about the relationship between y 
and the surface concentration and between the 
latter and ¢(¢). In this way the integral for M» 
can be evaluated in terms of D!, for the period 
of time which is necessary for equilibrium to be 
established. The values of 1 (=Gibbs’ T) at 
equilibrium may be obtained independently by 
application of the Gibbs’ equation to the ob- 
served relationship between the concentration 
and the equilibrium boundary tension. Since M 
is known in units of D}, it is then possible to 
calculate D. 

This method makes it possible to obtain diffu- 
sion coefficients from measurements of boundary 
tensions which vary with time. These will, of 
course, only be equivalent to the conventional 
diffusion coefficients if there is nothing in the 
nature of the physical system to invalidate the 
assumption of instantaneous equilibrium be- 
tween the surface and the sub-surface. If this 


assumption is not true, and there is an activation 
barrier between the surface and the sub-surface, 
the value of D calculated from surface tension 
measurements (which will be referred to as D,) 
will be a complex quantity characteristic of both 
diffusion and the crossing of the barrier. In such 
a case either of the processes may assume a 
controlling role. For a calculated value of D, 
which is of the same order of magnitude as the 
conventional D, the normal diffusion process is 
the slower of the two. In such a case the existence 
of a small additional barrier at the surface would 
become a trivial question, as no possible experi- 
ment could be devised to investigate the nature 
of this barrier. For a calculated value of D, 
which is smaller than the conventional D it is 
evidently the barrier at the surface which has 
the controlling role. If D, is very small, the 
diffusion process can be neglected in comparison 
with the activated adsorption. 

Quantitative conclusions cannot, however, be 
drawn from such very small values of D,, for 
it must be borne in mind that in such a case the 
assumption of instantaneous equilibrium between 
surface and sub-surface is invalidated and no 
great significance can be attached to the calcu- 
lated value of D,. The theory should therefore 
give the correct value of the diffusion coefficient 
(as compared with the conventional value) if the 
process is purely a diffusion process, but if the 
calculated value of D, is smaller than the usual 
bulk diffusion coefficient, it indicates that diffu- 
sion is not the rate-determining process. 

A series of short time-effects have been in- 
vestigated by Addison,*"!' who measured the 
surface tensions of freshly formed surfaces of 
aqueous solutions of fatty alcohols, at intervals 
of a small fraction of a second. He employed the 
oscillating jet method and measured changes of 
surface tension with time, up to about one 
second, by which time equilibriam had been 
established after a rapid initial fall. Experimental 
results were reported for n-amyl, iso-amyl, 
n-hexyl, n-heptyl, and n-octyl alcohols at 20°C, 
some of them at several concentrations. The 
theory developed in this paper was applied to 
Addison’s data and values of D, were obtained 
by the graphical integration method. Calculated 
values of D, are collected in Table I (column 6). 
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Fic. 2. Variation of calculated and observed diffusion 
coefficients with chain-length. (1) Observed D. (2) Calcu- 
lated D,. 


DISCUSSION 


Values of I have been obtained by the use of 
the approximate Gibbs’ equation (using con- 
centrations) since activity values for solutions 
of these alcohols were not available. For the lower 
alcohols at the lowest concentrations this was a 
sufficiently good approximation. It may be seen 
from Table I that for the highest concentration 
of n-hexyl alcohol the calculated value of IT 
(8.67 X10-") is obviously too high. This is 
probably caused by the divergence between 
concentration and activity. Since a high value 
of I yields a high value of D,, it appears that if 
accurate activity data were used the values of D, 
in column 6 would at any rate be no higher than 
those given. They may therefore be taken safely 
as an upper limit. 

The conventional diffusion coefficient (D) for 
iso-amyl alcohol in water has been obtained by 
Thovert!’ as 6.9 X 10~° cm? sec.~! at 15°C, corre- 
sponding to 7.9X10-® at 20°C. The calculated 
value of D, in Table I is 56 times smaller than 
this. Therefore diffusion is not the rate-deter- 
mining process here, and a further activation 
barrier must exist. Experimental values of diffu- 
sion coefficients are not available for the higher 
alcohols. But the probable values may be seen 
in Fig. 2, where the curve of experimental values 
of D for the lower alcohols (curve 1) is extrapo- 
lated. The decrease in D with increase in the 


18 J. Thovert, Ann. d. Physik [9] 2, 369 (1914). 
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size of the molecule would be expected from the 
Stokes-Einstein formula 


D=kT/6xrn. 


Curve 2 shows the values of D, calculated from 
Addison’s data for the lowest concentrations in 
each case. Up to hepty! alcohol the values of D, 
are all lower than the correct values of D. But 
the most convincing argument against the process 
being merely diffusion is not so much the fact 
that D,<D as the general trend of D, to increase 
with chain-length. This is the opposite of what 
would be expected for a diffusion coefficient. 
The fact that D, for octyl alcohol is of the right 
order of magnitude for agreement with D is 
therefore probably fortuitous. Thus the numeri- 
cal results indicate that the rate of adsorption in 
these cases is not determined by diffusion. 

Some doubt is thrown on the experimental 
values for the higher alcohols considered in 
Table I because of the probability that equi- 
librium was not attained. If $(¢) is plotted against 
the time, the curve should flatten when equi- 
librium is reached. These values are plotted in 
Fig. 3, for m-amyl alcohol (0.04544 molal) and 
n-octyl alcohol (0.0007962 molal). 

It will be observed that the curve for n-amy] 
alcohol flattens out satisfactorily within the time 
of observation. The curve for n-octyl alcohol 
indicates that equilibrium may not have been 
attained within 0.11 sec., the time of the final 
observation. 

Available literature seems to indicate that the 
time required to obtain equilibrium surface 
tensions increases very rapidly with chain-length 
for solutions of aliphatic acids. King,'® using a 
drop-weight method, found that drop-lives of at 




















TABLE I. 
(1) (2) (3) (4) (5) (6) 
Solute M:XDi M:XDt Yr x10 D, x10" 
(alcohols) Molality 106 X106 mole/cm? cm?/sec. 
n-amyl 0.01136 2.026 1.348 2.93 0.185 
n-amy! 0.02272 3.242 2.445 3.98 0.250 
n-amy! 0.04544 5.612 3.750 5.92 0.100 
iso-amyl 0.01384 2.410 1.620 2.97 0.141 
n-hexyl 0.00245 0.895 0.707 2.09 1.28 
n-hexy! 0.00735 2.620 2.290 3.77 1.31 
n-hexy| 0.01226 4.360 4.070 5.44 3.62 
n-hexy] 0.02161 5.975 5.240 8.67 1.39 
n-heptyl 0.001724 0.645 0.497 2.83 3.66 
n-octyl 0.0007962 0.298 0.127 4,98 8.48 








19H, H. King, Kansas State Agr. Col., Agr. Exp. 
Station, Tech. Bull., No. 9 (1922). 
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TIME-DEPENDENCE OF BOUNDARY TENSIONS 


least five minutes were required to attain steady 
values of surface tension with all acids above 
n-butyric acid. For solutions of nonylic acid 
30-minute drops had to be used and for decylic 
acid, the surface tension was still falling after 
30 minutes. Addison,” working with solutions of 
decylic acid, found that times of the order of a 
minute were required to reach a ‘dynamic equi- 
librium’”’ and that this was followed by a further 
slow fall of surface tension. These observations 
appear to throw more doubt on the existence of a 
real equilibrium, for solutions of long chain 
alcohols, after the lapse of a tenth of a second or 
less. Further experimental information on this 
point would be valuable. 

To explain the apparent increase in “migra- 
tional velocity’’ of the solute molecule with 
increasing chain-length, Addison has defined 
this velocity in quantitative terms and assumed 
that the main factors controlling the rate of 
adsorption were the free energy change on ad- 
sorption and the distance of the solute molecule 
from the surface. While it is not at all obvious 
why the velocity of a solute molecule traveling 
towards the surface should depend on its dis- 
tance from it, it is even less acceptable to sup- 
pose that this velocity should depend on the 
lowering of the surface tension by the adsorption 
of the solute molecule. In any physical or 
chemical change a system may undergo, the 
velocity of the change depends not on the differ- 
ence in the free energies between the initial and 
final states, but on any intermediate activation 
barriers. The numerical calculations of Addison 
are further invalidated by his use of the Gibbs’ 
equation for calculating surface excesses at stages 
prior to the establishment of equilibrium. 

The authors have been carrying out experi- 
ments on the time-dependence of interfacial 
tension between water and solutions of long 
chain compounds in hexane.” Full details of this 


20 C. C, Addison, Nature 156, 600 (1945). 
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work will be published shortly, but an example 
may be given now to support the conclusions 
reached above. For solutions of palmitic acid in 
hexane values of ¢(¢) were obtained from a 
complete set of equilibrium interfacial tension 
data. The surface concentration was obtained by 
means of an equation of state proposed by 
Langmuir.” For a concentration of 0.04101 molal 
the calculated value of D, so obtained was 
5.66 X10-!*, a value about 10’ times smaller 
than the expected magnitude of the diffusion 
coefficient. It is clear that in this system also 
there is an activation barrier between the sub- 
surface and the surface. 

We are much indebted to Dr. J. Howlett, Mr. 
W. K. Burton, and Mr. A. Howarth for several 
helpful discussions. We also wish to thank 
Imperial Chemical Industries, Ltd. (Dyestuffs 
Division) for a grant to one of us (L.T.) 
during part of the time in which this work was 
carried out. 


217. Langmuir, J. Chem. Phys. 1, 756 (1933). 
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Jones-Ray Effect, Wettability, and Zeta- 
Potential 


Hans M. CASSEL 
Speedry Products, Inc., New York 6, New York 
May 28, 1946 







XCEPT for an observation of Heydweiller! who no- 
ticed that chlorate solutions wetted the glass he used 
poorly, no consideration seems to have been given to the 
possibility that the presence of ions in solution may affect 
the angle of contact. In interpreting the Jones-Ray effect,? 
so far, perfect wettability of the capillary walls has been 
taken for granted. The explanation suggested by Lang- 
muir’ appears to be in accord with a great deal of experi- 
mental evidence.* However, it is difficult to reconcile his 
hypothesis with the Laplace equilibrium condition at the 
1 meniscus/film boundary in the limiting case of absolutely 
pure liquids where no difference of osmotic pressures. can 
be thought of. An alternative view is here presented, 
namely that an increase of the wetting angle may be in- 
volved in the apparent decrease of surface tensions as 
found with minute concentrations of strong electrolytes. 
The cylindrical capillary layer adjacent to a perfectly 
wetting meniscus is usually treated as if it were a ‘‘duplex 
film’’ of liquid in bulk. However, one should be aware of 
the fact that its vapor pressure is smaller than that of a 
hollow cylinder of liquid in mass of equal diameter whose 
vapor pressure equals that of a capillary film of twice the 
capillary radius. Therefore, the capillary film should be 
recognized as an “adsorbed phase.’’® Accordingly, one 
should not conclude that its ‘‘surface tension’’ is equal to 
that of the liquid, or that its “interfacial tension’’ at the 
film/solid boundary is the same as that of the liquid/solid 
interface. It appears also questionable that the solubilities 
of electrolytes or non-ionic substances are the same in the 
capillary film as in the liquid, because the dielectric con- 
stants may be different, that of the film probably being 
smaller. 

The following is based upon Gibbs’ treatment of ‘“Sur- 
faces of Discontinuity between Solids and Fluids’’® which 
implies that the ‘‘surfaces of tension”’ coincide closely with 
the physical surfaces of the solids, no matter what the 
thickness of the adsorbed films, provided that the radii of 
curvature are very large in proportion to the film thickness. 

According to Neumann’s equation the “apparent sur- 
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face tension’’ of the liquid equals the difference of the vapor 
solid surface tension, o’, and the liquid/solid interfacial 
tension, «. Whatever the value of the contact angle, a 
change in capillary height will result from the adsorption 
of a third component and the desorption of solvent by 
which the difference of the tensions is altered: 


d(o’ —o) = —(T2’ —V2)dy2— (T'3’ —T's)dys. (1) 


Here, the negative quantity of solvent adsorbed at the 
liquid/solid interface can be neglected in comparison to 
the excess adsorbed in the film. The excess of the third 
component in the film comprises the quantity adsorbed 
at the film/solid interface and the total amount dissolved 
in the film volume per unit area: 


(¢3'/c2')T 2’. (2) 


The molar potentials of solvent and solute are (approxi- 
mately) correlated by the (Duhem) equation: 


Codu2t+csdu;=0, (3) 


where ¢2 and ¢; are the numbers of moles per unit volume. 
Equation (1), therefore, assumes the form: 


C3 «C3 
d(a’—«) = jre(2-2) -ry+1.| dus, (4) 


which for very small concentrations reduces further to: 
d(o’ —o) =(T3—T's’) dys. (5) 


Hence, in the range of minute concentrations, an increase 
or decrease in capillary height would solely depend on 
whether the excess of adsorbed solute is greater at the 
liquid/solid interface or at the film/solid interface. 

Regarding the possibility that the dielectric constant 
of the film is smaller than that of the solvent in mass, the 
behavior of sugar solutions would appear as quite normal. 
In case of electrolytes, at first view, the same might be 
expected. However, Langmuir’s’? theory of the potential 
distribution between parallel solid surfaces as a function 
of their distance seems to open an understanding of the 
reverse Case. 

There is also experimental evidence in favor of the 
present interpretation. The diffuse double layer in the 
vicinity of an insulator of high Z-potential is obviously 
of the same character as that of a highly polarized metal 
electrode in contact with an electrolyte solution. Now, the 
contact angles of gas bubbles attached to metal electrodes 
have been the subject of extensive studies by Frumkin and 
his co-workers* who have established the fact that the 
contact angle decreases when the potential increases. 
Accordingly, it is logical to conclude that the Jones-Ray 
effect is caused by an increase of the wetting angle result- 
ing from a decrease in the Z-potential at the solution/solid 
interface. 


1A, Heydweiller, Ann. d. Physik (4) 33, 148 (1910). 

2G. Jones and W. A. Ray, J. Am. Chem. Soc. 59, 187 (1937). 

31. Langmuir, Science 88, 430 (1938). 

4L. A. Wood and L. B. Robinson, J. Chem. Phys. 14, 258 (1946). _ 
( 5 H. M. Cassel, J. Chem. Phys. 12, 115 (1944); J. Phys. Chem. 48, 195 

1944), 

6 J. W. Gibbs, Collected Works 1, 328 (1928). 

71. Langmuir, J. Chem. Phys. 6, 891 (1938). 

8 A. Frumkin, Couche Double, Electrocapillarité, Surtension (Hermann 
and Cie., Paris, 1936), p. 17. 
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Note on the Rigid Multi-Layer Postulate 
of Eversole and Lahr 


G. A. H. ELTON 
Physical Chemistry Laboratory, Battersea Polytechnic, London, England 
June 6, 1946 


VERSOLE and Lahr!’ have interpreted certain electro- 

kinetic measurements as indicating properties of rigid- 
ity in part of the double layer at solid-electrolyte interfaces, 
a view held by von Buzagh? and others. They have used! 
published zeta-potential figures to deduce and thicknesses 
of from 3A to 63A for the rigid layer in aqueous salt solu- 
tions, values which could be accommodated by classical 
double layer theory, although Bikerman’ is of the opinion 
that any rigid layer is caused by the inherent roughness 
of all solid surfaces. In a later experimental paper,* Ever- 
sole and Lahr produce evidence which, it is claimed, 
substantiates the previous calculations and indicates a 
rigid layer of thickness about 100A for water against 
quartz. 

In the experiments of Eversole and Lahr the “distance 
of separation,” ‘air, between a quartz lens and plate 
placed together in air was determined interferometrically. 
Water was then condensed between the lens and plate 
and the new “‘distance of separation,”’ ‘water, determined. 
Negative values both for ‘air and ‘water were obtained, 
but it was assumed that the subtraction of ‘air from ‘water 
eliminated the error to which the negative values were due, 
and gave the thickness of the rigid layer. The validity of 
this assumption is to be questioned, in view of the very 
large variations in the individual values of ‘air and ‘water. 
The results obtained by Eversole and Lahr (five in num- 
ber) are given in Table I. 

Although the deviation from the mean in each series is 
as much as 600A, the conclusion drawn is that the total 
layer thickness, ‘‘subject of course to considerable experi- 
mental error,’ is about 200A. Since this involves two 
interfaces, the layer thickness at one surface is taken as 
100A. 

If now Students “‘t-test’”’ for significance of differences 
between means of small populations, as described by 
Evans® is applied to the above results, there is indicated 
approximately a 75 percent probability that the difference 
between the means of the two series in Table I is not 
significant, because of the wide variation within each series. 

Using an interferometric principle similar to that of 
Eversole and Lahr and a more developed technique, details 
of which will be published later, the writer has found no 
evidence for the presence of these layers in water, methyl 
alcohol or toluene. Some of the results obtained are given 


TABLE I. 








“Layer thickness” 


fair A twater A (‘water —‘air) A 





—1210 —870 330 

— 950 —800 160 

— 90 +140 230 

—1160 —840 320 

S) — 350 —220 130 
Mean — 752 234 
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in Table II. (In the last column of this table, P indicates 
the approximate probability of insignificance in the differ- 
ence of the means, as determined by Students ‘‘t-test’’.) 


TABLE II. 








Liquid “xpt. fair A liquid A (‘liquid —‘air) A 
20 15 
- 5 5 
Water d — 5 0 
10 
0 
4 





5 

25 

Methyl 7 10 
alcohol —10 
—10 

4 


0 

Toluene 0 
15 

15 

4 








A further investigation of the subject is in progress, 
using air bubbles against quartz after the manner of 
Derjaguin,® who reports the presence of thick. interfacial 
films on glass and mica in water. In view of the recent 
work of Tolansky,’ any interferometric work on “smooth” 
mica cleavage planes must be viewed with suspicion, since 
variations in surface level of the order of several thousand 
angstroms may exist. 


1W. G. Eversole and P. H. Lahr, J. Chem. Phys. 9, 530 (1941). 

2 K.g., A. von Buzagh, Colloid Systems (Technical Press Ltd., London, 
1937), pp. 153-4. 

3 J. J. Bikerman, J. Chem, Phys. 9, 880 (1941). 

4W. G, Eversole and P. H. Lahr, J. Chem. Phys. 9. 686 (1941). 

5U. R. Evans, Chem. Ind. 14, 106 (1945); see also L. H. C. Tippell, 
Methods “ Statistics (Williams and Norgatt Ltd., London, 1941), 
pp. 112-3. 

6 B. Derjaguin and M. Kussakov, Acta Physicochimica 10, 25 (1939); 
B. Derjaguin and M. Kussakov, Acta Physicochimica 10, 153 (1939). 

7 E.g., S. Tolansky, Proc. Roy. Soc. A184, 41 (1945); S. Tolansky, 
Proc. Roy. Soc. A184, 51 (1945). 





Bond Hybridization in the Non-Tetrahedral 
Carbon Atom 


Joun E. KILPATRICK! AND RALPH SPITZER? 
Department of Chemistry, University of California, Berkeley, California 
June 17, 1946 


NASMUCH as Duffey’s calculations’ lead to the pre- 
diction of 180° H—C—H angles in cyclobutane and 

cyclopropane (and possibly ethylene) which we do not 
think reasonable, we are presenting a somewhat similar 
calculation which we had made and which gives more 
reasonable results. 

Duffey calculates s, p hybrids with maxima in the direc- 
tions of the required bonds. However, considerably stronger 
bonds can be formed by dropping this requirement and 
assuming that bond strength is proportional to the value 
of the orbital in the desired direction, whether or not the 
maximum of the orbital lies in this direction. 

The orbitals f; and fz were used to form the two carbon- 
carbon bonds whose angle is determined by the geometry 
of the molecule. The directions of the other two bonds are 
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taken in directions of the maxima of f; and fs. With no 
loss in generality and a considerable gain in the simplicity 
of the equations, orbitals f; and fz can be chosen with 
their maxima in the xz plane, equally inclined to the z axis. 
By symmetry, f; and f; then lie in the yz plane, also equally 
inclined to the z axis. 

Applying the usual ‘normalization and orthogonality 
conditions we have: 

fi=as+bp.+(1/2)'p., 
fz=as+bp.—(1/2)'p., 
fs=bs —ap.+(1/2)%,, 
fs=bs —ap.z—(1/2)4p,, 
b=+(4—a2)}. 

We now determine the value of a (and b by Eq. (5)) 
so that f; will have the maximum value in a particular 
direction @) in the xz plane by setting df,;/da=0. This 
requires that 


(6) 


1 ) 
wt 2+6 cos? 0)/ ’ 


3cos?@ \t a 
b=| —— . (7) 
2+6 cos? 6 


The angle 63, between the bonds formed with the 
orbitals f; and f, is given by the expression: 
—2(1+3 cos? 6o)4 


3 cos? A 





(8) 


tan 634= 


We have calculated the H—C—H angle, 63s, and bond 
strengths, fo: and fos, for several values of 269, the C—C—C 
angle (Table I). Assuming bond energies proportional to 


TABLE I. 





Strength of Strength of 
bonds (1) and (2) bonds (3) and (4) 


for O34 fos 


1.414 126° 52’ 1.982 
1.887 121° 58’ 1.989 
1.984 115° 23’ 1.995 
2.000 109° 58’ 2.000 
2.000 109° 28’ 2.000 
1.996 105° 50’ 1.998 

1.732 








square of bond strength and using! 59 and 87 kcal. for 
normal C—C and C—H bond energies, we obtain 10 and 
2 kcal. as the strain energies per CH2 group of cyclopropane 
and cyclobutane in fair agreement with experimental 
values of about 10 and 4° (estimated from cyclobutane 
derivatives). Probably about 2 kcal. should be added to 
the calculated values for torsional strain.6 The energy 
calculation does not apply to ethylene because the bond 
is much shorter than it is in normal paraffins. 

That the exterior angle in cyclopropane is not 180° as 
suggested by Duffey is shown conclusively by the large 
dipole moment of 1,1-dichlorocyclopropane.? Also, the 
calculation Duffey quotes from Herzberg’ is incorrect. 
The correct C—C distance for 180° angle is 1.616A not 
1.556A. The moment of inertia of cyclopropane is fitted, 


for example, by rce-c=1.530A (electron diffraction), 
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rc-H=1.070, XH—C—H=125°. No unique angle is ob- 
tained from one moment of inertia. Finally, Donohue et aj. 
find 120°+8° for <HCH in spiropentane. The electron 
diffraction of 1,1-dichlorocyclopropane would determine 
this angle. 

In conclusion, it should be pointed out that the ideas of 
s,p hybridization and proportionality of energy to f? are 
approximations, and the agreement of these calculations 
with experiment is to a great extent fortuitous. For this 
reason, we did not attempt any further refinements such 
as simultaneous maximization of C—C and C—H bond 
energies. 


1 Research Associate for the American Petroleum Institute Research 
Project 44. 

2 National Research Fellow in chemistry. 

3G. H. Duffey, J. Chem. Phys. 14, 342 (1946). 

4L. Pauling, The Nature of the Chemical Bond (Cornell University 
Press, Ithaca, New York, 1940). 

5G. E. K. Branch and M. Calvin, Theory of Organic Chemistry 
(Prentice-Hall, Inc., New York, 1941), p. 279. 

6K. S. Pitzer, Science 101, 672 (1945). 

7M. T. Rogers and J. D. Roberts, J. Am. Chem. Soc. 68, 843 (1946). 

8 G. Herzberg, Infra-Red and Raman Spectra of Polyatomic Molecules 
(D. Van Nostrand, New York), p. 440. The C-—C distance given for 
120° is also incorrect. 

9J. Donohue, G. L. Humphrey, and U. Schomaker, J. Am. Chem. 
Soc. 67, 332 (1945). 





Electronic Shifts at Moderate Temperatures 
in Metals of Group VIII 


W. R. Ham* anp C. H. SAMANS** 
American Optical Company, Southbridge, Massachusetts 
May 31, 1946 


URTHER evidence has been secured, supplementing 

that offered some years ago for pure iron,! of what 
are believed to be specific electronic shifts in iron within 
ordinary temperature ranges. Similar series have also been 
established for all the other Group VIII metals by measure- 
ments of either the electrical resistance of pure or alloyed 
metal wires or of the electrolytic conductivity of glasses 
containing these metals. For all the metals, the tempera- 
tures of the discontinuities follow, within experimental 
limits, the simple Ritz formula: 


Ti, =ctl(1/ne) —(1/n*)], 


where ct is a constant which is different for each metal; 
no has the value 3 for Fe, Co, Ni, 4 for Ru, Rh, Pd, and 5 
for Os, Ir, Pt; and », for each metal, has the successive 
values mo+1, m0+2, +--+ &. 

No change was found in the critical temperatures either 
with concentration or with the form of chemical combina- 
tion. Hence the effect is believed to be a purely atomic 
phenomenon. 

The fact that the Curie magnetic transition tempera- 
tures of Fe, Co, and Ni, all occur exactly at some definite 
term in the appropriate series indicates that the data may 
be of great importance in clarifying the theory of ferro- 
magnetism. 

Further details will be published shortly. 

* Consulting physicist. 

** Research metallurgist. 


1C. B. Post and W. R. Ham, J. Chem. Phys. 5, 913-919 (1937). 
W. R. Ham and W. L. Rast, Trans. A.S,M. 26, 885-902 (1938). 
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On Steric Effects 


TERRELL L. HILL 
Department of Chemistry, University of Rochester, Rochester, New York 
June 19, 1946 


TERIC effects have long been recognized in both 

equilibrium and rate studies. Quantitative experi- 
mental data are now available in a number of instances.! 
However, on the theoretical side, it appears that only 
qualitative interpretations have been attempted. The 
purpose of this note is to suggest a very obvious and simple 
(but apparently new) semi-quantitative theoretical ap- 
proach to the problem which might be found useful at 
least in some cases. 

The forces involved in steric effects are well known: 
(1) groups or atoms (within a molecule) may repel each 
other if too close together (just as two gas molecules do 
for sufficiently small inter-molecular distances); (2) in 
order to decrease this interaction, the groups or atoms will 
tend to move apart, but this will generally require the 
stretching or bending of bonds with a related increase in 
energy. The final configuration will thus be the result of 
a compromise between the two types of force, and will 
be the configuration of minimum energy. The above simple 
argument (which can be generalized) applies when we may 
neglect entropy effects, all configurations but the one of 
lowest energy, polar effects, solvent effects, etc. 

The suggestion to be made here is that it may sometimes 
prove worthwhile to attempt to formulate the above dis- 
cussion quantitatively. For example, we might express 
approximately the steric energy, E*, of a molecule as 
follows: 


Yo 


ro 12 6 
z= ud) —2uo(“*) 442,(1—1o)2+4k,(0—6)2, (1) 


where additional terms are necessary if there is more than 
one pair of interacting groups, more than one bond being 
stretched, or more than one bond angle being bent. In 
Eq. (1), 1 is a bond distance, @ a bond angle, k; and ky are 
force constants, J) and 60 are the values of / and @ in the 
absence of steric effects, 7 is the distance between the two 
interacting groups, and —xo is the minimum value of the 
interaction energy (at r=ro). Equation (1) assumes, as an 
approximation, that the Lennard-Jones form of the inter- 
action potential is appropriate here and that Hooke’s law 
holds (for large displacements in some cases). 

The minimum value of E* in Eq. (1) is desired. This is 
found using r=r(I, 0), 9E*/al=0 and dE*/3@=0. The con- 
stants up and 79 may be estimated from van der Waals 
constants and radii, or other data. The force constants 
are hypothetical force constants. They are not the force 
constants which might be determined experimentally, but 
rather those which would exist in the absence of steric 
effects. Hence, they must be estimated by analogy, guessed, 
or calculated backwards from experimental data. 

Calculations of this kind might also be of interest from 
the point of view of studying (classically) polar and steric 
eflects on molecular structure.” 

The best type of equilibrium experimental data for 
Comparison with theoretical calculations is of the sort 
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obtained by Brown': measurement of AH, AF, and AS 
for a series of gas phase reactions in which there is a sys- 
tematic change in substituents. The procedure then is to 
calculate E* for the reactants and products of the different 
reactions in the series (generally E* may be taken as zero 
for some of these substances), and compare calculated 
values of AE;* (where 7 refers to the ith reaction in the 
series) with the deviations of AH; from the values of AH; 
which would be expected in the absence of steric effects 
(e.g., in the absence of steric effects, AH; might be a linear 
function of 7). A corresponding treatment of steric effects 
on reaction rates involves calculation of E* for the reactants 
and the activated complexes in the series. 

Sufficient preliminary equilibrium calculations have been 
made for the gas phase addition reaction between tri- 
methylboron and the methylamines (and ammonia) to 
indicate that forces of the sort discussed here can give 
effects of the type and order of magnitude found experi- 
mentally by Brown. However, a discussion of results will 
be deferred until time is available to carry out more de- 
tailed and refined calculations (e.g., including polar effects). 

Note added June 25, 1946: It has come to the author’s 
attention that F. H. Westheimer has been developing in- 
dependently a similar approach to this problem, with par- 
ticular reference to reaction rates. 

1 See, for example, H. C. Brown and co-workers, Science 103, 385 
(1946); J. Am. Chem. Soc. 67, 374, 378, 503, 1452, 1765 (1945). 


2L. Pauling, Nature of the Chemical Bond (Cornell University Press, 
Ithaca, New York, 1945), pp. 78-80. 





Nitrogen and Stearic Acid Adsorption by 
Supported and Unsupported Catalysts 


HERMAN E. RIEs, Jr., MARVIN F. L. JOHNSON, AND JOHN S. MELIK 


Research and Development Department, Sinclair Refining Company, 
East Chicago, Indiana 


June 21, 1946 


T has been suggested earlier! that diatomaceous earth 

supports of low area (and relatively large pores) may 
be responsible for the high area small pore structure of 
certain supported catalysts. Considerable evidence sup- 
porting this viewpoint and bearing on the structure of 
supported catalysts has recently been obtained in these 
laboratories by means of nitrogen adsorption at low tem- 
peratures? and stearic acid adsorption from solution.’ The 
following observations strongly indicate that the unsup- 
ported catalyst precipitated in the absence of diatomaceous 
earth forms finely-divided essentially non-porous particles, 
whereas the supported catalyst precipitated in the presence 
of the low area diatomaceous support develops a small pore 
structure of much greater area. 


UNSUPPORTED CATALYST 


The nitrogen and stearic acid areas of an unsupported 
catalyst, U,‘ are, respectively, 91 and 93 sq.m per g 
(ignited weight). This good agreement indicates that the 
large stearic acid molecules reach all the surface available 
to the smaller nitrogen molecules and that the particles 





466 LETTERS TO 


are effectively non-porous. With the great majority of 
similar materials studied in the range of this and higher 
areas the nitrogen values are greater than the stearic acid 
values by a factor of two or more. Furthermore the 
nitrogen adsorption-desorption isotherm for catalyst U 
in the 0.4 to 0.6 relative pressure region does not show the 
hysteresis effect characteristic of the high area supported 
catalyst and generally attributed to small pores. Finally 
the small non-porous particle size required for an area of 
91 sq.m per g is qualitatively confirmed by the extremely 
slow sedimentation rate in benzene observed in the stearic 
acid experiments. Relatively large pores are, of course, not 
excluded by the above observations. 


CATALYST SUPPORTED ON DIATOMACEOUS EARTH 


The nitrogen and stearic acid areas of the catalyst 
supported on diatomaceous earth, H—G,‘ are, respec- 
tively, 270 and 46 sq.m per g (ignited weight), or expressed 
more significantly as area per gram of catalyst material, 
approximately 516 and 90 sq. m per g (see Table I). The 
increase in nitrogen area over that of the unsupported 
material is more than fivefold. Since the larger stearic acid 
molecules are adsorbed to the same extent as on the un- 
supported catalyst they evidently reach practically none 
of the internal surface of the high area supported material. 
The development of fine pore structure in the supported 
catalyst is also demonstrated by the marked adsorption- 


TABLE I. Nitrogen and stearic acid adsorption by supported 
and unsupported catalysts. 








Nitrogen 
hysteresis 
(0.4 to 0.6 
Area, sq.m/g relative 
Nitrogen (BET) Stearicacid pressure) 


U (unsupported) 91 (91)* 93 (93) 
H —G (on diatomaceousearth) 270 (516) 46 (90) 
T (on titanium dioxide) 72 (134) 45 (82) 





None 
Marked 
None 








* Values in parentheses are areas per gram of catalyst material. 
Catalyst material constitutes approximately fifty percent by weight 
of the supported samples. In computing these values the maximum 
possible contribution of the support is subtracted (although a consider- 
able portion of the support surface is undoubtedly covered). The 
nitrogen and stearic acid areas of the diatomaceous earth are, respec- 
tively, 23 and 2 sq.m/g and for the titanium dioxide, 10 and 8 sq.m/g. 


desorption hysteresis effect observed with nitrogen in the 
0.4 to 0.5 relative pressure region.® 


CATALYST SUPPORTED ON TITANIUM DIOXIDE 


In sharp contrast to the results obtained with a diatoma- 
ceous earth support are those found for the same catalyst 
supported on non-porous titanium dioxide. The nitrogen 
and stearic areas of this supported catalyst, T,‘ are, respec- 
tively, 72 and 45 sq. m per g or 134 and 82 sq. m per g of 
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catalyst material. No hysteresis is observed in the 0.4 to 
0.6 relative pressure region of the nitrogen isotherm. 
Evidently the non-porous titanium dioxide support does 
not produce the high area fine pore structure obtained 
with the diatomaceous earth. 


1 Ries, Van Nordstrand, and Teter, Ind. Eng. Chem. 37, 310 (1945), 

2 Brunauer, Emmett, and Teller, J. Am. Chem. Soc. 60, 309 (1938). 

3 Harkins and Gans, J. Phys. Chem. 36, 86 (1932). 

4 The unsupported catalyst, U, was prepared by the reation between 
sodium carbonate and a salt of a heavy metal whose oxide is readily 
reducible. The precipitation of the supported catalyst, H—G, was 
similarly performed in the presence of suspended diatomaceous earth. 
In the preparation of catalyst T, titanium dioxide was substituted for 
the diatomaceous earth. All catalysts were calcined at 350°C for two 
hours. Catalyst material constitutes approximately fifty percent by 
weight of the supported samples. The catalysts were prepared by 
Dr. Olson of these laboratories. 

5 Ries, Van Nordstrand, Johnson, and Bauermeister, J. Am. Chem. 
Soc. 67, 1242 (1945). 





The Interpretation of Pressure Dependence of 
Accommodation Coefficients 


Joun L. Morrison 


Department of Chemistry, University of Alberta, 
Edmonton, Alberta, Canada 


June 24, 1946 


N a recent paper, Amdur! interprets the pressure de- 

pendence of the accommodation coefficients of several 
gases on platinum. He attributes the pressure dependence 
to simple Langmuir adsorption of the same gas as that of 
which the accommodation coefficient is being measured. 

Probably this interpretation is the correct one for the 
gases He, Do, No, CO, and Ov. However an alternative 
interpretation is possible for the inert gases He, Ne, A, 
Kr, and Xe. 

In the paper giving the experimental results, Amdur, 
Jones, and Pearlman? report using neon of 99.9 percent 
purity. Thus in 0.1-mm pressure neon, there is impurity 
of 10-*-mm pressure. 

Morrison and Roberts,* using accommodation coefficient 
measurements, found that a clean tungsten surface adsorbs 
oxygen at pressures as low as 10-* mm and becomes com- 
pletely covered with molecular oxygen at a pressure of 
about 10-> mm, always in the presence of neon at 0.1-mm 
pressure. Moreover, Roberts! found that clean tungsten 
adsorbs a saturated film of hydrogen at 4X10~*-mm pres- 
sure in the presence of neon at 0.07-mm pressure. 

Since the amount of impurity in each of Amdur’s inert 
gases is directly proportional to the pressure of the inert 
gas, the adsorption of this impurity could account for the 
fitting of the Langmuir isotherm. 

11, Amdur, o Chem. Phys. 14, 339 (1946). 

—-> . M. Jones, and H. Pearlman, J. Chem. Phys. 12, 159 


3 J. L. Morrison and J. K. Roberts, Proc. Roy. Soc. A173, 1 (1939). 
4J. K. Roberts, Proc. Roy. Soc. A152, 445 (1935). 





